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1. INTRODUCTION 


1.1 The Need for a General Analysis ' D 

Silicon solar cells have been used for 22 years and have been a 
major space power sourcb from the very beginning. During the last two 
decades there has been a great deal of improvement in the basic design 
and technology of solar cells and this has resulted in the latest 15% 
AMO cell [1] as compared to the 6% cell of 1954 [2]. Despite recent 
theoretical analyses, which point to a pra<jtical 19-20% AMO efficiency, 
there is still a technological "gap" in achieving this high efficiency 
t3,4]. The lack of agreement between theory and actual conversion effi- 
ciency has befsn the basic motivating factor in the development of a 
complete solar cell numerical analysis program. 

silicon technology has reached a very high degree of development 
allowing meaningful comparisons between theory and experiments. Hence 
the present work emphasises the correlation of theoretical and experi- 
mental data In addition to the development of a complete solar cell 
analysis. It is believed that through a detalled comparison. It may 
be possible to reveal Che problem area which could evenually lead to 
performance Improvements and high conversion efficiency. 

1.2 A Brief Review of the Development of Silicon Solar Cells 
Although the discovery of the photovoltaic effects in an electro- 
lytic cell was made by Becquerd In 1839, the first practical solar cell 
was not made until 1954 by Chapin, Fuller and Pearson [2] . In the 
following years there was progress in the understanding of 


solar cell theory such as the spectral response theory, the p-n /^unction 
theory, the series resistance effect and the determination of the 

optimum band^ap for materials, etc. At the same time, |solar cell 

•( 

efficiency was steadily Increased to about 10 percent. However the 
major technological and theoretical breakthroughs have only been achieved 
In this decade. These accomplishments should be attributed to the 
extensive financial support of photovoltaic research activities through- 
out the world since 1972. This section reviews the progress of the 
silicon solar cell during the past seven years. The earlier development 
the silicon solar cell has been summarized elsewhere. 

Considerable efforts have been given to raise the efficiency of 
silicon solar cells In the past few years. The average cell efficiency 
Is about 15 percent under AMO solar Intensity. Practical high 
efficiency cells of 19 to 20 percent are expected to be achieved In 
the near future. 

1*2.1 High Efficiency and Short Circuit Current 

The high efficiency of silicon solar cells has been partly achieved by 
increasing the output current density in the last few j^ears. This has 
been due to novel technologies for increasing the surface collecting 
efficiency and Improving violet photon response. In 1973, a new 
generation of silicon solar cells was produced which utilized Ta 20 ^ oxide, 
shallow junctions and a fine grid pattern to enhance the collection 
efficiency . The metal coverage area has also been reduced from .10 per- 
cent to about 6 percent and the number of collecting fingers has been 
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increased from 3/cm to 30/cm. In 1975 a textured cell „as also fabri- 
cated and the surface reflectivity was reduced to the lowest limit of 
about 3 percent. Figure 1.1 shows the chronological progress of the 

average short circuit current density of cells ^ ^ 

ensicy ot cells produced during the past 

seven years. 

Efficiency and Open Circuit 

Efforts to ..mprove cell output voltage has been steadily going 
on during the past few years. These Include the trends of using a low 
resistivity substrate of 1.3 «-cm instead of the conventional 10 O-cm 
substrate and the fabrication of a high low bach surface function for a 
resistivity cell, in 1973. the development of this back surface 

field or BSP cell produced about a 50 mV Increase 1„ o . 

w mv increase in open circuit voltage 

for a thin 10 fi-cm solar cell. However, the highest voltage achieved 

from low resistivity cells is much less than that obtained from simple 

theoretical predictions. Hence, there is recently a great interest in 

the physical explanations of this discrepancy between theory and 
experiment. 

Hfgh Effi ciency and Curve 

The conventional 10 0-cm silicon solar cell has a low curve factor 
of 0.72. This value is much less than the theoretically predicted value 
of 0.82. The discrepancy has been identified as due to high space 
charge current density and high series resistance of a high resistivity 
cell. However, the curve factor of a recent 1.3 f!-cm cell has been 




by a lower series resistance and lower saturation current density of a 
low resistivity cell. Further Improvements of the curve factor can 
be obtained by using a low resistivity substrate and an Improved design 
of the collection junction and grid pattern, 

( - . • .7 

1,3" The Major Characteristics of the Analytical Technique 
The approach taken by the present analysis of solar cell operation 
Is a complete numerical solution of the semiconductor device equations. 
This system of equations Is quite well known,, although it has been only 
recently that solutions have been possible due to the advent of high 
speed computers and computer oriented numerical techniques. These 
equations include Poisson’s equation, the hole and electron current 
density equations, and the hole and electron continuity equations. Much 
Information regarding these equations is presented in the following 
chapters of this report; however , | some general remarks can be made as 

J - 

to the generality and completeness of the operation which they describe. 
It is known that the presence and interplay of both drift 

i 

and diffusion, current components must be present in the 
solution to allow for the appearance of such effects as high 
injection and bulk resistance. This is possible through the use of the 
general current density equations. The continuity equations include an 
internal net recombination rate which plays a significant role in solar 
cell operation. This phenomena is readily included through a selection 
of a recombination model with empirical models for lifetime parameters 
based on available experimental data. In a solar cell there is the 
additional factor of an external generation rate due to the full spectrum 
solar irradlance. This highly variable spatial factor is calculated 
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separately aed then Included point by point In the continuity equation 
through the external genetatfion term. Thus the analysis does not start 
with an assumed optically generated short circuit current, but begins 
with the generation rate Itself at each point within the device. This 
term has been calculated using the most current information available 
on spectral irradlance for various air mass conditions, antireflection 
films, and the relevant material properties for silicon. The inclusion 
of this term directly into the continuity equation avoids any 

a £riori assumptions with respect to the collection efficiency in various 
regions of the device. 


The surface region of contemporary solar cells is analytically 
quite complex. These regions are typically very heavily doped, 
diffused regions. The analysis Includes an erfc or Gaussian It^urity 
profile and the attendant electric fields due to this type of profile 
are ^ thus included m the analysis through a spatially dependent Impurity 
concentration in Poisson's equation. The doping of the surface region 
is usually high enough to cause significant heavy doping effects. Thus 

various models postulated for heavy doping effects are Included within ' 
the solution. 

Within the device equations themselves, there are several parameters 
such as mobility, lifetime, and diffusion coefficients which in a real 
device are not constant but vary with doping level and/or electro- 
static potential variation. These variations must be Included in an 
accurate analysis. Although not derived analytically, these 
variations are taken as empirical data from currently available 


measurements . 



. in general rhe analysis presented is ,„Ue free fros. the normal 
Simplifying approKinations s.de in semieonductor device analysis. A 
summary of the najor featnres of the analysis include the following: 

(a) Optical generation is calculated directly from available 
empirical measurements of a specific irradlance spectrum. Including 
the effects of antireflection films and wavelength dependent 
absorption and index of refraction coefficients. 

(b) I'he optically generated current is calculated directly from 
the interplay of the above generation rate and the device operation 
l.e. there are no assumptions pertaining to collection efficiency. 

(c) Recombination is included within the analysis, not only for 
the bulR regions, but also for surface and depletion regions. 

W drift and diffusion components of current flow are 

included. This allows the appearance of high injection effects, 
resistive loading, and the effects of any Dember type potentiall ^ 

(S) A diffused impurity profne^ is included in the surface region 
through an °r more complex impurity distribution. 

(f) Band gap shrinkage due to heavy doping effects is included 
In the diffused surface region. 

(g) A non-ohmic contact is Included at the irradiated surface 

through a finite surface recombination velocity. 
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1.4 General Device Equations 

The fundamental physical mechanisms lii semiconductor device 
analysis such as the existence of space charge regions, drift and 
diffusion currents, and carrier recombination are contained in the 
general device equations. In one dimension, neglecting any magnetic 
or thermal gradients, these can be expressed as: 


ll.-a [p - n + N(x)], 


E 



> 3 n 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 

, (1.4) 


3n 
3 t 


U + G + 
e 


1 

q 



3t 


U + G - 


, 3J 
i_JE. 
q 3x 


(1.5) 

( 1 . 6 ) 


These equations have general three dimensional forms; however, in this 
work only the one dimensional case has been considered. In the above 
equations, U is the Internal net recombination-generation rate, is 
any generation rate due to external physical processes, and the other 
terms take on their conventional meanings. The net ionized Impurity 

I 

doping is represented by N(x) and can be a complicated function of x 
! : 
in that it can represent the doping profile in a diffused region plus 

the changeover from donor to acceptor doping as an n-p interface is 

crossed. Equations (1.1) and (1.2) are Poisson's equation and the 
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defining relationship between electrostatic potential « and electric 
flild E. The current densities are expressed by Equations (1.3) and 
(1.4) and include both drift and diffusion terms for both species of 
carriers. Equations (1.5) and (1.6) are the continuity equations for 
holes and electrons. The InterMl net recombination rate of electron- 
hole pairs due to thermal processes as represented by U. is treated 
in greater detail in the following sections. Chapter 3 discusses Che 
extensive analysis behind the development of the electron-hole pair 
generation due to external soiirces such as Incident as 

represented by the term G . 

e ; , . 

In addition to the above equations several connecting expressions 
can be notg. Fbr one, the Einstein rilations can be used in relating 
the diffusion coefficient to mobility. For non-degenerate doping 
these can be expressed as: 


Another set of useful relationships can be obtained through a definition 
bf quasl-Permi potentials. In equikbrl™ the Fermi energy, E^, can be 
defined from Boltzman statistics as follows: 


n ' exp 

E -E 

P ' Hi exp [-^3, 


( 1 . 8 ) 


/ 
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where is the energy level corresponding to the center of the band 
gap. I The electrostatic potential nay be expressed as - -E^/q and the 
Fermi potential as ^ - E^/q. In terms of these potentials. Equation 

(1.8) may be expressed as 


n - n^ exp ] , 

P = Pi exp . 


(1.9) 


In non-equilibrium conditions, the quasi-Fermi potentials, and ^ , 

n p 

may then be defined from 


n “ n^j^ exp [- 


q(^-<l'„) 


p = n^ exp [- 




Note that under these conditions the np product is no longer necessarily 


equal to n^ ; in fact, 


np = n. exp [ 


2 9(<l'-'l'„) 

^ r p n . 


( 1 . 11 ) 


This defines a quasi-equilibrium as that condition in which the np 
product remains constant although differing from its equilibrium value. 
Thus the difference In quasi-Fermi levels can be viewed as a measure of 
the amount of displacement from equilibrium. 

The general set of six equations and six unknowns as expressed by 
Equations (1.1) thru (1.6) then represent all necessary phenomena for 



5 Plar cell analysis excluding of cpuirse varjUtlons in noblltty, choice 
of recomb^tlon ||)del, and weroal generation rate. * For the complete 
analysis however this set of ^ equations can be reduced to a set of 


three «o“pled eqijatlpna tlitpHgh, the quae leyel and electric 


1.5 Organization of the Report 

this report is divided into several major chapters. The first 
chapter Is a simple Introduction of a silicon solar cell and its 
current status of technology. The partlcdlar semiconductor transport 
equations which are of fundamental Importance In solar cell type devices 
are also discussed. Chapter 2 presents the results of the detailed 
correlation between theory and experiment. The detailed comparisons 
have Included the most important measurements such as photovoltaic I-V 
xharacterlstics, dark I-v characteristics and spectral responses, etc. 
The simuUtlons of violet and CNR cells are pfesented In Chapter 3. 
Chapter 4 discusses the temperature dependence of the short circuit 
current density. It Is emphasized In this chapter that the temperature 
dependence of the absorption coefficient cannot be neglected In the 
calculation of short circuit current. Chapter 5 discusses the 
heavy doping effects and Keldysh-Franz effects In order 


to explain the low open circuit voltage in, a highly doped cell. The 
characteristics of a | thin kllcon solar cell are presented in Chapter 6. 
The effects of the series resistance on concentrating solar cells are 


shown In Chapter 7. Due to the distortions of the Input spectral 
intensities of a concentrator, the effects of non-uniform illumination 



in llniitlng solar Cell conversion efficiency are discussed in Chapter 8. 
Chapter 9 presents calculations for a recently proposed hlgh'^low 
junction emitter solar cell. Finally, appendices 10 is divided into 
five sections. Appendix 10.1 discusses the detailed comparlsoh of 
theoretical and experimental solar cell performance. Appendix 10.2 
presents a two dimensional analysis of sheet resistance and contact 
resistance effects in solar cells. Appendix 10.3 discusses impurity 
gradients and high efficiency solar cells. A detailed two-dimensional 
model of a solar cell is presented in Appendix 10.4. Finally, 

Appendix 10.5 discusses the fill factor and diode factor of a solar cell. 
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2, BASIC SOLAR CELL CHARACTERISTICS AND DEVICE PARAMETERS 


Since the discovery of metal-semiconductor non-ohmlc behavior the 
rectifying effects between metal-semiconductor and semiconductor- 
semiconductor contacts have received a great deal of attention. In 1949 
Shockley proposed the modern p-n junction theory which established the 
Important role of minority carrier density and Its exponential behavior 
across the junction barrier [1]. The transport equations for minority 
carriers are particularly simple for low injection and uniform <ioped 

semiconductor regions. The minority carrier current density can be 
expressed as 


J - Jglexp(qV/nkT)-l], 


n 




n P 

where the saturation current density is a function of semiconductor 
parameters on both sides of the junction, and n » 1. Departures from 
Shockley’s simple I-V characteristic are usually observed in silicon 
at room temperature and further evolutions of the p-n junction theory v 
have modified and extended Shockley’s theory [2-5], 

In all silicon p-n junctions, several current transport mechanisms 
may exist simultaneously. The diffusion current density which is due 
to the injection of minority carrier over the junction barrier is, 
of course, the most important. Other mechanisms include recombination 
current within the depletion region [2], tunneling through the band- 
gap for highly doped semiconductors [5] and high injection effects for 
high resistivity semiconductors at large forward bias voltages [3,4]. 
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In addition to these fundamental limitations there are several other 
current mechanisms which are due to Improper fabrication processes 
and/or material Imperfections. Especially Important are the series 
resistance and any shunting resistance. Fortunately these can be 

minimised by using good contact metals, grid patterns and proper 
Sintering treatments. 

In general, the departures from the simple diffusion theory always 

lead to poor rectification In diodes and poor curve factors and low 

open circuit voltages in solar cell application. Figure 2.1 shows the 

Ideal I-v characteristics and some of the modifications at forward 

bias, voltages. As can be seen In the figure, the simple Shockley 

diffusion current has a diode n factor of 1 for all bias voltages. The 

space charge recombination current has an n factor of 2. Such an n 

factor may also be found at high currents caused.by high Injection. 

An abnormally large n value may be found in some devices at small 

voltages which Is caused by small shunting resistances. An n factor 

of 2 may also be present at high voltages and caused by a high sheet 

resistance. Curves (c) and (e) of Figure 2.1 show examples of these 
effects. 

The dark 1-V characteristics of a solar cell are as Important as 
the short circuit photocurrent In determining the efficiency and power 
output. The components of the dark I-v characteristics described above 
are discussed in detail In the following sections. 



CURRENT- DENSITY Cfl/2 CM^) 


KT* 


i(r» 


i(r* 


i(r* 


icr* 


10 -' 


a. Simple diffusion current 

b. Space charge current 
Shunting leakage current 

d. Modified by contact resistance 

e. High injection or modified 
surface sheet resistance 

f. I'^odified by combined contact 
and sheet resistance. 


1 o-“| J ( 1 L I I I H 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

iVOLTHGE (VOLTS) 

figm^ 2.1. Prototype of the dark I-V charaoteristioa of a solar cell. 
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2.1 Simple Diffusion Current 

The current density for minority carriers in Shockley’s model is 

(f 

J “ Jg[«?p(qV/kT)-l], (2. 

where the sayuratien current density is a function of aemlconductor 
parameters and the appropriate boundary conditions. First order 
analytical expressions for havebeen compiled in Hovel's book [6] for 
several models with different boimdary conditions. T"*" ' 

2.2 Space Charge Recombination Current 

I 

! The generation-recombination current of the Sah-Noyce-Shockley 
model is given by Equation 2.3 which assumes that the recombination is 
through a center located in the vicinity of the center of the bandgap. 

pn-n . ^ 

J = q f ^ j 

Y o T (n+n )+T (p+p,) U) 

po 1' no'*^ ‘^ 1 ' 

2 

qn^ W^[exp(qV/kT)-l] 

^poV'"noPl-^+^'^po+'^no^"i®''P<'5V/2kT) » <2.3) 

X +T ®^p(qV/2kT) for medium voltages • 
po no 

It is clear that for a silicon p-n junction the space charge 

recombination current has a diode n factor of two at medium voltages 
and at room temperature. 

2.3 High Injection Current 

High injection occurs when the minority carrier density on one side 
of the junction becomes comparable with Che majority carrier density, 
e calculation of the high injection current Indicates an exp(qV/2kT) 
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beljavior 17], H^gh injection is likely to occur for low base doping 
densities near the junction or fx>r silicon solar cells operating in 

1 . 

multi-sun environments. For normal silicon solar cells with 
^®®^®^^^^ties in the range of 10 0*cm to 0.1 0*cm which operate under 
one sun power intensity, high injection is unlikely to ocdhr. 

2.4 Tunneling Current t 

, i , ' 

A tunneling current may exist in heavily doped junctions with a 

resistivity of less than 0.01 0*cm. The tunneling current takes the 
form [8] 

“ KN^ exp(ev) (2.4) 

where is the density of energy states available fqr an electron or 
hole to tunnel into, and K and 8 are functions of semiconductor 
parameters. The n factor for tunneling currents lies between 1.3 and 

2 at room temperature [8]. 

2.5 Leakage Current 

I ' : 'V 

Since a solar cell is a relatively large area device, there is a 
great chance of a leakage channel existing through the imperfect 
junction, especially under the metal contact [9] . The leakage current 
can be modeled by a shunting resistor Rg^ across the junction and the 
current form is quite simple 
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Where is the junction voltage at the impe^jfectlon location. As a 
result of this leaky current, the diode n factor may he very high with 
values of 3 to 5 usually being observed at voltages of less than 0*4 yolts 

2.6 Current Voltage Characteristics Modified by Series Resistance 
Series resistance becomes important as the current density " 
increases and/or junction depth decreases. The series resistance comes 
from two sources: the surface sheet resistance and the metdl- semicon- 

ductor contact resistance. For the contact resistance which appears 
in series with the cell, the exponential dependence of current on voltage 
can be modified by replacing V in the exponential with V-R^I. For the 
sheet resistance, however, the two dimensional distributed nature of 
the current flow does not allow one to define a purely lumped resistance. 

In this case at large currents the equation becomest 

I =/ATl^exp(qV/2kT), (2.6) 

I - 1 

I T *' q ps^ * A* (2.7) 

where A is the total area of the solar cell, J is the saturation 

j : S 

current density cf the simple diffusion theory, is the surface 
sheet resistance and h Is the total perimeter of the contact grid 


pattern. , 

The parameter has the physical significance that it is the 
current level at which the characteristic makes a transition from an 
exp(qV/kT) dependence to an exp(qV/2kT) dependence. In a practical 


This effect is very similar to current crowding which occurs In 
bipolar transistors. A discussion is contaiued in Appendix 10.2. 
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silicon solar cell both sheet resistance and contact resistance may 
e^lst simultaneously and the diode n factor may be as high as 5 at 
"iroltages greater than 0.5 volts. In this combined case one cannot 


tiwdel the cell correctly by a lumped constant resistance. 


ir/.v, . ■ ii 


2.7 Parameters for Device Model 

r ' ' 

Diffusion Length and Lifetime. 

The lifetimes of electrons and holes are of great-’ Importance In 
understanding the electrical and optical behavior of a semiconductor 
device. For indirect bandgap semiconductors such as silicon » the 
carrier lifetime Is generally high and basically determined by 
recombination through intermediate centers within the bandgap instead 
of direct band-to-band recombination. The minority carrier lifetime 
has been developed by assuming a single Shockley-Read center as [10 ] 


[(1 + ^ exp[-(E(,-E^)/kT]) + exp[-(Eg-E^)fcT]] , 


P cr, 


p th R 


n no 


( 2 . 8 ) 


v*xere Xp ® hole lifetime in an n-type semiconductor with doping density 

^no* ^R ** recombination centers, a , a = hole and electron 

p n 

capture coefficients. A similar form can be written for electrons 
in p-type material. Equation 2.8 indicates that the minority carrier 
lifetime is lower in general for a higher doping density. This behavior 
of minority carrier lifetime with doping density has been experimentally 
observed [11,12]," 


Some representative curves of measured diffusion lengths as a 
function of doping density are shovra in Figure 2.2 [li]. Comparisons 


DIFFUSION- „ LEN8TM («m) 
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Figure 2.2. Model of the diffusion length and comparison to the experimental 
data. All experimental data are for p-type material ekcept 
thpse marked by the triangle. 
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k : 


r> >-> 


to Other experimental data are also shown on the same graph [11-15), 
At high doping densities band^to-band Auger recombination may become 
the dominant recombination process. This gives a decreasing life- 
time which Is Inversely proportional to the square of the doping 
density. The experimental curve of Ljj(MAX) In Figure 2.2 has a llfe- 
^^^^^^^^tlme dependence on doping density with an exponent of -1.7 at doping 


I 

t 


densities greater than 10^^/cm^. This value Is close to the 
theoretical band-to-band Auger lifetime model which has an exponent 

Hence It Is a fairly good approximation to define Ljj(MAX) and 
Ljj(MED) of Figure 2.2 as the upper and lower bounds for electron 
diffusion length In p-type silicon. Similarly Ljj(MED) and Ljj(MIN) 
of Figure 3.1 could he considered as upper and lower bounds for hole 
In n-type silicon because of the lower hole mobility. 
In a practical silicon solar cell, the actual diffusion length may 
vary between some upper and lower bound depending on the material 
perfection and the fabrication processes. In a solar cell the density 
of recombination canters is generally much smaller than the doping 

majority carrier lifetime equals that of the 
minority carriers [l6]. 


Although the diffusion length data of Figure 2.2 was measured in 
bulk material, it is assumed valid for the shallow diffused layer of 
solar cells. If the lifetime is a function of total doping density 
only, such as in the Auger process, this will be a good approximation. 
However, this may not be valid if lifetime is dominated by deep level 
Impurity recombination. 
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2 . 7.2 Surface Recombination Velocity 

' (/ 

In addition to bulk recombination, aurfaoe recombination le 

another loss mechanism which Is modeled by a surface recombination 

velocity S. The minority carrier current flow toward the surface Is 
given as 


where 4p^ is the surface excess minority carrier density. 

Ihe value of S is basically determined by surface conditions such 

as the density of Interface states, any antl-reflectlon oxide layer and 

surface treatments. Very low S values of 10^ cm/sec can probably be 

achieved only through the use of high temperature oxidation processes 

which may cause a drlve-ln of the surface diffused layer and may not 

be^compatlble with solar cell technology ( 17 ) . Hence, a value of 

10 ^ cm/sec may be a lower limit for SRV of typical oxide coated solar 
cells [ 18 , 19 ]. 

2 . 7.3 Diffusion Doping Profile ' 

It has been found that shallow diffusions (< 1 of phosphorous 
In silicon result in considerable deviations from the simple diffusion 
theory of an erfc function [ 2 o. 2 l). For short diffusion times 
(< 1 hr) and temperature below 1100»c. It has been found that a constant 
concentration layer exists near the surface of about 1/3 v 1/4 of the 
junction depth and the electrical active phosphorous concentration In 
this layer Is about half of the solid solubUlty limit at the particular 
diffusion temperature. Beyond the constant concentration region, the 
diffusion profile can be reasonably well represented by a complementary 
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error function. One of the typical diffusion profiles is reproduced 

here in Figure 2,3 [21], The parameters which characterize this 

'particular diffusion profile are the surface doping density C , width 

s 

of the constant doping layer X and the doping density C_ at the 

c B 

boundary of the constant doping and the erfc doping profile. 


2.7,4 Mobility 

Two major contributions to mobility are phonon scattering and 
impurity scattering. These effects make mobilities a function of 
doping density, temperature and internal electric field intensity. 
The general empirical equation developed by Gummel [22] was used in 
this work. 


0.2 


(~) 
4 


1 + 


I (E/A)^ 

N E/A + F 




( 2 . 10 ) 


f 

This equation has been confirmed by measuring the relations between 
drift velocity and electric field [2 3], 


2.7.5 Heavy Doping Effects 

Heavy doping phenomena occurs in silicon for total doping densities 
19 3 

above about 10 /cm . The high doping effect on minority carriers can 
be represented by a bandgap reduction where the empirical expression of 
Equation 2,11 has been used in this work. 

AEg = -0.45y/N/(10^^cm~^) (eV), (2.11) 

In this expression, N is the net doping density. 
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2.7.6 Absorption Coefficient 

The ability of a semiconductor to absorb light of a given wave- 
length Is characterized by the absorption coefficient a. Values of the 
absorption constant a were taken from Dash and Newman [24] up to about 
0.95 pm wavelength. At wavelengths above 0.95 pm there Is considerable 
variation In the reported absorption coefficient values. Several 
reported values at 0.95 pm, 1.0 pm and 1.1 pm are shown In Table 2.1 
[24429]. The data of Dash and Newman is seen to be larger than most 
of the reported data at 0.95 pm and 1.0 pm. Good agreement In the long 
wavelength spectral response of solar cells could not be obtained by 
using the data of Dash and Newman. Best results have been obtained 
by using the values of the last line Which are Intermediate between the 
largest and smallest reported experimental values. Thus the absorption 
data which has been used is that of Dash and Newman with the data at 
0.95 pm, 1.0 pm and 1.1 pm modified to the values of Table 2.1. The 
Index of refraction as a function of Incident wavelength was taken from 
Phillip (1972) [30]. 


Table 2.1 Comparison of reported absorption coefficient values at long 
wavelengths 




1.1 pm 

1.0 pm 

0.95 pm 



Absorption Coefficient 

(cm ) 

[24] 

Dash & Newman 

7 

100 

220 

[25] 

Vedam 

- 

- 

270 

[26] 

Runyan 

- 

67 

170 

[27] 

Vol'fson & Subashiev 

- 

64 

150 

[28] 

Macfarlon 

3.9 

61 

- 

This work 

3.9 

74 

204 
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2* 7* 7 Spectral Response 

I The spectral sensitivity of a solar cell to incident photons is 
measured by the spectral response or the quantum yield. For a practical 
solar cell the quantum yield is always less than unity because of 
surface reflection losses and internal recombination losses. 

Internal quantum yield can be defined as the ratio of the collected 
short circuit current density to the input current density which is 
generated by the incident photons assuming 100 percent transmission 
through the surface^ l.e. 


QY(X) 




qF(X)(l-R(X)-Aa)) (l-e“''<f‘ * 


( 2 . 12 ) 


where F(X) « incident photon flux (proportional to input power density), 
R(X) “ reflection at surface, 

A(X) “ absorption in AR layer if any, 

Wj “ device thickness. 


Another practical parameter is the external quantum yield QY (X) ,^hlch 
includes losses due to surface reflection and antireflecting layer 
absorption: 


QY ,(X) 


sc 


(2.13) 


^ qF(X)(l-e-V(^>) 

The spectral response is represented by the ratio of collected current to 
input power density as: 


SR _(A) „ 

ext ' P 




input 


(2.14) 
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The spectral response theory of •Pirince and Wolf [31] shows that the 
overall spectral responses can be considered as made up of somewhat 
independent responses from the surface and base layers. Hence It is 
sometimes useful to specify the spectral response from the surface 
region, depletion region and base region respectively as 

SR^^^(X) ■ SR^^^(X, surface) + SR^^^(X, depletion) + SR^^^(X,base) . (2.15) 

Some of the parameters and results of the spectral response analysis are 
shown In Table 2.2. 

The calculated reflectance R in Table 2.2 is in direct agreement 
with Phillips data of oxide free silicon [30] , although it is well 

O . .... 

known that a thin layer of oxide of about 20 'v« 35A in thickness may be 
grown on an exposed bare silicon surface. The correction on R due to 
such a layer is less than 1 percent for photon wavelengths of 0.4 to 
1.'0 micrometer (This also agrees with Ref . [30].). 

For Tantalum oxide calculations, a reflection index of 2.20 was 
used which is based upon ellipsometry measurements performed at a wave- 

O . 

length of 5461A [32]. This value of reflection index is in general 
agreement with reported literature values [33], 

The calculated transmission and reflection coefficients for Ta20^ 
are shown in Figure 2.4.. The data indicates a much better surface effi- 
ciency at short .wavelengths as compared to a bare. Si, SiO or Si02 coated 
surface. Since the AMO power spectrum peaks between 0.4 ym and 0.6 ym, 
Ta«0_ is superior to the other oxides studied. 
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Table 2.2. Parameters for spectral res ponse calculation at various wavelengths. 
0,4 0.45 0.5 0.6 0.7 0.8 0,9 0.45 1 


Absorption Coeff. 
(cm“^) 

8.70E?'* 

2.62E4 

1.23E4 

4.56E3 

2.10E3 

9.64E2 

3,67E2 

2.04E2 

7.42E1 

Photon Energy 
(eV) 

3.09 

2.75 

2.47 

2.06 

1.77 

1.54 

1.37 

1.30 

1.24 , 

Transmission for 
Bare Si 

0.521 

0.583 

0.615 

0.647 

0.663 

0.672 

0.679 

0-pBl 

f-'- 

4 % ' 

0.683 

Reflection for 
Bare Si 

0.478 

0.416 

0.384 

0.352 

0.336 

0.327 

0.318 

0.318 

0;316 > 

INPUT POWER ■ 
- (mW)' 

10 

10 

10 

10 

10 ■ 

10 

10 

10 

10 

Surface Rate 

9.15E20 

3.47E20 

1.91E20 

8.91E19 

4.91E19 

2.61E19 

1.13E19 

6.85E18 

2.55E18 " 

INCIDENT CURRENT 
(mA/cm^) 

3.226 

3.629 

4.032 

4.839 

5.646 

6.452 

7.259 

7.662 

8.065 
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SI + 5S5A Tq^O^ j / 

I / // 

I / //Si4IIOOi 


Si4 600A SiO 


liOOA SiOs 


WAVELENGTH (uM) 

Figure 2i4. Transmission coefficients for bare and antl-reflectlon film 
coated silicon surfaces at wavelength 0.35 pm to 1.1 pm. 
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8 General Comparison of Theore^.ical and Experimental Results 
\ The solar cells which have been studied can be broadly divided into 
th^^ categories based .upon the origin of the bulk material and the cell 
design. The first category consists of n’^p cells on 0.1 n*cm base layers 
''ith a finished thickness of about 6 mils. The second category consists | 
of n p cells on 10 cm base layers with thicknesses of about 10.5 mils. ' 
Both types of cells were made in a standard 2 cm^ area and use a NASA 

'• ' 1 ■ 'A 

Lewis Re^arch Center 10-finger grid. No anti-reflection layers were 
present on the silicon surface* 

. t "■ 

The third type of cell is the Aluminum BSF cell on 16 J2*cm substrates 

with a finished cell thickness of about 6.5 mils. These cells have a 

O 

TajOj coating about 595A In thickness and a 5 mil "Teflon" FEP cover on the 


Tantalum Pentoxide. On these the nine finger grid pattern of Spectre 
Lab was present. 

The n-type surface layers were phosphorous diffusion, using POCl^ at 
the NASA Lewis Research Center. The temperature and the duration of the 


diffusion process are described for each type of cell in the following 


sections. The top and bottom contacts were made using metal masks and 
by evaporating a thin layer of Aluminum (200'V/500A) followed by the 
evaporation of about 3 to 5 micrometers of silver on the surface. The 
contacts were then sintered at temperatures of 550 to 650®C in H^. 

Al-Ag contacts have been found to have less degradation in the cell 
electrical characteristics than that which occurs for sintered Ag-Ti 
contacts, although Ti makes a better ohmic contact than Al. For BSF 



ii 

; : ■=> 92 

cells the Aluirilnum was alloyed at about 800®C for one hour or less and 

thi^ produced a high-low Junction of 0.5 to 1.0 \m in depth [34]. 

The detailed comparison of the theoretical and experimental results are 
shown as Appendix 10.1 and 10.2. 
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3. HIGH EBTICIENCY SILICON SOLAR CELLS 


3.1 Introduction 

Although silicon solar cells were first wade during the mid- 
fifties, several Important technological breakthroughs were not 
achieved until recent years. During this time. Improvements In solar 
cell -.technology have lead to the production of high efficiency 
silicon solar cells. Higher efficiency of a silicon solar cell has 

been mainly achieved by increasing the output voltage and/or short 
Circuit current density. 

The discovery of the back surface field concept raised the output 
wltage of a lightly doped solar cell m 1973 [1]. Then the violet 
cell, which utilises a very shallow Junction, fine grid and Improved 
antl-reflectlon film to enhance the optical responses of short wave- 
length photons, was introduced to produce higher short circuit current 
density [2]. later this cell was further Improved by preferentially 
etching to give a serrated surface and higher surface efficiency was 
achieved. This novel cell was named a CNR (COMSAT non-reflective) cell 
due to the non-reflecting properties of the serrated surface [3] . 
Calculations of a CHS cell are difficult because the Incident light 
la not in a normal direction to the cell surface and several reflections 
and refractions usually occur on the serrated surface. 

It is the purpose of this section to model these high efficiency 
solar cells Including the violet and CNR cells, it Is also shown that 
the limitation td cell efficiency due to the surface optical efficiency 
can be predicted from the calculations. 


3. 2 Optical Reflection and tranaioission 

Traditional solar cells have a flat polished surface and normal 
illumination can be assumed. In this case the surface reflection and 
transmission can be straightforwardly calculated. However, the 
calculation for oblique, illumination is far more difficult. A brief 
.review Is outlined In Appendix (A). 

In a CNR cell, light illuminates normally at the serrated surface 
where it strikes the face of a pyramid at an angle of 54.75”. The 
reflected rays then strike the face of an adjacent pyramid at a smaller 
angle of 15.75” as shown in Figure 3.1. it is known that the reflected 
light will be elllptlcally polarized whenever the Incident angle is' 
other than the normal direction. Also, the Fresnel reflection 
coefficient is different for TE or TM polarized light. Therefore, the 
total surface reflection and transmission coefficient must be deduced 
by the superpositions of the decomposed TE and TM components of the 
incident light which is assumed here to be uniformly polarized. 

Figure 3.2 shows the surface transmission coefficients of the 
incident light with angles of 54.75” and 15.75” for TE and TM modes 
respectively. The refractive Indea of silicon is taken from reference ■ 
[4] and is Shown in Figure 3.3. The calculated surface transmission 
coefficients of the violet and CNR cell with and without antl-refleciion 
fllmsare shown in Figure 3.4 and Figure 3.5 along with reported 
experimental values. The comparison to the reported values is very good 
which substantiates the above calculations [3.5]. Although the compari- 
sons are not satisfactory for wavelengths greater than 1.0 pm, there 
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are several other reported values which do follow the calculations In 
these wavelength ranges [6-8]. Several calculated surface optical 
properties are shown In Table 3.1 and Table 3.2. 



Figure 3.1. The serrated surface of a CNR cell. 

O 

The optimum thicknesses of Ta20^ are found to be about 600A and 

O 

700A for a solar cell under AMO and AM2 solar spectrums respectively., 

A solar cell under an AM2 solar spectrum has a higher optimum thickness 
of anti-reflection films. This is because the AM2 spectral intensity 
peaks at higher wavelengths. The average surface reflection losses 

O 

weighted to the AMO spectral density are about 3.5 percent for a 600A 
thickness of T&2^y This value is very close to the reported value of 
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WAVELENGTH (uM) 

Figure 3,2. The Qbliq^ue iurface transmission Coefficient at a Bare 
Silicon Surface. 
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WAVELENGTH (uM) 


Figure 3,4. 


Transmission coefficient for a 
wavelength 0.3 urn to 1.1 uW, 


plane silicon surface at 
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Figure 3.5. Transmission coefficient for a serrated silicon surface 
at wavelength 0.3 uitt to 1.1 um. 
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o optical Generation 

/A\ 



2.5 tc 3 percent [31. It eheuld be noted that the above option, thick- 
ness of anti-reflection oalde refers to the perfect collection of the 
input light intensity e 
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Table 3.2. Optical properties of CNR cell with an ant i'-ref lection film 
of Ta 20 ^. 


Spectral 

Condition 

Thickness 
of o 
Ta20j (A) 

Surface 

Efficiency 

(X) 

Available 

Optical 

Current 

(mA/cm^) 

Surface 

Generation 

Wavelength 
of peak 
Transmission 
(pm) 

Si AMQ 

0 

86.35 

46.40 

2.93E22 

1.05 

AMO 

500 

96.18 

51.68 

3.62E22 

0.45 

AMO 

600 

96.48 

51.84 

3.29E22 

0.53 

AMO 

650 

96.46 

51.83 

3.24E22 

0.57 

AMO 

700 

96.36 

51.78 

3.26E22 

0.61 

Si AM2 

0 

87.84 

30.13 

3.84E21 

1.05 

AM2 

500 

^ 96.98 

33.26 

5.00E21 

0.45 

AM2 

600 

97.93 

33.59 

4.58E21 

0.53 

AM2 

700 

98.16 

33.67 

4.27E21 

0.61 

AM2 

800 

97.80 

33.54 

4.20E21 

0.71 

AM2 

900 

97.05 

33.29 

4.40E21 

0.80 
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as 


3.3 Optical Generation 
In a solar cell, the number ofabsorbk photons can be 


Va<*> ■ '*phl<^>'l - e*P(-<>a)x)] 


represented 


(3.1) 


*’phi^^> 'he Input photon, density at wavelength X within the 
wavelength Interval AA and e(A) is the absorption coefficient at 
wavelength A. if it is assumed that each absorbed photon creates 
n^(A) electron-hole pairs, then the resulting current density is 

Jg(x,A) - q.T(A).Npl^^(A)-n (A).[l-exp(-a(A)x)], (3.2) 

where 1(A) is the surface transmission coefficient at wavelength A. 

The actual generation rate is 

2 ^^ 3Jg(x,A) 


G (x, ) > - 
e * q 3x 


(3.3) 


or 


G^(x.A) - T(A)-n^(A)-Np^^(A)a(A)exp(-o(A)x). (3.4) 

For a full spectral Irradlknce, the generation rate at each spectral 
point must be added over the Incident wavelengths to give 


°e^*^ “ yT(A)nA(A)Nphi(A)a(A)exp(-a(A)x). (3.; 

If the incident light is in other than the normal direction, the 
actual traversed path is lengthened by a factor of 1/sin and 1/sln e 
of the device width for the first and second refracted light rays 
respectively. These factors are shown In Figure 3.6 as a function of 


V/AVELENGTH (uM) 

Figure 3,6, Factors of l/sin9^, and l/sinB^ as^^a function of wavelength. 

silicon surface (solid lines), 600A Ta_0 antl-yef lection film 
(dashed lines). 




47 


photon wavelength. Therefore the absorbed photon density 
becomes 

. * - 

(3.6) 

and the resultant generation rate becomes 

r 

1.1 

f{ 

where T^, are the coefficients of transmission and reflection with 
respect to the i-th refracted light, 
j The calculated generation rates are shown in Figure 3,7 and 3.8, 

It is seen that a CNR cell has a higher generation rate near the 
illuminated surface within about 20 pm, and has a lower generation rate 
at distances deeper than about 20 pm. This is because more carriers can 
be generated near the surface from the oblique transmission of the 
incident light. The comparison of the generation rates for the plane 
and serrated surfaces is presented in Figure 3.9 at one particular 
wavelength of 0.9 pm and 100 mW power density without any anti-reflection 
film. The crossing point of the generation rates for both cases occurs 
around 40 pm from the illuminated surface. 

The maximum possible collection efficiency of a solar cell is 
plotted against the cell depth4n Figure 3.10. The CNR cell has a 
steeper slope of the collection efficiency vs. width curve. This implies 
that more carriers are generated near the irradiated surface. Therefore 


GENERATION RATE C#/CC/SEC) 






^ DEPTH CUM) 

Figure 3.8. Generation rates of AMO and AM? solar spectrums for a plane 
surface cell or a CNR cell. 
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Figure 3*9. 


depth ( uM ) 

Comparison of the generation ratp<5 itf j 

surfaces at a wavelength of Q 9 and serrated 

serrated surface (dashed). * ’ ^ne surface (solid). 
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Figure 3.10. 


WIDTH (uM) 

Collection efficiency as a function of cell depth for CNR ceil (solid) and 
plane surface cell (dashed) at a wavelength of 0.9 pm. 


cn 



a CKR cell u .ore reaUrion resUtent a. ccperea to a o„avear.ooal 
cell. This poaar „U1 aiscossea 1„ .ore derail later. Table 3 3 
SHOWS the ceu widths where 90 or 95 percent oT the available current 
densities can be potentially collected for both type cells. 

^le 3.3. Widths for 90 and 95 percent collection efficiencies. 


Width (urn) („„) ( , 

Collection Bare S . 

waciency Plane He AnS-Ref . ‘ 600^X^05 


(um) 

Bare Plane 
Si 

^=0.9 jjm 


90 

90 

80 

60 

65 

95 

250 

240 

210 

80 

Available 




Current 

Density 

34.2 

46.40 

51.8 



(Mm) 

CNR 

A=0,9jjm 


45 

60 


3.4 Violet and CNR Cells 

The para.eters of the violet and CNR cells used In the present 
calculations are listed In Table 3.4. The para.eters of this table are 
reasonable When compared to reference [2], A CHR cell differs fro. a 
Violet cell mainly in the surface transmission efficiency as a result of 
®“>^face texture. Also, a CNR cell has a larger surface 

-Won current density and space charge recombination current density 

ue to the enlarged surface area of 2 83 nmo «-u 
. tl.es Chat of the original flat 

surface area. 
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Table 3.4, Models 6i the violet and CNR cell. 


I 


Structxire 

n'*'p (Cell A) 

j 

n‘‘‘pp''‘(Cell B) 

Base resistivity 

2 R* 

) 

1.3 0* cm 

Surface junction 

0.13 urn V 


Surface doping density 

10 /cm , erfc 



Device width 

250 pm 



Device area 

4 cra^ 



SRV 

3 

10 cm/sec 



Two-way reflection 

Yes 



Anti-reflection film 

o 

Ta^Og, 600A 



Sheet resistance 

550 JJ/fa 



Contact resistance 

-4 2 

10 fl* cm 



Sun power density 

AMO, (135.3 mA/ 

2. 
cm ) 


Temperature 

27*C 

^ . - i . 


Surface diffusion length 

Ljj(MED)+Ljj(MIN) 

2 



Base diffusion length 

Ljj(MAX)+Ljj(MED) 

«. 

} 

210 pm 

Grid pattern 

2 

Two tapered bus bar 

20 ]im to 200 um) 

Number of fingers 

60 



Width of fingers 

20 pm 



Thickness of fingers 

1.9 pm 
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The predicted results of the computer analysis are listed in 
Table 3.5. Here the one-dimensional analysis refers to the solution of 

the fundamental device equations in a one-dimensional model. However, 

! ; ■ 

it is more realistic to simulate a practical solar cell by the general 
two-dimensional program of Appendix 10.4. The calculated efficiency 
values can be compared to reported Valpes of 13 and 15 percent for 
violet and CNR cells respectively. The predicted spectral responses 
of violet and CNR cells are shown in Figure 3.11. The experimental data 


Table 3.5. The predicted resulte of the violet and CNR cella. 


Model 

Violet -A 

CNR -A 

Violef-B 


1-DIM 

2-DIM 

1-DIM 

2-DIM 

2-OlM 

2 

I (mA/4cm 
sc 

164.4 ■' 

^ 151.5 

185.6 

171.1 

154.0 

V (volts) 
oc 

0.582 

0.582 

0.585 

0.582 

0.590 

CFF 

0.821 

,0.802 

0.821 

0.779 

0.803 . 

(mA/4cm^) 

153.9 

143.7 

177.2 

161.4 

146.1 

(volts) 

0.510 

0.492 

0.504 

0.481 

0.500 

Pj, (mW) 

78.5 

70.6 

89.3 

77.6 

73.0 

EFF (%) 

14.5 

13.1 " 

16.5 

14.4 

13.5 
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xgure 3.11. Spectral responses of violet and rKn? tt- , 

cell (solid liees) . o.d are rhe reported^aW 


of tl^e spectral response are taken from reference U. The quantum yield 
is then deduced from Equations 2.12 and 2.14. It is clear from 

Figure 3.11 that there is a general agreement between the calculation 
and the experiment. 

It has been pointed out that there are three different forces 
acting on the charged particles in the diffused surface region [ 9 ] . 

They are the electrostatic force, an additional force due to the 
position-dependent band structure and the diffusion force from the 
optically generated carriers. In a solar cell, the short wavelength 
photons can generate large carrier densities near the surface and the 
diffusion force is qu& large in this region. Our calculation shows 
that the minority current is in the direction of the collecting junction 
despite the large retrogate field force from the position-dependent 
band structure in the heavily doped surface layer. This suggests that 
the diffusion force can counteract the internal retrogate field force. 
Therefore, most of the carriers generated in the surface region can be 
collected instead of being lost by recombination at the surface. 

Within the base region, the calculation shows a much shorter 
collection width than the device width. Collection widths of 85 ym 
and 95 pm are calculated for type A violet and CNR cells respectively. 

This is because the back surface ohmic contact is in competition with 
the p-n junction for the collection of the light generated carriers. 

The available device short circuit current density over the above 
collection width can be readily calculated. These calculated results 
are 42.3 and 47.6 mA/cm^ respectively, which are about 90.5 and 92 percent 
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of the total available current density for violet and CNR cells 
respectively. These values are close to the predicted short circuit 
current densities of Table 3.5 Indicating only a small percentage 
of the availabU. current density within the collection width is lost 
through bulk and depletion region recombination. 

3.5 Models of the New Generation Violet and CNR Cells 
Since 1977 violet and CNR cells have been made by combining the 
technologies of the back surface field and the early violet cell of 
1973. The difficult back diffusion process can then be controlled 
and the back contact has also been made highly reflective in order to 
Improve the long wavelength response [10] . 

Typical parameters for such cells are listed in Table 3.4 as cell B. 
The base resistivity of a violet cell is nominally between 1 to 3 fl-cm. 
However, measurements indicate that the majority of cells have a 
resistivity around 1.3 fi-cm. The base diffusion length is taken from 
the measurement of reference [H] where a typical value of 210 pm is 
given for the preradiated cell. The radiation damage coefficient 
for the diffusion length is assumed to be 2.5x10"^® at a doping 
density of 10^^/cm^. 

The short circuit current density under radiation has been 
calculated and compared to the reported measurements for a typical 
violet and CNR cell. These are the type B cells of Table 3.4. The 
results are shown in Table 3.6. 
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Tables 3.6, Comparison of the short circuit current density for 

typical violet and CNR cells* The experimental data are 
taken from reference [10], 



Radiation 

Dose 

0 

xo« 

10^** 

10^® 

5xl0“ 


Experiment 

(mA/cm^) 

39.11 

38.31 

36.45 

31.8 

28.3 

Violet 

Cell 

Calculation 

(mA/cm^) 

40.62 

39.56 

36.46 

31.57 

27.28 


Agreement 

(%) 

+2 

+3.2 

+0.02 

-0.7 

-3.7 


Experiment 

(mA/cm^) 

46.98 

44.43 

41.9 

37.27 

33.83 

CNR 

Cell 

Calculation 

(mA/cm^) 

45.85 

44.81 

41.81 

37.08 

32.92 


Agreement 

-2.4 

+0.8 

-0.2 

-0.5 

-2.8 


(%) 
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The calculated short circuit current density has been corrected for 
the 6 percent metal coverage and the magnitude of the simulated AMO power 
density of 140 mW/cm^. The small discrepancy at low radiation dose is 
perhaps due to the unknown base diffusion length for the particular cell 
measured. The discrepancy for heavily radiated cells may be due 
to an annealing effect before taking the measurement. However, it has 
been found from a comparison of baloon-flown cells and ground measure- 
ment that solar cell response is a few percent less in the space measure- 
ment [12] . This is consistent with our calculation and may suggest 
that there are problems in the simulated solar spectrum in the short 
wavelength region. 

The calculation of spectral response under radiation is shown 
in Figure 3.12 for a typical violet and CNR cell. This is in very good 
®S^®®ment with Figure 8 of reference [11] 


' 3.6 Conclusions 

The following conclusions can be drawn based upon the work in this 
section . 

(1) High efficiency silicon solar cells can be made by increasing the 
surface optical transmission efficiency through a textured silicon 
surface. Such textured surfaces have been found to produce the highest 
surface transmission efficiency, and it has also been found that the 
thickness of the anti-reflection film is not very critical for near 
optimum surface transmission. 



SPECTRAL RESPONSE (MA/MW) 





WAVELENGTH (uM) 

Figure 3.12. The c^culated spectral responses for violet (dashed lines) 
and CNR (solxd lines) cells under one Mev electron radiation 
doses, d - ideal spectral response for infinite cell depth, 
e — ideal spectral response for actual cell depth. 
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(2) It has been found that a high efficiency slliooh solar cell can be 
nade by combining the technology of the BSF structure and the shallow 
surface doping density. Good collection efficiency can be obtained 
with a practical base region diffusion length which has a value In the 
range of the device depth. 

(3) CNR cells are found to be more radiation resistant than a standard 
planar surface cell. This is due to the fact that more carriers are 
generated near the collecting junction in a textured surface cell. 

Hence it degrades slowly with high energy particle irradiation. 

3.7 Appendix A 

The oblique optical reflection and transmission coefficients from 
an anti-reflection film on a silicon surface [13’].’ 
a) Reflection at the boundary of an absorbing medium 



Figure A.l The oblique optical reflection and transmission 
on a silicon sur.face. 
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Let the medium have the complex index of refraction N-nfik and the 
complex propagetion vectcr of the refracted wewebe 1^. h ^ ^ 

ref le'ctlhnT houndai^i^^ 


7" i 

k •Y ■ k’ *Y 
o' o' 


HqTY ■ K*y ■ (k+ia)*Yj 


MiJ tiv* ' : 


(3.8) 

(3.9) 


k sin6 ■ ksin^. 


The wave equation can be written as 


2 .2, 


A . El 3 E 

V t ■ -r — - . 


2 2 
C 3 t 


For plane harmonic waves we have V « -ik, ^ -iw, so 


nV. 

° o 


By equating the real and imaginary parts, we have 


, 2 2 , 2 2v 2 

k -a » (n -K- )k^ , 


k*a « kacosi^ « nKk . 

o 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


After some algebraic manipulation of Equations (3*12) and (3.13), it is 
found that 


kcos(J)+ia « kA^?IsliA 


(3.1A) 


Snell law of refraction can be written 


as 


n^ sine ■ Nsln^ 


(3.15) 
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Now (|) Is a complex number and has no direct correlation to the angle (|i. 
By Equations (3,14) and (3.15) 


jj . kcos<t)+la 


k COS(J» 

■ f 

The boundary condition of continuity of the 


(3.16) 


tangential components of 


electric and magnetic fields leads to Fresnel’s equation for TE and TM 
polarization. 




nj^cos0-Ncos(j) 
^TE n^cosG+Ncos(J) ’ 


(3.17) 


-Ncos0+n^cos<j) 

^TM Ncos0+n^cos(J)* 

If there is an antireflection film, the transmitted and reflection 
amplitudes need to be added as Fig. A2 shows. In this figure, 



Figure A. 2 The multiple reflection and transmission across an 
anti-reflection film on a silicon surface. 
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Where y'"-Y and t'=-t from the Fresnel coefficients of Equations (a. 17) 
and (3.18). Hence the total reflected amplitude is given by 




Tl+Y;e~^% 


Where is the change in phase of the beam on transversing the films 
and 


'^1 * nld;^cos♦J^. 

Finally, it is the energy instead of amplitude which is needed and the 
resultant total reflection and transmission coefficients are 


(3,20) 

I-Real(^)TT* (3 2 ^, 

The exact coefficients of reflection R and transmission T could be 
solved from Equations (3.10), (3,12), (3.13) and (S.16)-(3.19) using the 
apppriatfi complex refractive index of and n 3 for TE or TM modes. 

The complete equation, however, is cumbersome, and hopelessly complicated. 
The desired values of R and T however can be calculated numerically 
using complex algebra. 
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4. SHORT CIRCUIT CURRENT DENSITY AS A FUNCTION OF 
TEMPERATURE AND RADIATION INTENSITY 


4.1 


Introduction 


Although solar cells were conventionally measured at room tempera- 
ture, the operation of solar cells is usually under varied temperature 
and solar Intensities. The cell temperature may vary from 60®C to 90*C 


in the earth's synchronous orbit, and there -are extreme temperatures 
of -120^C and +140“C near Jupiter and Mercury respectively. In the 


terrestrial operation of solar cells, the corisideration of temperature 
is also Important because of the large temperature variations with 
respect to locations and seasons. Especially in a multi-sun concentrated 
system, a high working temperature is usually inevitable from the higher 
input solar intensity. 

The subject of the temperature dependence of short circuit current 
density is tlie main topic in this section. The material parameters as 
a function of temperature will be discussed, and a comparison between 
theory and experiment will also be presented. ' 


4.2 Material Parameters as a Function of Temperature 
Several important material parameters, which determine solar 
cell operation as a function of temperature, are (a) intrinsic carrier 
density (b) diffusion length and lifetime (c) absorption coefficient S 
and refractive index ^d) mobility and diffusion coefficient. 

(a) Intrinsic carrier density n^. 

The intrinsic carrier density n^ is a function of temperature and 
bandgap energy. The general form is [1] 



-E CT)/2kT 

n, ■* e ■ * 


c V 


4.9-10 *C- 


,15 3/2 -EgCT)/2kT 

•C 5“1 *T ve * 


n 
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(4.1) 

C4.2) 


“de “dh ®^^® *^^® ‘J“»8tty-of-state effective masa of electron 
and hole and E^(T) is the bandgap energy at temperature T. The empirical 
form of the temperature dependent intrinsic carrier density » which 
satisfies the above function of temperature, is used in this work [2J 


n^ = 3.73‘10^^*T^/2.g-"7014/T^ 


C4.3) 


The open circuit voltage of a solar c^ is found to decrease 
with temperature. This is due to the strong temperature dependence of 
the intrinsic carrier density. In the first order model the open circuit 
voltage can be defined by 


kT 


■J 


s 


(4.4) 


where 1^(1) is the saturation current density which is proportional to 
2 

®i ”i in Shockley’s diffusion current model or Sah’s recombination 
current model respectively. The diode factor n is 1 or 2 for the above 
models respectively. Hence^ the reduction of open circuit voltage as 
a function of temperature can be derived from Equations 4.3 and 4.4 with 
the same result for both cases. >».. 




(4.5) 


or 


AVq^( 273"K,T) -1.9 mV/K“ 


(4.6) 
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However, both current mechanisms may exist simultaneously in a practical 
silicon solar cell. In this case AVq^(T) can be derived assuming an 
abnormal diode factor A and the dominance of diffusion current near 

''oc®- 

* VoC«2J-''ocf^l^ - A- ^ ITilogWogT^J C4.7) 

where A is between 1 and 2. Therefore the AV^^fT) v^lue is usually 

lower than -1.9 mV/K“ depending on the A value of a practical solar cell. 

The experimental measurements always show a negative temperature 

coefficient of V._(T) with, a value of -2.2 to -2,3 mV/K® for 10 0-cm 
oc 

and 2 fl-cm silicon cell respectively [3]. These values are in agreement 
with other reports. Luft reported a value of -2.13 to -2.29 mV/*K for 
10 ft-cm cells [4]. A value of about -2.2 mV/K® was reported by Yasi [5J 
for a 2 0-cm cell. For 6aAs, the reported data is between -1.9 and 


-2,2 mV/®K [6,7] . 

The negative temperature coefficient of the open circuit voltage 
can partly be compensated by the small positive temperature coefficient 
of short circuit current. However, this compensating effect is small due 
to the very small temperature coefficient of the short circuit current 
density. 

(b) Diffusion length 

The magnitude of the diffusion length is determined by the carrier 
lifetime and the diffusion coefficient as. 


L = VD t 
n n n 


C4.8) 


A similar formula can be vnritten for holes by substituting sub-index p 
for n. The diffusion coefficient is calculated from Einstein relation 


' .for non-degenerate semiconductors as 
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I. . 


kT 

D * ~ u * 

n q 


C4.9) 


where the mobility la generally determined for silicon by the 

acoustical phonon scattering mechanism and the impurity scattering 

mechanism. The component of lattice scattering has the theoretical 

“3/2 

temperature dependence of T , While a temperature dependence pf 
—2,5 "2.7 

T * and T * is experimentally observed for electrons and holes 
in a lightiy doped material [8J . The component of impurity scatter- 
ing has a temperature dependence of 

In this work total mobility is calculated from lattice scattering 
mobility and impurity scattering mobility as 


- . -V ; ^ -t 


JL 


(A.io) 


where 


17 

M, - 4.7x10^' 


p3/2 


and 


'eff Npln(l+4*5xlO®T^/Njj^^^) 


■T,V 


(4.11) 


Ln 


.9„“2.5 2 


- 2«lxl0 T * cm /V» sec (electrons) , 


y^p = 2.3xlO^T ^*^cm^/v-sec (holes). 


(4.12) 

(4.13) 


For an indirect bandgap material such as silicon, lifetime is primarily 
controlled by defect levels located in the forbidden bandgap. If the 
single defect level of Shockley and Read's model is assumed, the lifetime 
of the minority carriers can be shown to be 


■"p “ ''po^^ ®*PCEt“Ef^)/kTj 
in the n-type region and similarly 


(4.14) 
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i 

in the p-type region, where t and x\ are equal to l/o v and 

po no p th t 

respectively. Hence the carrier lifetime depends on the 
relat^ive location of the trap level with respect to the quasi-Fermi 
level and the intrinsic energy level. In early reports, the lifetime 
was found to increase exponentially with temperature as Equation 4.14 
and 4.15 for a shallow trap level [9]. However; more recent measurements 
indicate the dominance of a deep recombination level near the" intrin- 
sic energy level and this implies that the above exponential terms 
of Equations4.l4 and 4.15 can be neglected [10]. The lifetime will then 
be approximately constant for operation near room' temperature. These 
considerations suggest a relatively constant diffusion length as a 
function of temperature for the deep recombination level model of a 
silicon solar cell* The published results of the measured diffusion 
length as a function of temperature demonstrate the almost constant 
diffusion length except at very low temperature [11]. 

(c) Absorption coefficient and refractive index 

The temperature dependence of the absorption band edge in a silicon 
material is a well established physical phenomenon. The temperature 
dependence of the absorption coefficient can, in general, be obtained 
by the appropriate displacement along the horizontal and vertical axes 
in the absorption-wavelength curve. The horizontal displacement is 
basically due to the bandgap shrinkage as temperature Increases . The 

empirical value of the bandgap temperature coefficient for silicon is 
-4 

about -2.4x10 eV/“K at temperatures above 150K* [1]. 
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Since silicon is sn indirect bandgap jnaterial, the photon 
excitation of the absorption mechanism requires phonon participation 
for the conservation of momratum and energy. Therefore the temperature 
functional dependence of the vertical shift can be deduced from the 
phonon statistics assuming McLean’s model of two phonon interaction 
with energies at temperatures .of 212*K and 670®K respectively I12J . 

The absorption coefficient will then increase or decrease with 
temperature beyond the 300’K curve by a factor f in the vertical 
direction of [13] 

_ ^ + 712 ^ + :6757iJ 

' 1^ -1 . 1 • (4.16) 

^12/300_^ * ~67o7300 “ 

Since the absorption coefficient Is determined by the direct bandgap 
for photon energy greater than the direct bandgap energy, vertical 
movement of the a-A curves Is negligible for these high energy photons. 
Several absorption-wavelength curves at different temperature can 
then be obtained by the appropriate horizontal and vertical movement 
with respect to that at room temperature. These are shown In Figure 4.1. 

Measured data of Bash and Newman [14] on the photon absorption of 
silicon are replotted In Figure 4.2 as versus («v-E ) at 

temperatures of 300»K and 77"K [15] . The straight line matching of the 
1 / 2 ' 

a -(-liv-Eg) relation is a well known characteristic of the indirect 
optical transition. Our analysis predicts a good fitting to the 
experimental data ext-ept in the short wavelength range, where, the 
direct bandgap photon transition becomes dominant over the indirect- 
bandgap photon transition. 


BBSORPTION COEFFICIENT (/CM) 



WfiVELENGTH OlM) 

Figure 4.1. Absorption coefficient as a function of temperature by 
the proper horizontal and vertical movements In the 
a-X plot. 




ENERGY (eV) 


ulated absorption coefficient 
rement at 77®K. 
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The general correlation between calculation and experliiient 
shows the ability to determine the absorption coefficient as a function 
of temperature by suitable horizontal and vertical translations of the 
a— X curves. Since the major portion of the solar spectral intensity is 
located at photon energy below about 2.5 eV, the deviation of the 
predicted absorption coefficient for higher energy photons has only small 
effects on the predicted short circuit current density. 

The temperature dependence of the refractive index is taken from 
refierence [16] for the wavelength range from 0.18 pm to 1.1 pm. 

However, the temperature coefficient of the refractive index has a very 
small value of 5xl0“^ 


4.3 The Temperature Dependence of Short Circuit Current as a 
Function of Base Doping Density and Radiation Dose 

The variation of the short circuit current as a function of 

temperature can be conveniently defined by the normalized temperature 

coefficient K^. of the short circuit current 



^SC«o> 


«SC<"o> 

dT 


(4.17) 


where T^ is a convenient reference of temperature such as 0“C or 27“C. 
Several reported experimental results for conventional silicon solar 
|j:ells h^ve been compiled in Tables 4.1 and 4.2 [3,4,17,18]. 


Thtioreti-cal calculations were made to compare with this experimental 
data, ite diiffusion length was kept constant with respect to the 
tempemtuta ^farlatlon in the calculation. (A few calculations using con- 
stant lifefTLme ,:re also attempted. However they have given identical 


results.) A modified generation rate, which is calculated from the lin- 


vertically shifted absorptipn coefficient, was used 
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A Violet-cell model has been assumed to be a typical structure of 250 pm 
thickness, 0.12 pm Junction depth, and 0.5 pm high-low junction with a 
600 A thickness of Ta20^ for an anti-reflect^ion film. The base layer 
has been assumed to have a doping density of 10^^/cm^ and the surface 
layer has a doping density of lO^O/cm^ with a profile pf a complementary 
error function. The model of the conventional cell has been assumed to 
have^a structure of 300 pm thickness and 0.25 pm junction depth with an 
800 A thickness of SiO. The base and surface doping density are assumed 
to l^ve values of l.SxlO^^/cm^ and 2xlo20/cm^ respectively. The one MeV 
electron radiation damage coefficient is assumed to be 9xl0"^^ 

e/cm^ and 2.5x10 e/cm^ for 10 fi-cm and 1.3 fi-cm p-type silicon 
substrates respectively. 

Several characteristics of Che calculated normalized temperature 
coefficient of short circuit current are found to be in agreement with 
experimental measurement. The I^^-T curves were found generally not 
to be a linear function of temperature. The nonlinear behavior 
intensified at larger radiation doses. This characteristic has been 
reported previously [19]. Also it was found that the normalized 
temperature coefficient of short circuit current has a larger value at 
a heavier radiation dose. This is due to more photon-generated carriers 
available near the surface through the absorption coefficient change 
with temperature. Hence the rate of Increase of short circuit current 

with temperature is higher for the solar cells with a shorter diffusion 
length after radiation bombardment. 

The compiled data are shown in Figure 4,3 along with the calculated 
results. In Tables 4.1 and 4.2, the temperature coefficient of short 



Table 4*1. Normalized temperature coefficient of 
silicon solar cell. 
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Figure 4.3, Normalized temperature coefficient of short circuit current density as a function 
of one MeV electron radiation. 
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circuit current was measured within the temperature range AT and 
normalized to IggCT^). Despite the minor differences of these measnre»- 
ments. the mean valoe and standard deviation of the normallred 
temperature coefficient were calculated and shown in Figure 4.3. 

On Che other hand, the numerical calculation was performed in two 
temperature ranges of 27“C to 77»C and 27"0 to 127»C respectively. 

The normalized temperature coefficient Kj was then calculated within 
these temperature ranges and normalized to 1^^ (27°C). Therefore, a 

mean value and standard deviation of can be calculated with respect 
to the above temperature ranges. 

The predicted value of is higher than the measured value at the 
pre-irradiation condition. This discrepancy could be due to an Improper 
initial diffusion length used in the calculation. However the comparison 
between measured and calculated values is pretty good at high radiation 
dose density where the diffusion length is primarily determined by the 

radiation damage coefficient irrespective of the initial value 
assumed. 

4.4 Discussions and Conclusions 

Although the small positive temperature coefficient of short circuit 
current has been a well known experimental result, thd physical' 
explanation has been unsatisfactory. Wysockl and Eappaport discussed 
the dependence of diffusion length on temperature [20]. Luft suggests 
that it is the Increase in the minority carrier lifetime with tempera- 
ture which improves the ted response of the long wavelength photons. 
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He also surmised that the higher normalized temperature coefficient 
after heavy radiation bomhardment is due to the stronger tempera” 
ture dependence of the carrier lifetime for the heavily radiated cells 
than fot the lightly radiated cells. However, the temperature dependent 
absorption coefficient and generation rate were found to have the 
dominant influence in this work and these were not considered by 


these previous authors. 

One factor, which may affect the temperature coefficient of short 
circuit current, is the reduction of the forbidden bandgap and 
consequently its effect on long wavelength response. However, the 
chances of collecting these long wavelength photons are very small due 
to the low absorption coefficient of 10/cm or less. The short circuit 
current contributed from this spectral range can be calculated by the 
Equation [21] as , 






1 + aL 


D 


-W/L -o/Lj, 

j. e -g De , 

*■ (l-aL^)coshW/Ljj 


(4.18) 


The calculation shows a negligible effect because of the small solar 


spectrum density and small absorption coefficient in this long wave- 
length range. Another factor, which has been suggested as responsible 
for the temperature dependence of short circuit current, is the high 
temperature annealing effect of the minority carriers [55] . However, 
the above cited experiments except reference [3] were carefully performed 
near room temperature to avoid any annealing effects. 
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5. KELDYSH-FRANZ EFFECT AND SILICON SOLAR CELLS 
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5.1 Introduction 


The effect of high doping density on limiting the 
open circuit voltage of a silicon solar cell has been 


efficiency and 
recognized for 


many years. Experiments have shown that a low resistivity cell does not 
produce high open circuit voltage. This is not consistent with the 
simple theoretical predictions. However it is generally believed that 
Increasing the output voltage is a key factor for obtaining a high 
efficiency solar cell [1,2]. Therefore, the real physical limitations of 
low open circuit voltage need to be identified for low resistivity 


solar cells. 


Many hypotheses have been given to explain the inconslstancy 
between the predicted and the achieved voltage performance of a low 
resistivity silicon cell. The heavy doping effects have received much 
attention [3]. There are also other explanations and speculations. For 
example, it has been suggested that an excessive recombination current 
occurs in a heavily doped cell and the open circuit voltage is thereby 
reduced. However, tl. \ measured short circuit current and diffusion length 
of low resistivity cells has shown that this is not the primary factor in 
limiting the cell performance [3]. Recently the so-called Keldysh-Franz 
effect has been postulated and examined [4], 

The purpose of this section is to investigate the low output voltage 
performance of a low resistivity silicon cell and to assess the relative 
influence from the I|eldysh-Franz effect. The Keldysh-Franz effect is 

'jl 

modeled by an equivalent bandgap reduction in the depletion region from a 

I'l 
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high Internal electric field. Therefore the temperature dependence of 
the forward and reverse dark I-'V characteristics can be calculated and 
used to Identify the real limitations by several proposed physical 
effects. 

5.2 Keldysh-Franz Effect 

In 1958 Fr^nz [5] and Keldysh [6] Independently calculated and 
predicted the so-called Keldysh-Franz effect; which Is tunneling assisted 

i] 

light sbsorption. The essential idea of this theory is that the elec- 
trons tunnel Into the forbidden ^bapd before optical absorption can occur. 
This .results in a broadening of absorption and can be treated like a band 
gap shift with the electric field of 

2 

( 5 . 1 ) 

This shift of the fundamental absorption edge has been observed for many 

semiconductors such as Ge; SI, GaAs, etc. The space charge recombination 

current density will thereupon Increase by a factor of exp(AE (E)/2kT) as 

S 

qn. 2slnh(qV,/2kT) 

' " ( ' v^-Vj)q/kT • f(b)»exp(qAE^(E)/2kT). (5.2) 

Therefore the recombination current becomes important for high bandgap 
and heavily doped materials and at low temperature because of this 
Keldysh-Franz effect. These phenomena are generally consistent with 
the experimental observations. 

In this work an empirical value of the bandgap reduction as a 
function of internal electric field has been taken from the work of 
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Brltsyn and Smirnov [7]. Their results have confirmed the functional 
dependence of the Keldysh-Franz effect on the electric field strength. 

5.3 Junction Analysis and the Electric Field in the Depletion Region 
At a p-n junction there exists a strong Internal electric field from 
the potential barrier. It is knovm that such an Internal electric field 
can alter the bandgap through the Keldysh-Franz effect. Thus the forward 
or reverse dark current density will be increased by the enhanced deple- 
tion region recombination current density. From Equation (5.2) it is 
clear that the recombination current density is higher at low resistivity 
and low temperatures. In order to predict the importance of Keldysh- 
Franz effects on a silicon solar cell, the complete Poisson equation 
must be numerically solved. The effective electric field is calculated 
for the unperturbed or perturbed bandgap along the p-n junction. The 
forward and reverse dark current densities are next determined. The 
effective electric field in the depletion region Is plotted in Figure 
5.1 for the model cell of Table 5.1. 
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Table ,5.1. Material and . dlmenalonal parameters . 


Table 5.1 Material and Dimensional parameters 
Cell Thickness 
Cell Structiire 
n^ Thickness 
n^ Doping concentration 
. p Doping concentration 
Diffusion length in p region 
Diffusion length in n region 
Sxirface recombination velocity 


250 van 
+ 

n p 
1 ym 

20 3 

4x10 /cm , erfc 
5x10^*^/ cm^, (0,1 fi-cm) 
100 v>m 

(Lp(MED) + Lp(MIN))/2 
3 

10 cm/ sec 


The characteristics of the electric field in the depletion region 
is changed after considering the Keldysh-Franz effect. The peak 
electric field is reduced; while the strong field region is broadened. 
This is a direct consequence of the bandgap variation with the electric 
field in the depletion region. The peak electric field is located near 
the n"^p metallurgical boundary where there is the greatest bandgap 
reduction. Therefore the peak electric field will be decreased due to 
the lower forbidden barrier height. 

The dark current density is presented in Figure 5.2 and 5.3 for the 
forward and reversed biased conditions respectively. The higher current 
density after considering the Keldysh-Franz effect is basically from the 
enhanced recombination current density in the depletion region. It is 
estimated from Figure 5.2 that the reduction of open circuit voltage due 
the Keldyshr-Franz effect is approximately 6, 8 and 11 mV at temperatures 
of 350*K, 300*K and 250“K, respectively. 


Figure 5.1. 
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5.4 Heavy Doping Effect 

It has been recognized that the heavy doping effects on minority 
carriers can be treated by either an effective intrinsic carrier 
concentration or by an effective doping density [2,8). The reduced 
effective surface doping density of a conventional silicon solar cell 
is found to be responsible for the low voltage output of a low resistivity 

cell [1], This also leads to the reduction of the emitter efficiency and 
current gain 3 cxf a transistor [9], 

Figure 5.4 shows the effective surface doping profile as a function 
of temperature. The characteristics of the effective doping density can 
be used to explain the influence of emitter doping on the temperature 
sensitivity of the emitter efficiency of a transistor [9]. However, this 
point will not be discussed here. The reduction of the output voltage of 
a highly doped cell is a consequence of the lower effective surface 
doping density which greatly increases the surface back injection cur- 
rent density. At the same time short circuit current density is decreased 
by the retrograde field* near the surface. 

The percentage of the surface injection current density to the total 
current density is shown in Table 5.2 at forward voltage of 0.5, 0.6 and 
0.7 volts for temperatures of 300°K and 350°K respectively. It is clear 

that the back Injected current density has an abnormally large value for 
a low resistivity cell. 



EFFECTIVE SURFACE DOPING (#/CM' 


Figure 5 




Table 5.2 Percentage of surface Injection current to the total current' 
density as a function of temperature . 


V( volts) 


350°K 


300°K 


5.5 Discusslonsf and Conclusions 
Although the so^-called Keldysh-Franz effect was originally 

investigated for the optical absorption properties of a semiconductor 
or insulator, this effect on the efficiency of a heavily doped silicon 
solar cell is not clear. Our calculation shows that this effect is not 
very important as compared to the heavy doping effect for a cell with 
0.1 n-cra base resistivity. On the other hand, the Keldysh-Franz effect 
is found to be important in a first order calculation which assumes 
an unperturbed bandgap structure. However our model takes into account 
the perturbed bandgap with an effective electric field. The effective 
electric field is thereupon reduced because of the smaller bandgap in 
the depletion region. Therefore, these calculations predict a small 
influence due to the Keldysh-Franz effect for a silicon solar cell. 
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6. THIN SILICON SOLAR CELLS 
6.1 Introduction 

For years there has been a great interest in experimentally pro- 
ducing very thin solar cells in the range of 50 ym. Such cells can be 
used to minimize power loss in the base region under large solar concen- 
tration. They are also required to produce lightweight space solar 
power systems. The interest in thin solar cells is also enhanced by 
the economical consideration that there is less kerf loss from the ingot 
cutting process. The cost from the silicon wafer can thereupon be 
minimized. Recently thin silicon solar cells have been fabricated using 
an improved low kerf loss ingot cutting process [1] and by non- 
preferrential etching in a NaOH solution [2], 

These new technological breakthroughs are inevitably leading to new 
interest in a' better understanding of thin cells, where there are more 
interactions between the incident photons and the back surface. High 
injection effects and the high-low junction leakage current can be 
studied in such structures. There are also reflections in the optjipal 
spectrum between the front and back surfaces in such a thin cell. The 
purpose of this section is to investigate the physical behavior of a 
thin solar cell. The conversion efficiency of a thin solar cell under 
one MeV electron irradiation is also Included for cells of several base 
resistivities. 
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6.2 Some Characteristics of .a. Thin Solar Cell 
A thin solar cell is a good structure for investigating the physical 
properties of the interactions between the front and back surface, due 
to the incident photons as well as the injected carriers. 

BaS‘'Srfacr^''*''^°''^ Between the Long Wavelength Photons and the 

A good BSF cell requires a high optical reflectance at the back 
surface and it shou|d also provide a low surface recombination velocity 
in order to improve the long wavelength response. Unfortunately, these 
two requirements are experimentally found to be in conflict with each 
other [3]. An optically absorbing alloy interface is usually found to 
provide a higher electron barrier and a lower surface recombination ‘ 
velocity than an optically reflecting Interface. The optical 
reflectance for alloy contacts has been measured in the range of fifty 
to eighty percent depending on the alloying time and temperature for 
wavelengths from 1.3 ym to 1.4 ym [3]. 

It is known that the available optical current density is a function 
of device thickness and the optical reflectance at the alloyed back 
surface. Hence it is possible to calculate an' ideal caHectipn effi- 
ciency as a function of device thickness and the back surface reflectance 
for a particular solar energy spectrum. Table 6.1 shows this character- 
istic for AMO and AM2 solar spectrums. 
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Table 6.1. Available current density as a function of back surface 
reflectance for a device thickness of 50 ym. 
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Table 6.1 shows that the improvement in the available current den- 
sity is only a few percent even for a perfectly reflecting back surface. 
On., the other hand, it has been, shown that the minority carrier collection 
depth is a strong function of the back surface recombination velocity. 

The collection depth is usually found to be much shorter than the device 
width due to recombination at a poor high-low junction. The reduction of 


the collection depth results in a smaller short circuit current for a 
thin solar cell. Also a poor high-low junction enhances the interaction 
between the injected carriers and the back alloyed surface and produces 
a low open circuit voltage and low conversion efficiency. Therefore, it 
can be concluded that a good high-low junction is much more Important 
than good optical reflectance in improving the conversion efficiency of 
a thin solar cell with cell thickness of 50 


ym or more. 
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6.2.2 Injected Carriers and the Back Alloyed 

It has been shown that the high-low junction theory is an adequate 
model for a BSF cell [4,5].^ It has also been found that the high-low 
Junction barrier height can be lowered for high Injection operation [6]. 
Hence the high-low junction leakage current cannot be neglecte^d^for a 
solar cell trith long base diffusion length and short cell thickness [7]. 
Figure 6.1 shows the prediction of open circuit voltage as a function of 
the back surface recombination velcoity of a 10 J^-cm cell with the para- 
meters of Table 6.2. Curve (a) is calculated using an ideal SRV which is 
a function of the material characteristics on both sides of the high-low 
junction. Curve (b), (c) and (d) are calculated by increasing the SRV 
values by two, ten and one hundred times the ideal SRV valuer. Reductions 
in the SRV can occur in a real solar cell from low lifetime or from a 
high doping density at the high-low junction. It can also occur if there 
is local ' surf ace damage, such as dislocations, etc. 

Since open circuit voltage is a direct indication of the interaction 
between the injected carriers and the back alloyed surface, it can be 
concluded from the result of Figure 6.1 that a good high-low junction 
barrier is important for obtaining high output voltage of a thin BSF cell. 

In Figure 6.1 an exact numerical calculation of V is also 

oc 

presented which shows good agreement with the predictions of first 
order models. 


OPEN CIRCUIT VOLTAGE (VOLTS) 



Figure 6.1. 


DIFFUSION LENGTH (uM) 

yoltap as a function of base diffusion length 
and_high:low junction SRV for thin solar cells with a base 

resistivity of 10 fl-cm. HL junction theory( ), HL juncti 

plus leakage current theory ( ) , exact numeri^l calcula 

tion under condition a (A) . 
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6.3 Numerical Calculations 

Since there are several^ complicated physical interactions between 

incident photons, and minority carriers with the front and back surfaces, 

an exact numerical calculation is necessary to accurately predict the 

performance of a thin solar cell. In order to determine the best 

material parameters for optimum operation, a series of calculations ' 

have been performed. The model parameters are listed in Table 6.2. The 

radiation damage coefficient is deduced from experimental measurements 

and shown in Figure 6.2. The results of calculation are shown in 
Figures 6,3 .to 6,6. 




Table 6.2 



Material p arameters of thi n solar cells 
Structure ~ ^ 


n'*‘density (#/cm^) 

•f* 

n thickness (yM) 
p density (S2«cm) 

P thickness (yM) 

density (i/cm^) 
p thickness (yM) 

Base diffusion length. 
Surface diffusion length 
SRV (cm/sec) 
Anti-reflection oxide • 
Double reflection 
Radiation coefficient 


20 

l.SxlO^^'CerfG) 

0.15 

100 , 10 , 1 , 0 . 3 , 0,1 
^49.35 
19 

10 CGaussian) 

0.5 

(Ljj(MAX)+L^CMED)/2 

(L^(,MED)+Lj^(.MIN)/2 

10 ^ 

Ta^O^, 600A 
YES 

Figure 6.2. 




(%) A0N3I0 




OPEN CIRCUIT VOLTAGE (VOLTJ 
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Figura 6.3 shows the conversion efficiency as a function of radia- 
tion dose. The optimum material is found to have a base resistivity from 
1 n-cm to 0.3 fl.cm for a pre-irradiated cell. However, an optimum base 
resistivity is located between 0.3 «.om to 0.1 «-cm after irradiation. 

The conversion efficiency is found to be essentially the same for 
different material parameters after a radiation dose of 5x10^= e/cm^ 

Figures 6.4 and 6.5 show the short circuit current and open circuit 
voltage as a function of radiation dose. It is clear from Figure 6.5 
that high injection and an ineffective high-low junction are the causes 
Of low open^, circuit voltage for a 100 n-cm cell. ' 

Figure 6.6 shows the ratio of the peak power densities to unlrradi- 
ated power density as a function of radiation dose in order to determine 
the radiation degradation for a thin solar cell. The results show that a 
low resistivity cell has the highest ratio of peak power despite its 
higher radiation damage coefficient. Also the degradation rate is found 

to be smaller for high resistivity cells as indicated by the smaller 

slope of the curves in Fip'iirp» fi r u-. 

Mgure 6.6. However, the conversion efficiency is 

lower at high radiation doses for the^P v, • u.- . 

j.or xnese high resistivity cells. 

6.4 Conclusions 

The characteristics of thin BSF silicon solar cells have been 
analysed In this section. It is found that a good high-low j,n>ctlon Is 
more important than an alloy junction of high optical reflectance for a 
thin solar cell with cell thickness of 50 pm or more. The optimum cell 

resistivity is identified to be from 1 fl*cm to 0 q o k ur 

om 1 cm to 0.3 fi-cm before irradiation 

and from 0.3 fi»cm to 0 1 o.r-m 

• cm after a medium dose of one MeV electron 
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radiation. The radiation degradation is found to be smaller for a high 
resistivity cell due to the smaller radiation damage coefficient. How- 
ever, the conversion efficiency is lov/er for these high resistivity 
cells after heavy radiation. 

// 
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7. OPTIMUM „ SOLAR CELL DESIGNS FOR CONCENTRATED SUNLIGHT 
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7 . 1 Sunmiary 

A computer program has been developed to perform a two-dimensional 
calculation of silicon solar cell performance in order to evaluate the 
importance of series resistance on the conversion efficiency of 
concentrator solar cells. It is shown that' the optimum concentrator 
power density for a specific grid design or the optimum grid design 
for a specific concentrator power density can be quantitatively predicted 

O 

by this model. 

7.2 Introduction 

For a successful exploxtatxon of terrestial solar energy'^ tKe cost 
of a photovoltaic system must be drastically reduced from today'*'S‘ 

P^^oe. Multi— sunlight concentrator systems provide a possible, approach, 
to“solving this problem. Hence it is essential that an optimum cell 
design is achieved in order to economically utilize the solar energy. 

It is known that the collection efficiency is a major factor in 
limiting solar cell efficiency at low solar concentrations. However 
series resistance becomes more important at high solar concentrations. 
Therefore the optimum cell design depends on the solar intensity at 
which a solar cell is operated. Similarly, the optimum solar concentration 
should be used for a specific solar cell designed to be operated under 
concentration. 

This section describes the application of a distributed resistance 
and current model to the optimum design of concentrator cells and the 
grid pattern for multi-sunlight operation of a solar cell. A comparison 
between model and experiment is also attempted. In -addition, the model 
is compared to a first order model of lumped series resistance. 
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7.3 Solar Cell Model 

A detailed two-dimensional model has been reported elsewhere 
(Appendix 10.4). 

7.4 Comparison Between Model and Experiment 

A numerical model of a 2 cm^ bac?Pi • 

X d ^ cm oaseline silicon solar cell has been 

■«ade and compared to experimental values [1] . _ .ihe ceU is reotanaulat 

in shape with dimensions of 1 oni hy 2 cm. The front grid pattern 

consists of 39 fingers. 0.005 cm wide and 0.98 cm long spaced evenly 

across the 2 cm cell dimension. The fingers are connected at one end 

by a common bus bar. The front metal of silver has a resistivity of 

^ about 1.59xl0-^fl.cm" and a thickness of 3 pm. Baseline cells with base 

layer resistivities of 0.3 «. cm and 10 8- cm have been simulated at solar 

intensities ranging from one to ninety AMI power densities. The dark 

current density is described by a single exponential voltage dependence 
vrith a saturation current density of S.SxlO-^^ 

and an equivalent diode factor 1.05 and 1.04 for 0.3 8-cm and 10 8.cm 

cells respectively. The surface sheet resistance is assumed to have a 

value of 120 Sl/0 for a junction depth of about 0 . 25 'um to 0.3 pm. The 

metal contact resistance is assumed to have a small value of 19 -“ n-cm^. 

The short circuit current densities under one AMI solar spectr™ are 

28.0 mA/cm= and 30.5 mA/cm= over the active areas for 0.3 n-cm and 10 fl-cm 
cells respectively. ' 

The comparison of model and experiment is shown in Figure 7 . 1 . It 
is seen that the predicted optimum concentrated solar intensities are very 
close to the experimental values where the peak efficiencies occur at 
about 20 . 25 AMI and 5 AMI solar intensities for 0.3 8.em and 10 8.cm 
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— 5(— EXPERIMENt-v 
MODEL ■ 


I ' To 20 30 40 50 60 70 80 < 

CONCENTRATION (SUNS) 

Figure 7.1 Osteen experiment and two-dimensional model for baseline cells with base 

^2^0^ r "^^^P^rtively. Surface sheet resistLci is 

»d ir-‘*-..crforeirieTair“ ^ data . 
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cells respectively. Figure 7.2 shows the absolute total power 
from metal coverage, sheet resistance, finger resistance, base resistance 
and contact resistance as a function of solar intensities for a 0.3 n-cm 
cell. Figure 7.3 shows the relative percentage of -the power losses from 
each component of the total series resistance. It is shown in this 
figure that the optimum concentrated sunlight occurs at the' point where 
a compromise occurs between the three major series' resistances of sheet, 
finger and base resistances. Also the relative power loss from the grid 
coverage area is lower when the power losses from the series resistances 
become higher at large solar concentration. 

v\ 

7.5 Comparison Between the bumped Resistor Model and the Exact Numerical 

The first order model usually assumes a constant lumped resistance 
connected in series with a solar cell. Therefore the photovoltaic 
current-voltage characteristic becomes 


1 = 1 [expO 


q(V-lR^) 


o ^ A kT ) "'IJ - , 


(7.1) 

where is the saturation current density and is the short circuit 
current density. 

It can be shown that the lumped series resistance can be represented 
by the following equation (Appendix A) 


r. _ _M 

s “ WF.T 


LF" 


R 


3-Aa IQ-A^y-^+A^ 


‘SC 


R 


'BC 


A 


B 


(7.2) 



POWER LOSSES (%) 



t^osMCEMTRATION ( SUMS ) 

Figure 7.2 ^solute power losses of a 0.3 fi-cm solar cell by the gri 

resistance, 

""""" line)/Lu.ped resistance .ode, 


POWER LOSSES (%) 



CONCENTRATION (SUNS) 

Figure 7.3. Relative power losses of the four components of the series 
resistanoe, R . 
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The above five components of Rg are the' finger , surface layer, 

base layer, surface contact and bottom contact resistances respectively. 

Rjj and Rg^ are metal and sheet resistivities respectively, Rg^ and Rg^ 
are the surface and bottom contact resistivities respectively. SF, WF and 

LF are the spacing between the fingers, the 'width of the fingers and the 
length of the fingers, Rg is the base resistivity and Wg is the 

cell thickness. A- and are the surface and bottom contact arenas 

O D 

respectively. A^ is the total active area. 

The comparison of the lumped resistance model to the distributive 
resistance and current model i^ shown in Figure 7.2. It is found that 

the lumped resistance model overestimates the power loss in a nonlinear 
behavior. This discrepancy is due to the fact that the total current 
density is assumed to flow through the total lumped resistor of a 
solar cell in the lumped resistance model. In fact, the resistance 
and current elements of a solar cell are distributive parameters and a 
distributed resistor and current model is a more realistic model. 

In order to determine the discrepancy of the lumped resistance 
model, a more detailed comparison between both models is necessary. 

The numerical calculation of model cell D-7-1 with parameters as 
given in Table 7.1 is shown in Table 7.2. 

Model (a) assiimes that there are voltage drops across fingers, 
sheet surface and base bulk regions. The overestimation of the efficiency 
loss by using the lumped resistance model is found to be about 6 percent 
at fifty suns concentration le\^el. Model (b) assumes there is only base 
layer b\alk resistance. The overestimation of the efficiency loss is 

found to be about 4.4 percent. Hbdel Cc) assumes an additional base 

t 

5 contact resistance, and the over estimation by the lumped resistance 

'"i 

model' is- about 4.3 percent at fiftyysuns concentration, level. 
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Table 7.1 Material parameters of c^H D-7-1. 


Structure 

Surface doping density 
Surface junction depth 
Base doping density 
Base width 

Surface Diffusion Length 

Base Diffusion Length 

Surface Recombination Velocity 

Syp power densi'ty 

Anti-reflection film 

Baseline structure 

Width of fingers 

Thickness of Fingers (T) 

Sheet Resistance (R„ ) 

b r 

Surface Contact Resistance 

Base Contact Resistance (rSa) 

BC 

Base bulk resistance 


+ / C' 

n p Cl.O S2-cm) 

20 

1 . 5x10 i^/ om3 , erf c function 

V 

0.3 uM 

1.6x10 H^/cm (1.0 J^-cm) 

250 yM 

(L(MED)tL(MlN))/2 

(LCMAX)tLCMED))/2 

3 

10 cm/sec 
AM2 C74 mW/cm^) 

Ta^Og, 60.0.A 
39. fingers 
50 yM 
3 yM 

92 (calculation) 

-4 9 

10 fi-cm 

-4 2 0 

10 J2'cm or 0.05 Q'cm^ 


0.025 n 
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Table 7,2. Comparison of the Lumped and Distributive Resistor Models. 
Model a (T = 3 ^m, = 92 , R^^ = lO""^ n-cm^). 


Model b (T = 30,000 ym, R„ =0.01 

b r 

Model c (T = 3 ym, R^^ = 92 , 

^BC “ 

=10“^ J2*cm^ 
05 fi'cm^). 

). 

Sun Numbers • 

1 

5 

50 • 

Models 


Efficiency (%) 

Exact Model a. ' 

15.44 

16.60 

17.48 

Lumped R Model a. 
s 

15.23 

16.43 

16.48 

Overestimation (%) 

1.4 

1.0 

5.7 

Exact Model b. 

15.48 

16.74 

18.53 

Lumped R Model b. 
s 

15.26 

16.56 

17.71 

Over estimation (%) 

1.4 

1.1 

4.4 

Exact Model c. 

15.41 

16.42 

15.49 

■■Lumped Rg Model c. 

15.20 

16.26 

14.83 

Overestimation (%) 

1.4 

1.0 

4.3 


Therefore, it can be concluded that the overestimation of the effi- 
ciency loss by using the lumped resistance model comes from the over- 
estimation of the voltage drop across the base bulk region, surface sheet 
and fingers resistances. However the lumped resistance model is still a 
good model at low series resistance and low solar concentration. 
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7.6 The Optimum Concentrator Intensity for Some Grid Designs 

ff 

Since the peak efficiency of a high- intensity solar cell is a trade- 
off between good collection efficiency and small series resistance, there 
exists an optimum concentrator intensity for a specific grid design . 

The purpose of this section is to calculate the optimum multi-sun power 
density of some standard grid patterns. The resiilts are summarized in 
Figures?. 1 and 7.4 to 7.7. A summary of the grid designs is shown in 
Table 7.3. 

Figure 7.1 shows one example of the efficiency as a function of 
solar concentration level for the baseline cell. The other baseline 
structures have similar characteristics. In these calculations, the 
grid coverage areas are kept constant while the number of fingers are 
varied. Also two different sheet resistances are used in the calculations 
The optimum solar concentration level is higher for cells with increasing 
numbers of fingers as shown in Table 7.3. 

If the voltage drop across the series resistance can be found 

at the current density of maximum efficiency under the optimum solar 

concentration, a design equation can be derived in a first order model. 

Table 7.4 shows the lumped series resistance, the voltage across R and 

s 

the ratio of the peak -efficiency current density to the short circuit 
current density at the optimum solar concentration for several cell 
models. 

The average voltage across the series resistance is about 46 milli- 
volts under the optimum solar concentration . The ratio of the peak- 
efficiency current density to the short circuit ciirrent is about 0.95. 
Hence the design equation can be given as 


0.046 = 0.95.Ig^.Ng-Rg, 


(7.3) 



Table 7.3. Sunimary of grid design 

for the 

various 

cell 

conf igToratlons . 







a 

a' 

b 

b‘ 

c 

c* 

d 

d’ 

e 

e* 

f 

2 

Cell geometry (cm ) 

1x2 

1x2 

1x2 

1x2 

1x2 

1x2 

1x2 

1x2 

2x2 

2x2 

5 cm 











Dia. 

No. of Fingers 

39 

39 

30 

30 

20 

.20 

10 

10 

60 

60 

240 

Width of Fingers (ym) 

50 

50 

65 

65 

97.5 

97.5 

195 

195 

20 

20 

13‘'^56 i 

M 

Thickness of Fingers (inn) 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

i 

4 ;■ 

Space between fingers (ym) 

462.82 

462,82 

477.82 

477. 

82 510.32 

510.32 

1805 

1805 

313.3 

313.3 

varied * 

Base doping density (Q-cm) 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0 3 H 

it 

Device Thickness Cym) 

250 

250 

250 

250 

250 

250 

250 

250 

250 

250 

250 1; 

Sheet Resistance (fl/lP) 

45 

120 

45 

120 

45 

120 

45 

120 

120 

120 

120 ] 

Width of Bus Bar (ym) 

200 

200 

200 

200 

200 

200 

200 

200 

20-»- 

20^> 

1800 t 



* 







200 

200 


Thickness of Bus Bar (ym) 

3 

3 

3 

3 

3 

3 

3 

3 

3 

15 

4 ;i 

No. of Bus Bars 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

1 ;! 

Concentration Design 

40 

25 

30 

19 

22 

10 

7 

2 

7 

25 

20 


Level 


M 

H 

03 




^le 7.4. The voltage across the lumped series resistor at the optimum solar . concentration . 


Cells 

a 

a’ 

b 

b' 

c 

uuiict 

c' 

snxraxion . 
d 

d* 

lumped series resistance 
(fl) 

0.02748 

0.04087 

0.03303 

0.05566 

0.04999 

0.10090 

0.14164 

0.34524 

voltage drop (mV) 

46.8 

47.4 

45.9 

48.8 

46.3 

46.5 

45.8 

44.7 c 

o . . ^MAX 
Ratxo 

■^sc 

0.957 

0.957 

0.956 

0.953 

0.953 

0.952 

0.955 

0,949 
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120 

where Ig^ is the short circuit current density at one sun power density, 
Ng is the designed solar concentration level and Rg is the lumped series 
resistance of Equation (7.2). A very similar equation has been reported 
elsewhere [2J. . However the numerical interpretation is different 
despite the closeness of both equations. 

Figure 7.4 shows efficiency as a function of illumination 
intensities for the violet-type cells. These cells have two tapered 
common bus tars and 60 fingers evenly spaced across the cell dimensions, 
in this case the one-sun efficiency is much higher than that of the 
baseline cells due to the smaller grid coverage area and the better 
collection efficiency of the violet cells. It also shows the peak 
efficiency can be greatly increased by reducing the conducting resistance 
of the tapered bus bar. For curve e' of Figure 7.4, the violet cell is 
assumed to have a conducting bus bar with a thickness of 15 ym instead 
of 3 ym. 

If the voltage across the common bus bar is plotted against the 
voltage at one terminal of the bus bar, the curve looks like the 
illuminated I-V characteristics. In fact, the voltage difference between 
twO' terminals of the bus bar is proportional to the current density and 
the bus bar resistivity. Figure 7.5 shows these characteristics for the 
simulated violet- t3q)e cells under several solar concentrations . It is 
found that the optimum solar concentration level occurs when the voltage 
across the bus bar has a value of about kT/q. However, the cell efficiency 
is reduced at a higher solar concentration or a larger voltage difference 
across the bus bar^ This value is also a good_paramfeter for the design 
of optimum solar concentration for a violet cell. 




ition. level for a violet type cell. 
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( mV) 



EFFICIENCY (%) 
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Figxire 7.6 gives an example of a solar cell with very low sheet 
resistance, thicker- conducting metal and lower base resistance. The 

optimum solar concentration level is found to be close to 70 AMI solar 

> 

power density-, -'-v 

Figure 7.7 shows a comparison between calculations and the 
experimental values for the large area, high intensity silicon solar 

cell of Sandia Labs Cl] This cell has an effective area of 

2 2 
15.2 cm , and the active area is about 13.38 cm . The finger width is 

56 vim at the bus bar and tapers to 13 ym near the center of the cell. 

The photolithographic defined metallization consists: of approximately- 

4 ym of silver and a thin aluminum underlayer. The predicted optimum 

solar concentration level is about 20 'v 25 which, is very close to the 

measured values. 

7,7” Conclusions 

I 

The effects of series resistances on the multi-sunlight operation ' 
of a solar cell have been examined by a distributive resistor and 
current model. A number of conclusions can be made based on this work. 

(1) An exact numerical calculation can be made to predict the optimum 
illumination Intensity of a specific cell design without fabricating 
the cell. Similarly, this technique can be used to design the optimum 
high intensity solar cell with respect to the. specific solar concentrator. 

(2) High efficiency solar cell operation can be achieved with a low 
sheet resistance design at a high concentration level. A first-order 
equation is also given which can be used to design the optimum concentra- 
tion level at low series resistance and/or low solar concentration levels . 



EFFICIENfcy {•/&) 



Figure 7.7. 


CONCENTRATION ( SUNS ) 


omparison between experiments Ca) and the 
area, high intensity silicon solar cell, 
contact resistances are 10 and 0.03 fi-cm 


two dimensional calculations of the lar 
burface sheet resistance is 120 n/o. Ba 
for models b and c respectively , 
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(3) The lumped resistance model is found to overestimate the power 
losses. This discrepancy is shown to be due to the overestimation of 
the total series resistances of the- surface sheet, finger and base bulk 
in the first order model. 

. ^ 

(4) Ah exact solar cell model provides a valuable analysis technique 

■for deslghing the optimum solar cell under multi- sunlight operation. 
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7*8 Appendix (B). The LuTnped Series Resistance 



V SF 


The lumped series resistance usually consists of five components 
which include finger, surface sheet, base, surface and bottom contact 
resistances. If it is assumed that the space between fingers is much 
less than the length of the fingers and there is negligible voltage 
drop across the bus bar, the lumped finger and sheet resistances can be 
derived as follows. 

Let the voltage drop along a finger be AV(sf) at position y. Then 


y 

AV(y) = / i(LF-y) 
o 


SF 

2 


R, 

1 

WF 

2 


M 


— dy = 


iR,,.SF 2 

“WF ^LF-y 


where i is the current density. The overall average voltage drop across the 
finger beomes 


AV(y) 


Lp iR • SF 

/ AV(y)dy = — S 



I- 


MF.T-A, 

A 



Hence the effective finger resistance is equal to 


Rm'SF 
M 

WF-T-A^ 

A 




similarly, the voltage drop along the sheet resistance is 


AV(x) = / i(^ - x)*R dx=i*R (SF 

0^0 ST j- • X - — ) 

The overall average voltage drop across the sheet resistance becomes 

-rv— 1 Rqp-SF^ 

/ AV. .dx = ^ = I. -Sf, 

- ^ (x) 12 12-A. 

2 " 

„ . R^t’SF 

Hence the lumped sheet resistance is where Rg^ is the surface 

A 

sheet resistance and is the total active area. 
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8. NON-UNIFORM ILLUMINATION EFFECTS ON SOLAR CELLS 

8.1 Abstract 

A computer program has been developed to calculate non-uniform 
illumination effects on a solar cell. The non-uniformity of illumination 
is found to change the conversion efficiency by modifying the resistance 
losses of the top surface layer. It is shown that a non- uniformly 
illuminated solar cell can be operated near its peak power density by a 
suitable design of the top surface sheet resistance and grid pattern. 

8 . 2 Introduction 

• Recently , solar concentrator systems have become popular as a means 
of economically utlizing terrestrial solar energy in the near future. 
However there are some problems which are not experienced in a non- 
concentrated solar system. For example, the input spectral Intensity 
is* usually distorted from the reflection and refraction of the different 
wavelengths in an optical concentrator system. At present, it is not 
clear how the distortion in spectral intensity affects the conversion 
efficiency. It is also unclear as to how intensity variations across 
the concentrated solar cell changes the solar cell characteristics . 

In fact, there are two different kinds of non-uniform illumination 
in the operation of a solar cell. The space— flight cells generally have 
an abrupt light- dark boundary caused by the satellite body or antenna, etc. 
On the other hand, a concentrator cell usually has a steeply varied 
illumination across the solar cell. The degree of non-unifOrmity generally 
increases with the concentration level. 
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A computer program has been developed to calculate the characteristics 
an abruptly illuminated cell with the illuminated boundary pe,rpendicular 
to the finger contact. Two extreme cases of calculation are possible for 
the illuminated area on the near-side or far-side of the bus bar. For a 
concentrator cell, the computer program has been used to calculate two 
non-uniform illuminations in order to determine the effect of a varying 
intensity across the solar cell on the cell conversion efficiency. The 
first model assumes a cosine intensity prqfile with the maximum intensity 
farthest away from the btis bar and zero intensity at the bus bar. The 
second case uses a similar cosine intensity profile but with the maximum 

intensity at the bus bar and zero intensity at the opposite edge of the 
cell. 

This section describes the application of a modified, distributed 
resistance and current model to calculate the characteristics of such 
non-uniformly illuminated cells. In addition, it is shown that the 
efficiency losses due to the non-uniform illumination can be minimized 
by a proper design of the concentrator cell. 

8 #3 Solar Cell Modeling 

A detail two-dimensional model of a solar cell has been reported 
elsewhere (Appendix 10 .4) . 

8.4 Comparison Between Uniform and Non-Uniform Illumination 
The characteristics of a non-uniformly illuminated cell can be 
demonstrated for a 1x2 cm^ baseline solar cell. This cell has a grid 
pattern of 39 fingers and a common bus bar at one edge of the cell. 
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The contact fingers are assumed to be made from silver- titanium 

o 

with approximately 3 ym of silver deposited over a 500 A titanium a^he- 

>V 

sion layer. The 50 ym wide fingers are 0.98 cm- long and evenly spaced 

across the 2 cm cell dimension. This cell has a base resistivity of 

-12 2 

3 S2»cm, and the saturation current density is 5.6x10 A/cm with a 

diode factor of 1.05, The surface sheet resistance is’ 'aBbub'’l5’0 “n/O 

for junction depths of 0,25 ym to 0.3 ym', and the metal contact 

-4 2 

resistance is assumed to have a negligible value of 10 fi»cm . The 

2 

short circuit current density of an AMI solar spectrum is 28.0 mA/cm 
in the active area. 

The cell efficiency as a function of solar concentration for the 
general non-uniform illumination is shown in Figure 8.1. The cosine 
illumination, which has a maximum intensity near the bus bar, produces 
the highest efficiency due to the lowest power losses in the top layer 
and contact grid resistance. On the contrary, the cosine illumination 
which has a maximum intensity farthest from the bus bar has the lowest 
efficiency at all illumination levels . The uniformity of the illumination 
influences the cell efficiency by modifying the resistance power losses 
in the top layer and grid resistance. Therefore, the efficiency of a 
solar cell depends on the current distribution on the cell surface due 
to a non-uniform illumination. 

In Figure 8.1, a calculation is also shown for the same baseline 
cell which has a grid thickness of 6 ym. It is interesting to note 
that the uniform illumination produces a higher efficiency at high 
concentration than the case of maximum illumination near the grid contact . 
This is because the sheet resistance losses counterbalance the reduced 
grid resistance losses in such a favorable cosine non-uniform illumination 
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More evidence of the competition between grid resistance losses and 
surface sheet resistance losses can be obtained from the comparison of 
the two cases in Figure 8.1. The power losses in the top surface 
layer consist of the two components of sheet resistance and grid 
resistance losses. The power losses in the surface sheet resistance 
dominate over that in the grid resistance at high intensities of 
illumination. Therefore the total power loss of a non-uniformly 
illuminated cell is higher than that of a uniformly illuminated cell, 
despite the smaller grid resistance loss in this case. From Figure 8.1, 
it is clear that the calculation of a non-uniformly illuminated cell 
can be used to show the dispersion of the conversion efficiency as a 
function of the oonoerfetation level and the non-uniformity of illumination. 

Figure 8.2 shows a solar cell design for multi-sun operation. The 
effect of non-uniformity of illumination is shown to be reduced to a 
negligible minimum with a design of low surface sheet resistance. 

Figure 8.3 shows the results of the efficiency as a function of the 
active area for two cases of abrupt illumination. The results are 
plotted as a ratio of efficiency versus the ratio of the active area, 
where the ratios of efficiency and active area are referred to the fully 
Illuminated solar cell. It is clear that the cell efficiency is nearly 
proportional to the active area and is not a strong function of the 
Illuminating profile. This is because the cell efficiency is limited 
for low intensities by the collection efficiency instead of the power 
losses in the senes resistance. The power losses in the grid resistance 
is small and the conversion efficiency is only a function of the active 


area. 
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RATIO OF EFFICIENCY 


0.25 



05 0.75 |.( 

RATIO OF ACTIVE AREA 

Figure 8.3 Ratio of the conversion efficiency as a function of the 
ratio of the illuminated area for the abrupt illuminated 
cell. (0) dark area near bus bar; (D), (X) light area near 
bus bar. The surface sheet resistance is assumed a value 
of 500 fi/O for (X) points instead of 120 ^/a for (0) and 

(o) 
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o,b Conclusions 

Some of the oheraoteristios of non-uniform illumination on a solar 

ceu have been analysed in this ohapter. Several important conolusions 
can be made as follows: 

(1) The conversion efficiency can be analysed as a function of the 
ceu structure, grid pattern. iUumination level and the non-uniformity 
of Illumination by the two-dimensional solar cell model. This phase of 
work is important in the design of solar cells under concentration and 
in the design of concentration systems. 

(2) The power losses in the surface sheet resistance, bulk resistano 

contact resistance and grid resistanno k 

g a resistance can be accurately predicted with 

the model developed in this work. 

(3) The effects of non-uniform iUumination on cell conversion 

efficiency have been analysed for a few solar ceU designs. It has also 

been shown that the effects of non-unifon. illumination can be reduced 

tc a negligible amount by an appropriate Resign of the surface and grid 
resistances. 

(4) Perhaps, the most important conclusion is that the non-uniform 
illumination is found to decrease the cell efficiency at high concentration 
levels irrespective of the lliuminatlon profile, if the surface sheet 
resistance is the dominant loss factor. 
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9. HIGH-LOW JUNCTION EMITTER SOLAR CELLS 
9.1 Abstract 

This section discusses the physical characteristics of a recently 
proposed solar cell — the HLE solar cell [1] which has been predicted to 
give a substantial increase in the output voltage and the converstion 
efficiency of a highly doped junction solar cell. However, our calcula- 
tion predicts the negative results. The discrepancies have also been 
Identified as high injection effects and heavy doping effects in the 
emitter-high- low junction of the HLE solar cell. 

9.2 Intr oduc t ion 

High efficiency has been predicted for solar cells with low base 
resistivity near about 0.1^. cm [2]. Unfortunately, the measured conver- 
sion efficiency is substantially less than the theoretical expectation 
for these highly doped cells. This has been found to be due to the 
discrepancy in the values of the output voltage and open circuit 
voltage. In fact, the measured value of open circuit voltage is about' 
150 tnV less than the predicted value from the simple .Shockley diffusion 
theory for a 0.1 n* cm cell. 

Physical studie,s have recently shown that the discrepancy 
of the performance of a highly doped solar cell is due to heavy 
doping effects in the heavily doped surface layer [3]. Therefore, there 
is excess surface recombination current which is usually much larger than 
the base layer injection current density for a highly doped solar cell. 
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In order to reduce the abnormally large back injection current 
density, a high-low junction emitter solar cell has been proposed and 
predicted to achieve high following a similar reason of a high-low 
junction base solar cell [1]. 

This section will discuss the physical operation of a high-low 
•emitter sili.con solar cell. The results of the numerical calculation 
will be presented for several structures of the proposed high-low 
emitter cells. 

■ 9.3 Numerical Calculations 

A HLE cell has an additional high— low junction located in the sur- 
face region of a conventional cell. Figure 9.1 shows a schematic struc- 
ture of a HLE cell with a structure of n‘*'np. The back injection current 



density is expected to be reduced from the minority carrier confinement 
in the surface n-type region due to the small SRV at the emitter high-low 
junction. The open circuit voltage will thereupon be increased from the 
resultant smaller junction saturation current density. Before discussing 
the physical fundamentals of the emitter hlgh.-low junction cell, the 
numerical calculations will be presented for several proposed 
models of Table 9.1. 



Table 9.1. Structures of 
solar cells. X , and 
200 Pm, resnectivelv. 

models a to 
Xg are equal 



Ng (*/cm^) 

Ng ( /cm^) 

Model a 


2x10^^ 

10^^ 

Model b 


2x10^^ 

10^^ 

Model c 


2x10^° 

10^^ 

Model d 


2x10^° 

lo"-^ 

Model e 


2x10^® 

10^^ 

Model f 


2x10^^ 

10^"^ 

Model g 


2x10^^ 

10^^ 

Model h 


2x10^® 

10^6 

Model 1 


2x10^® 

■ 10^® 

Model j 


2x10^*^ 

10^' 


(pm) 

0 

1 

10 

0 

1 

10 

1 

5 

5 

10 
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Efficiency 

16.59 
14.77 
12.23 
16.25 
14.58 

11.60 
15.40 
14.94 
13.70 
12.78 


Cells a-j are assumed to have a 595 Af&lokness of Ta^Oj at one AMO 
aclar Intensity. The SRV of the diffused surfaee Is assumed a value of 

10^ cm/sec. The base diffusion length is assumed to be 100 ym, which is 

close to L(MAX) at a base resistivitv of 0 1 rrir 

y cm. The surface diffusion 

length is taken to be (L(MED) + L(MIN))/2. Therefore, the minority car- 
riers of the surface epitaxial layer have been assumed a diffusion length 
of 55 um. Ihxs value is very close to the experimental value at an 
n-type doping density of 10^^/cm^ [4]. 

The results of the numerical calculations are summarised in Figures 
9.2 and 9.3. It is seen in Figure 9.2 that the cell conversion 
efficiency is lower at a wider width of the constant doped n-layer. This 
la a direct consequence of the high injection effects in the lightly 
doped n-layer. In Figure 9.3, it is clear that the cell with a wider 
n-layer has a high slope in the dark I-V characteristics. Therefore, the 
cell curve factor is degraded and the conversion efficiency is lower. At 
the same time, the back injection current density is larger for cells 
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Figure 9. 
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VOLTRGE (VOLTS) 

Figure 9.3. Dark I-V characteristics of n'^np cell c, h, j, and k of 
Table 9.1 ’ 




4. ‘to 

with a wider epitaxial layer. This is due to the ineffectiveness of the 
SRV at the EHL junction and the low carrier lifetime or diffusion length 
of the n-type material. The back injection current density can be 


conveniently represented as 


f ^Aqpdx 


where the first tern, is the recombination current density in the con- 

etantly doped layer. It is clear that the back injection current density 

is higher for a wider epitaxial layer of if the second component of 
Equation 9.1 can be neglected. 


9.4 Physical Mechanisms of an Emitter High-Low Junction Solar Cell 
It has been established that the carrier transport through a conven 
tlonal p-n junction can be described by the Shockley diffusion theory or 
Sah's recombination current model. However, a high-low junction can be 
described by the HL junction theory or the HL junction theory plus the 
junction leakage current model [5.6]. In these models, a high-low junc- 
tion has the advantage of confining the minority carriers in the lightly 
doped region where the carrier lifetime is supposed to be higher than 
that of the highly doped layer. Therefore, the junction saturation cur- 
rent density Is lower and the corresponding will be higher than Its 
conventional counterpart. 

In order to determine the validity of the high-low junction model. 

a comparison between experiment and model is attempted and shown in 

Figure 9.5 for a high-low junction base solar cell of Figure 9.4 and 
Table 9.2 [7]. f 



Table 9.2i Parameters used in the calculation of V variation with 
substrate resistivity. 


Substrate 


Resistivity 

SJ’cm 

0.01 

0.1 

1.0 

10 


Resistivity 
* cm 

10 

10 

10 

10 


Diffusion ^ 

Diffusion Length 

Width 

Coeff.Ccm /sec) 

(wm) 

(urn) 

3.2 

10 

250 

9.3 

50 

250 

22 

)) ■ 125 

250 

33 

C\ 150 

250 

Epitaxial Layer 


Diffusion 

Diffusion Length 

Thickness 

Coeff.(cm /sec) 

( m) 

(Um) 

33 

80 

10 

33 

80 

10 

33 

80 _ 

10 

33 

80 

10 




The calculations of Figure 9.5 are based on the following eguatlons 
for Shockley's model and HL junction model: 
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(9.4) 


S =: 


V V ^ ^ ^ 


L. + 


P 


s a + Jg-- %■ 

o H 


'3* (9.5) 


P 


J,iL Jhll Mturatlon current densities for ShccUey's 

"■°del. HL Junction , nodal and HL junction plus leakage current n,odel, 

p Ivoly. s and are the normalized and unnormalized SRV at Che 
HL junction, respectively. These values are functions of the material 
parameters on both sides of the HL junction, while and N„+ are the 
majority carrier densities on both sides of the HL junction. 

It Is clear from the calculation of Figure 9.5 that the high-low 
Junction theory is a good model for a high-low junction cell. The con- 
finement of the minority carriers In the lightly doped Wer requires a 
hlShly reflecting HL surface or a low SRV at the HL junction. This can 
be easily achieved In a high resistivity BHL junction solar cell. On the 
other hand, this may not be true for a EHL junction solar cell, since 
the EHL junction Is located near the Illuminated surface. will then 
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Figure 9.5. 


exp 
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be increased by the optical generated majority carrier density in the 
lightly doped n-type region. This is shown in Table 9.4. Therefore, the 
carrier confinement in the n-type region is deteriorated. On the other 
hand, for a BHL solar cell, the interaction between the base HL junction 
and the optical generated carriers is negligible due to the depth of the 
BHL junction from the illuminated surface. In this case, S is pri- 
marlly a function of the material parameters at both sides of the HL 
junction. 

The predictions of from the first order model are shown as Tables 
9.3 and 9.4 for cell (c) of Table 9.1. These are compared to the value 
of 0.565 volts from the exact numerical calculation. 

“I" 

9,5 P PN Emitter High-Low Junction Solar Cell 
Since the previously proposed n np solar cell is not more efficient 
than the conventional junction solar cell, we now propose a similar EHL 
junction solar cell with a p^n structure. The new models are shown in 
Table 9.5. The structures of these cells are similar to those shown in 
Figure 9.1. In the numerical calculations, the solar cells are assumed 

O 

to have a 595 A thickness of Ta^O^ at one AMO solar intensity. The SRV 
of the top surface is assumed to have a value of 10^ cm/sec. The base 
diffusion lengths are assumed to be 30 ym and 60 ym, which are close to 
Ljj(MED) for the n-type base doping densities of 5x10^^ /cm^ and 10^^ /cit? 
respectively. The surface p"^ and p layers are assumed to have a diffu- 
sion length of (Lj^CMAX) + Lp(MED))/2. 

The results of the numerical calculations are shown in Figure 9.6, 
which shows the dark I-V characteristics for cells of Table 9.5. The 
dark current densities are found to be nearly equal at high forward volt- 
age irrespective of the doping density in the p-type layer. This implies 
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Table 9.3. Parameters used in the calculation of V for cell (c) of 
Table 9.1. 


EPITAXIAL LAYER 



Doping 

Diff . 

Diff. Length 

Width 


Density 

Coeff . 

(yra) 

(ym) 

Shockly's model 

o 

r-i 

12 

55 

10 

HL junction model 

10^'' 

12 

55 

10 

HL junction plus 
leakage current 
model 


12 " 

SURFACE LAYER 

55 

10 

• 

Doping 

Diff. 

Diff. Length 

Thickness 


Density 

Coeff. 

(V«n) 

(ym) 

Shockly's model 

10^^ 

2 

2.5 

0.25 

HL junction model 

10^^ 

2 

2.5 

0.25 

HL junction plus 
leakage current 
model 

10^^ 

2 

2.5 

0.25 


Table 9.4. Calculated from the first order model. Model A uses the 
extrinsic doping density of the n-type region in the calcula- 
tion of Vqc while Model B takes into account the light- 
generated carrier density in n-type region. 


SRV 

MODEL A 
at HL junction 
(cm/sec) 

, Voc 
(volts) 

MODEL B 

SRV at HL junction 
(cm/sec) 

Voc 

Cvolts) 


S 

o 

S 


'o S 


Shockly's Model 

- 

- 


— _ 

0.416 

HL junction Model 

0.80 

(3,41) 

0.504 

8.0 1.29 

0.563 

HL junction plus 
leakage current 
Model 

0.80 

(5.7) 

0.503 

8.0 12.8 

0.562 


In Table 9.4, model B takes into account the light generated carrier 
density which is equal to Gt. The calculated value of the light generated 
carrier density is about 10 /cm where G has a value of about 10^^/cm3 
in the n-type region. 



Table 9.5. Structures of models a' to f- of p+p„ junction solar cell ' 
rispe§tt^l"^ 200y». 


Model 

a' 

b' 

c’ 

d' 

e’ 

f ! 

Ng at/cmh 

5x10^^ 

5x10^^ 

5x10^^ 

5x10^^ 

5x10^^ 

r 

,10^^ 

Ng (///cm^) 

- 

10^^ 

10« 

io“ 

5x10^^ 

10^' 

Xg (#/cm^) 

0 

10 

20 

10 

10 

10 

Efficiency (X) 

14.53 

15.76 

15.84 

15.66 

15.00 

15 . 30 


Table 9.6. The photovoltaic characteristics of emitter hleh lo« • 

solar cell of Table 9.1 and 9.5. nigh-low junction 



^ h ■ i b> 

^ ucture _n^ Q-^np n^p p+pn n-*-no 

Ig^CMi/Zcm^ 79.82 77.94 62.62 75.21 74.92 

^oc 0*565 0.630 0.660 0.678 0.678 

Ij^^/2cm^) 73.84 74.19 59.84 72.39 72.09 

Pj^(MW/2cm ) 33.08 40.43 34.59 42.64 42.39 

OFF 0.733 0.823 0.837 0.836 0.834 

EFF (%) 12.22 14.94 12.78 15.76 15.66 


-P. Pn 
71.24 
0.678 
68.20 
40.58 
0.841 
15,00 


+ + + 
P pn n npp 


77.70 85.16 
0.649 0.566 
73.03 80.25 
41.41 37.15 


0.821 0.771 




CURRENT DENSITY Cfi/CM'*) 




^ wju 


that the suppression of the dark emitter recombination current density 
and the total saturation current density m primarily determined by the 
base recombination current density In a pV solar cell. Therefore, the 
open circuit voltage of these cells are relatively constant as shown In 

Table 9.6. which contains the detailed photovoltaic characteristics for 
some cells of Tables 9.1 and 9.5. 

Finally, a novel n\pp+ solar cell Is modelled which has both 
emitter and base high-low Junctions. Both n-type and p-type layers are 
assumed to have a doping density of I.IxIqI^ /cml The p-type base layer 
has a diffusion length of 300mm and the n-type surface Is assumed to have 
a diffusion length of (l^(med) + lj,(MM))/2. The emitter and base high- 
low junctions have doping densities of IxlO^O /cm^ and 10« /cm^ and 

depths of 0.25:pm and 0.5pm. respectively. The calculated results are 

shown as Model k in Tablp Q t*- , 

Table 9.6. It is found that this cell does not show 

higher efficiency despite Its high short circuit current density. 

9.6 Discussions and Conclusions 
This section will discuss the validity of the calculation of V 
rom the p-n junction saturation current' density and the corresponding 
Short circuit current density. In fact. Is a series combination of 
the voltages across the p-n junction, surface region, base region and HL 
junction at the open circuit condition, among the above components, the 
Dember voltages across surface and base regions are usually very small. 

The voltage across the HL junction Is also small unless under high Inten- 
sity illumination where reaches the junction built-in voltage [8] . 
Therefore. It Is a good approximation to calculate from the junction 
saturation density and short circuit current density for a solar cell 
under one AMO solar density. 
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It can be concluded that the proposed n^np junction solar cells do 
not yield a higher output voltage as expected from the incorporation of 
an additional high resistivity emitter layer in a conventional solar 
cell. The results show that the highest open circuit voltages are 
obtained for low resistivity emitter layers. On the other hand, a p pn 
junction solar cell can produce a higher open circuit voltage than an 
n^np cell irrespective of the doping density in the surface layer. This 
is due to a more effective high-low junction and a higher diffusion 
length for the p^p high-low junction. 
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10. APPENDICES 








Appendix 10.1 


CpMPARISO^f OF THEORETIC Ar. AND 

expekim ental sodar cell performmce* ’* 

c. R. Fang and J. R. Hauser 


ABSTRACT 


North. Carolina 
Raleigh, 


This paper presents the results of a detailed 
comparison between the theoretically predioS^ 
perfo^ance of silicon solar cells Ld the e^eri- 
performance. The comparisoS has 
ncluded the dark I-V characteristics, the soec- 
tral response characteristics and the I-V chLac- 
teristics under AMO illumination. In generafit 
as een found that the agreement between theory 
and experimental behavior is very good. This TOod 
agreement has been obtained for a variety of cEK 
fabricated with junction depths ranging ft«m 0 1 

pm to more than 1 pm and for base ifyer Sstivi- 
ties from 10 R* cm to 0,1 R.cm. resistivi 


INTRODUCTION 

interest has increased in solar ceUs in 
recent years, considerable progress has bepn made 

sol2°rn understanding the physics underlying 
XSe Sories^S°"' discrepancies betweL ^ 

Sfanf stud-f n f in?er- 

Study [1] . The presence of high dooine 

band pp reduction effects has generally beLmf 
accepted as the reason for the Iqw nnon ■ 

llth" blsM^y^r^^^ola^^cSl. 

this effect^CS]^'^"’® =onb^^oversy over 

ho„ rTl!! "as undertaken to determine 

i?th ®srees with experiS? 

known performance. It is 

known that theoretical predictions are in general 
agreement with experimental results. However no 
- tai_ed comparison has been reported between’ 

ANALYTICAL APPROACH AND DEVICE PARAMETERS 

has been performed using a 
solution of the semiconductor 

£Snrn/°h applied to solar cells, 

stalls of uhe approach and analysis have been 

useu is quite free from the normal simplifying 
pproximations made in semiconductor device 

^is section has. been' published in 
the Proceedings of ' the Thirteenth 
Photovoltaic Specialists Conference. 


state University 
NC 27650 

SSS™ Including th. effects on„S- 

tlcn 2rL£n^%“ff"S?-Sett^^ 

generation rate and the device £eLtio£i T iLr, 
pertaining to collection-effi- ® 

Recombination is included within the 

end tne Sffe^S=lf.“" SS ’ 

in e impurity profile is included 

(g) 'A non-ohmio contact is included at thp 

b£n£So£1e!oS;! ^ --- 

tance in thi £Slinrof™solaJ Tellf. "JhrniSr'" 
important of these are reviewed here to dL£L 
the values and models used in the analysis. The 
diffusion length is perhaps the single 
device parameter. The diffusion 
5° decrease in general with 
incr easing impur>ity concentration. However widelv 
rying values as shown in Figure 1 have been ^ 
diffusion length It Zy gj£n 

(Sw Sl7f^3 doping LfsitLs 

^below 10 /cm3) measured diffusion lengths rouom 

order! of magni- 

£e sca£e£ift£^“® densities (lol9/cm3 or Sove) 
ine scatter in the experimental values aDDears i-n 
be much less as seen in Figure 1, This i! liLw 
to to the do.i„„c. of Auger recombtotto 
large doping densities. 



For solar cells which have good efficiency 
values, the base layer diffusion length must be 
above the center curve labeled L_(MED) in Figure 1. . 
Diffusion length values about midway between the ; 

and LACKED) curves of Figure 1 are typical 
"bF the val,ues^for the cells studied in this work. 

For the base layer of the cells studied here, 
measured diffusion lengths were available. In most 
oases it was found that the calculated solar cell 
properties agreed Well with the experimental data 
when the measured diffusion length or a value close 
,to the measured value was used in- the analysis. 

' -fhe BSF cells studied were found to be an exception 
to this rule. For these cells good agreement 
between theory and experiment could only be obtained 
using diffusing length values considerably larger 
than those measured by the X-ray technique. This is 
discussed in more detail in a later section. 

For the diffused surface layer experimental 
data was not available for the diffusion length so 
the data of Figure 1 was used as a starting point 
in the analysis. The surface layer diffusion \ 
length was adjusted to obtain the best agreement 
between theory and experiment in the short' wave- 
length region of the spectral response. Specific 
values used in the analysis are discussed in connec- 
tion with the results.. In general, the best agree- 
ment between theory and experiment was obtained 
when the surface diffusion length was taken some- 
where between the LACKED) ai^d Lp(MIN) curves of 
Figure 1. Since the surface layer is heavily 
doped, the region of Figure 1 of importance for the 
siupface is around lO^^/cra^ or abova* In this region 
the difference between the solid cxu'-Ves is not 
nearly as great as at lower doping densities. 

The properties of the surface layer of any 
silicon solar cell plays ^an important part in the 
performance of the cell. In this work, the analy- 
sis has been found to be sensitive to the doping 
profile used in the surface layer. 'As an approxi- 
mation a Gaussian or erfo profile is frequently 
used for a diffused layer. Experimentally it is 
known that a phosphorous diffused layer has an 
impurity profile which differs significantly from 
the simple Gaussian or erfo profile [9]. Near the' 
surface a region of near.l.y constant doping appears. 
This is followed by a rapid drop in carrier concen- 
tration to a region deeper into the diffused layer 
which follows fairly closely an erfc profile Tsai 
L9] has given curves from which the surface doping 
density C , the width of the constantly doped layer 
X and thi doping density Cg at the boundary of the 
constantly doped and erfo doped region can be esti- 
mated using the diffusion time and temperature. 

In the early stages of this work simple erfo 
or Gaussian doping profiles were used in the theo- 
retical calculations. With such profiles it was 
found to be very difficult to obtain good agreement 
between theory and experiment in the spectral 
response data. In most cases good agreement was 
never obtained with the simple doping profiles. 

After an impurity profile such as that found by Tsai 
[9] was used in the calculations, it became much 


easier to obtain goqd agreement between theory and 
experiment. 

Not only the doping profile, but also heavy 
(doping bandgap reduction effects appear to be 
important in the diffused surface layer. It has 
now been pointed out by several workers that the 
bandgap reduction effect appears to be the physical 
origin of the low open circuit voltage of low 
resistivity base layer cells [2,10J. As discussed 
previously, bandgap reduction effects were included 
in the analysis performed here and appears to be 
essential for obtaining good agreement between 
theory and experiment unless one uses very small 
diffusion lengths within the surface layer. The 
combination of the correct doping profile and 
bandgap reduction effects leads to a thin stirface 
layer whii-h behaves almost identically for carrier 
collection as the "dead layer" model of Lindmayer 
[ 11 ]. 

The Keldysh-Franz effect has recently been 
proposed as an important effect in solar cells 
[3]. This effect was included in some of the 
theoretical calculations using experimental data 
for the bandgap reduction due to an electric field 
[12] . The major effect has been to increase the 
calculated dark current from a cell at low volt- 
ages by a factor of 2 to 5. The effect has typi- 
cally been very small for voltages above about 0.5 
volts forwtrd bias. The Keldysh-Franz effect can 
only be of major importance in solar cells domina- 
ted by depletion region current and the cells 
studied in this work did not show this character- 
istic. 

An accurate modeling of the optical processes 
is important to obtaining good agreement between 
theory and experiment in spectral response data. 

In general the optical processes associated with 
surface reflection and bulk absorption are well 
known and characterized for silicon. One region 
where some uncertainty still exists is the long 
wavelength absorption coefficient around 1.0 pm. 
This is difficult to experimentally measure 
because of lopg wavelength free carrier absorption 
in silicon. At the beginning of this work the 
data of Dash and Newman [13] was used for o in the 
long wavelength region. However, it was soon 
foimd that good agreement between theory and 
experiment could not be obtained at 0.95 pm and 
1.0 pm using this data. The reported range of 
absorption coefficient values in this wavelength 
'range is shown in Table 1. As can be seen there 
are very significant reported differences with 
Dash and Newman's values being near the largest 
reported. After using several potential values, 
good agreement between theory and experiment was 
obtained using the values in the last row of 
Table 1. These values can be seen to be well with- 
in the range of reported values and we believe 
these values are more accurate than most of the 
reported experimental values since they give good 
agreement b'etvfeen theory and experiment for a wide 
variety of solar cells. 

In'iiaddition to absorption coefficient an 
accurate modeling of surface reflection is needed. 




Soma cff the cells studied 
while others had a Ta_0_ 
thickness. At the begiiln 
son was made bet^wan theory and experimental- data' 
tor the transmission coefficient of a bare surface 
and a Ta„0 coated surface. The agreement was 
within aBovSt 1% except in the 0,9 to 1.1 ym region 
where larger differences occvirrad. The origin of 
the discrepancy here is not known. However, the 
excellent agreement elsewhere, give confidence in 
the accurate modeling of surface reflection and 
transmission properties . 


had bare silicon surfaces 
coating of about S95A in 
ing of this work a eomnar»^-. 


A complete characterisation of solar colls 
requires that account be token of the two-dimen- 
associated with the contact grids. 
Spectral response measurements made at low light 
levels are not greatly affected by the contact 
grid sxoept for the area reduction factor". Ho\i- 
and light I-V characteristics are 
influenced by the sheet resistance, contact 

two-dimensional pronerties. Values 

infSnSrs^ factor are especially 

intluencod by t\vo-dimensional effects. In this 

Srid moddl was used to 

accwataly model these affects, Details of the 
analysis and model are discussed olsevvhera [ 19 ]. 

A tw-dimensional modeling of the contact effects 
IS necess^y if good agreement is to be ob?Sned 
between theory and experiment for effieiancy, 
open circuit voltage, and curve factor. 


RESULTS OF COMPARISON 

A few general comments are needed to exolain 
the manner in which the comparison between theo- ' 
retinal and experimental results has been made. 
Tie experimental data which was used consists of 
three different typos of measurements. These are 
1) dark current 1-V ohoraotaristics, 2) short 
circuit spectral response data and 3) light I-V 
characteristics under simulated AMO, 1 sun illu- 

ssneral breakdown, 

u-io tpes of doi'k current data was available. 

T V simply measures the terminal 

IV characteristics of a solar call in the dark, 
me second technique measures ooen circuit volt- 
circuit current I^„ at various 
intensities of illumination and plo?l the result- 
ing pairs of numbers, Those two techniques will 
be referred to as the dark 1-V technique and the 

ideal ona-dimenaional 


^ cell with no series resistance both teoh- 
niquas should give the same curve, since at the 
open circuit volta;’ ' condition the internal dark 
currant exactly equals the short circuit solar 
generated ourront. The two curves differ whan 
senes resistance and t^-fc-dimonsional effects 
are considered , At high ourrent densities the 
.sq'^oq icchniquq tends to be much closer to the 
laeai ona-dimonsional theory since the voltage 
IS measured under conditions of saro net current 
a<5lar cell, The differences betaraan 
the t^-(o measurement techniques will become clearer 
after the discussion of the specific examples. 

theory. and experimental data 
n t all device parameters are of major importance 


^ each sat of measurements . For example the spec- 
tral response depends very little on sheet resis- 
tance, contact resistance and other two-dimensional 
properties of the coll. Also as just discussed the 
^sc“^oo "’a^suremant is relatively insensitive to 
tn-io-dimenslonal effects. The device parameters of 
Mst importance in these characteristics are the 
device diffusion lengths, doping profile, surface 
recombination velocity, and heavy doping models. 
Thus the data and the spectral response 

data has been compared with the results of the enu- 
dimensional calculations. The parameters of sheet 
resistance, contact resistance and grid structure 
were than included along with the one-dimensional 
results in the two-dimensional computer program an-^i 
compared with the dark X-V characteristics and the 
1-V characteristics. For almost all of the 
devBaes' a combination of sheet resistance and con- 
tact res-istance was- required to accurately describu 
the t\-fo-dimansional nature of the calls. 


^ The solar colls studied can be broadly divided 
into ttoee categories of cells based upon the 
resistivity of the base layer and the type of cell 
design. The first category consists of n‘*'-p cells 
with 0.1 n- cm base layers and a finished thickness 
of either 6 mils or 10.2 mils. The second categorv 
consists of n -p colls with 10 fi*cm base layers and 
a thiCKTiQSs of about 10. S milsi The third typo of 
cell is a back surface field (BSF) cell with a 
16 n*cm base layer and a thickness of about 6.S mils. 


The first ti-fo categories of cells were 1 cm hv 
2 cm and used a ten finger grid pattern wlulo the 
third category of calls were 1 cm by 2 cm and used 
a nine finger grid pattern. Table 2 gives parame- 
ters for the t^-Jo different gird patterns used on 
the cells. The first two categories had a bare 
silicon surface while the third category cells were 
coated with S9S A of Ta 0. and a S mil Teflon FEP 
cover . * ® 

The solar cells ware fabricated and measured at 
the NASA-Lowis_Rosoaroh Center, Cleveland, Ohio. 

The n-typo surface layex-s wore phosphorous diffused, 
using POGlg, The tamparaturo and duration of the 
diffusion process are described for each type of 
call in the following sections , The top and bottom 
contacts wore made using Sq metal mask with evapo- 
rated aluminum (200 "v 500 A) followed by evaporated 
silver of a few pm, The contacts ware sintered at 
temperatures of 550 to 6So°G in For each cell 
the base layer diffusion length wSs measured at 
NASA Lewis by the N-ray method. 

Within each class of colls the measured char- 
acteristics were consistent and certain cells con- 
sidered to be typical of each type were selected 
for detailed study. The major differences in any 
given category of cells were in the dark current 
-values at low voltages. These differences do not 
greatly influence solar cell performance since 
operation is at larger voltages where the differ- 
ences between cells is reduced. 

CELLS WITH 0.1 0* CM BASE LAYERS 

Results for two different colls within this 
general category are presented. These cells were 


15 V 



fabricated on 0.1 n*cit Boron doped 

wafers. Various device parameters for the two 
cells which are D-1 and D-2 are listed in Table 3. 
The major differences between the two cells are the 
cell thickness, 6 mils and 10.2 rails, and junction 
depth, 1.0 pm and 0,72 pm. The base layer diffu- 
sion lengths of 60 pm and 105 ym are the experi- 
mental values which were also used in the compu- 
ter calculations , 

A comparison of the calculated and experi- 
mental dark I-V characteristics is shown in Figure 
2. The comparison between theory and experiment is 
very good for voltage values of 0.45 volts or 
larger. There is a large difference between 
theory and experiment at low voltages and this is 
typical of all the solar cells studied. It might 
first be suggested that this difference is due to 
the Franz-Keldish effect which was neglected in 
the calculations. However the current is much 
larger than that calculated including the Franz- 
Keldish affect. Also the voltage dependence of 
this excess current is not ooniPistent with the 
Franz-Keldish effect. _ | 

The triangle points in Fig^vre 2 show the cur- 
rent which would result from a pUJie resistor of 
2.57 Kfl in parallel with the solarj^cell. As can 
be seen from the figure this voltage dependence 
fita the excess current almoat exactly. The pres- 
ence of such a shunting resistance has been seen 
in all cells studied. The magnitude of this resis- 
tance has been observed to vary greatly from cell 
to cell within a particular category. This shunt- 
ing resistance has not been studied in any detail 
in this work, but it has also been observed by 
previous workers [20]. The physical origin of this 
resistance remains somewhat of a mystery. It may 
be due to shunting precipitates for example with- 
in the junction depletion region. 

The dotted cur* in Figure 2 illustrate the 
sensitivity of the calculated dark I-V character- 
istic to the diffusion length (or lifetime) in 
the diffused surface layer. Values calculated 
using the Lj^(MED) curve of Figure 1 are definitely 
below the experimental curve. The best fit was 
obtained using values midway beween b (MED) and 
Ljj(MIN) of Figure 1. The final selection of base 
layer and surface layer dif^nsion lengths must be 
made from comparing not only the dark current data 
but also the spectral response data. A small sur-- 
•face layer lifetime cannot be distinguished from 
a small base layer lifetime using dark current 
data alone since both lead to large dark currents. 
However, differences Can be seen between these teo 
cases in the spectral response calculations. 

A comparison bat\'feen the theoretical and 
experimental spectral response calculations is 
shown in Figure 3 . The parameters wliich give good 
dark current calculations are also seen to provide 
good spectral response calculations. Curves are ■ 
also shotm of the relative contributions to the 
spectral response from the base layer, depletion 
layer and surface layer. The spectral response 
calculations prove to be a sensitive test for 
accurate device parameters. 


The spectral response at 0.9 - 1.0 pm is 
determined almost entirely by the base layer 
properties. Values at this point are very sensi- 
tive to the base layer diffusion length and long 
■wavelength absorption coefficient. The corrections 
to the long wavelength absorption coefficient dis- 
cussed earlier ware identified studying the long 
wavelength spectral response. Before the absorp- 
tion coefficient values of the last line in Table 
•1 were used, it was not possible to accurately 
match the dark current data and long wavelength 
spectral response, not only of this cell, but of 
all the cells studied. 

The spectral response at 0.4 - 0.45 pm is 
determined almost entirely by the surface layer 
properties, which include not only surface layer 
lifetime, but also surface recombination velooi-ty, 
doping profile and bandgap reduction effects. 
Because of the interaction of all these effects 
it is difficult to aftribute a given spectral 
response value to any one of these effects alone. 
Initial calculations were made for the cells 
using an erfc’ doping profile in the surface 
layer. With this type of doping profile it 
l^roved difficult to obtain a good match of spec- 
tral response at all wavelengths . In general the 
spectral response would be too large at 0.4 - 
0.45 pm if a good fit was obtained at larger wave- 
lenths . Reasonable combinations of diffusion 
length, sui-face recombination, surface doping ate. 
could not bo found to match the data. However, 
when an impurity pr\3file modeled after the experi- 
mental work of Tsai was used, good agreement such 
as shown in Figure 3 became relatively easy to 
obtain. The parametet's describing the surface 
doping profile are listed in Table 4 for this 
device and the other devices studied. 

In the final analysis the only device param- 
eters which were adjusted somewhat to obtain the 
good agreement beti'feen theory and experiment seen 
in Figures 2 and 3 were the surface recombination 
velocity and the surface layer diffusion length. 

The final values used and reported in Table 3 are 
.quite consistent with the diffusion length data 
,of Figure 1 and other reported values of surface 
recombination velocity [7]. 

In calculating the light I-V characteristics 
•and the overall efficiency at AMO , additional 
parameters describing the contact grid (Table 2) 
and values of surface layer sheet resistance and 
contact resistance must be introduced into the 
tiwo-dimensional model. The D-1 solar cell and 
other devices with the 1.0 pm junction depth 
were foxind to have very low values of sheet 
resistance and contact resistance . Values used 
in the analysis are given in Table 5 for this and 
the other cells studied. Figure 4 shows a com- 
parison Between the solar cell power quadrant I-V 
characteristics and the theoretical calculations . 
Again the agreement between theory and experiment 
is very good . 

The second cell in this category for which^ 
data is reported (D-2) is similar to that just 
discussed except for the shallower junction depth 
and, increased base layer thickness as shown in 
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^ show that essentially 

the s me parameters ware used in the theoretical 

^ diffusion 

length of 105 jira was measured in this call and this 
value was used in the calculations. 

Figures 5, 6 and 7 show a comparison bet^^een 
theory and experiment for this call. In Figure 5, 

experiment 

^ dmk I-V and !„ -V measurements at voltages 
.of 0.6 and 0.7 volli. °?he I -V data agrees 
one-dimensionai°ca£oulations while 
the wo-dimansional calculations including sheet 
resistance and contact resistance is required to 
describe the dark 1-V data. 


Several curves are shown in Figure 5 to 
illustrate several different theoretical models. 
The solid curve gives a good fit at large voltages 
using the parameters of Tables 3 and 4. A dashed 
curve shows the effect of including the Keldish- 
Franz effect in the device. Above about 0.5 volts 
the Keldish-Franz effect is seen to have little 
curves. The curve including the 
Keldish-Franz effect combined with a shunting 
resistance of 5 kO does appear to give the best 
rit to the low current data for this particular 
cell. 


The dotted curve shows a calculation using the 
same p^araeters as the solid curve but neglecting 
heavy doping bandgap reduction affects. The cur- 
rent can be seen to be much too small at high 
voltages. Without the heavy doping effects a very 
low surface layer lifetime would be required to 
get agreement between theory and experiment . 

The parameters which give a good fit to the 
.dark current, also give good spectral response 
calculations as seen in Figure 6. Without heavy 
oping effects the spectral response is also seen 
to be much too large at short wavelengths. 

A comparison beti'reen theory and experiment 
for the light I-V characteristics is given in 
Figure 7, In this case it is seen that the 
inclusion of the shunting resistance and the 
Keldish-Franz effect gives the best fit to the 
experimental data. For both of the t\ra-dimensional 
calculations sheet resistance and contact resis- 
tances of 21 n/p and 0.12 0 respectively ware 
used in the calculations. 


ci^rents be-^een the T _-V characteristic and 
the standard dark I-V SRarSSteristic. This indi- 
cates a very large sheet resistance and/or contact 
resistance. The difference between the curves is 
accurately described by a sheet resistance of 
1500 fl/p and a contact resistance of 1.12 n. 

These values are both very large and much larger 
than can be tolerated in high efficiency solar 
cells. No attempt was made on this cell to des- 
wibe the low current region by either the Keldish- 
Franz effect or a shunting resistance since this 
has little effect on the characteristic above about 
0.4- volts. 


The comparison bati'feen theory and experiment 
for the spectral response is show in Figure 9 . 
Good agreement is again seen over the entire wave- 
length range. Included for comparison are calcu- 
lations for an erfc doping profile at various sur- 
face recombination velocity (SRV) values. A good 
fit to the short wavelength values could not be 
obtained for the erfc profile even with very small 
values of SRV. 

A reasonably good overall fit to the light I-V 
characteristic is shown in Figure 10 . The agree- 
ment in this case is not quite as good as for the 
0.1 n-cm devices. The largest error betoeen theory 
and experiment occurs in V which has about a 2% 
error. This is probably due to a theoretical dark 
current which is slightly too large. This in turn 
could be due to a diffusion length used in the base 
or surface layers which was slightly too small. In 
general, however, the agreement between theory and 
experiment is good considering the large values of 
sheet resistance and contact resistance for these 
cells . 


CELL WITH 16 JJ*CM BASE LAYERS 

The final general class of cells studied were 
BSF cells on 16 n*cm base layers. Relevant device 
data as given in Tables 3 and 4. The devices had 
a very shallow junction depth of about 0.2 pm. 
Because of the very shallow junction depth, the 
diffusidn profile was modeled by a single erfc 
profile. ^ The back surface field region was made 
by alloying Al. This p’*' layer was modeled as a 
Gaussian doped region of doping density 

and 0.5 pm in width. A Ta_0_ antirefleoting layer 
of 535A was present on the cells and this was 
included in the computer calculations. 


CELLS WITH 10 n- CM BASE LAYERS 

The second general class of cell studied were 
fabricated^on 10 Jl* cm material of 10.5 mils thick- 
ness. Additional data on the cells is contained 
in Tables 3 and 4. The junction was diffused at 
SSO C for 30 mins, resulting in a junction depth 
of about 0.57 pm. Parameters used in the calcula- 
tions to describe the doping profile are given in 
Table 4. 


Basically the same models and device param- 
eters were used for these cells as for the 0.1 
SJ" cm cells. A detailed comparison beti'jean theory 
and experiment is shorn in Figures 8 , 9 and 10, 

In Figure 8 a large difference is seen at large 


The terminal I-V characteristics are shown in 
Figure 11 . Two theoretical curves are shown. The 
o^lo'^a-'ted using the experimental 
diffusion length value of 160 pm. This is seen to 
give a current much, larger than the experimental 
values. The dark current could not be made to 
agree with experiment until the base layer diffu- 
sion length was increased to about 460 pm. The 
resulting curve is the solid one in Figure 11 . 

For the non BSF cells studied the measured diffu- 
sion length was found to give a good theoretical 
calculation. However, the failure of the X-ray 
measurement technique to give a reliable diffusion 
length value is not too surprising in the present 
case. The cell thickness of 6,5 mils corresponds 
to 165 pm which almost exactly equals the measured 


diffusion length. The X-ray technique is known to 
become unreliable when the diffusion length becomes : 
l^ger than the ceU thickness [21]. The 460 pm dif- 
fusion length value is also consistent with the 
expectation that 16 n*cm material should have a 
larger diffusion length than lo n-cm material. This 
inconsistency between the measured diffusion length 

performance was observed for all 
thin BSF cells studied in this group. 

dal-;, of spectral response and light Hv 

data are shown in Figures 12 and 13 . The agreement 
IS again very good. For this cell a sheet resis- 
t^ce value of 380 f5/Q and a negligible contact 

® 8°°'^ to both 

the efficiency and the dark current I-V data. 

SUMMARY AND CONCLUSIONS 

involved a detailed comparison 
of the experimental performance of three types of 
^ designs with tk theoretical results 
calculated with a detailed numerical analysis of 
solar cellperformance. The analysis included a 
two-dimensional analysis of sheet resistance and 
contact resistance effects. 

general it has been found that an accurate 
modeling of all the physical effects present and 
of the device doping profile leads to theoretical 
in very good agreement with 
xperimental results. The most Important physical 
effects which were found to be necessary in order 
to obtain good agreement betoeen theory and 
experiment are: 

1. The base layer and surface layer lifetimes 

or diffusion lengths. •‘■etimes 

2 . The doping profile within the diffused 
•surface layer. 

3. The presence of heavy doping bandgap 

reduction effects . ® ^ 

4. _ An accurate modeling of the absorption 
coefficient and surface reflection. 

5. An accurate modeling of the two-dimen- 
sional nature of the surface sheet resistance and 
metal contact resistance, 

correlation of theory and experiment 
■obtained here give confidence that the theoretical 
calculations can be extended to other device 
obtaii^^^ device designs and accurate results 
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Table 1. Comparison of reported absorption eoeffi- 
cient values at long wavelengths [18]. 


Absorption Coefficient (cm“l) 
1.1 (am 1.0 pm 0,95 um 


[13] Dash S Newman " 

7 

100 

220 

[14] Vedam 

- 

_ 

270 

[15] Runyan 

- 

67 

170 

[16] Vol'fson S Subashiev 

- 

64 

150 

[ 17 ] Macfarlon 

3.9 

74 

204 

This work 

3.9 

74 

204 



Table 2. Parameters for the 

two grid patterns. 


1 cm by 

1 cm by 

Parameter 

2 cm 

10 finger 
cell 

2 cm 

9 finger 
cell 

Length of fingers (cm) 

0.9736 

1.949 

Width of fingers (cm) 

0,02 

0.0055 

Thickness of metal (cm) 

0.00019 

0,00019 

Spacing beti^reen fingers (cm) 

0.18 

0.1056 

Bus bar width (cm) 

0.0264 

0.051 

Resistivi-ty of metal (n«cm) 

l.S6xl0~b 

1.56x10-6 


Table 5. Re^stance parameters for two-dimensional^^^ 


Parameter 

Cell 

number 


D-1 d-2 

D-3 D-4 

Sheet resistance 

n/a 

9 21 

1500 380 

Contact resistance 

lO"*^ 0.12 

1.12 lo"** 

' I 


Table 3. Device parameters. 


Parameter 



D-1 

Base layer resistivity 
Thickness (mils) 

0.1 0* cm 
6.0 

BSF Cell 

NO 

Surface profile 

empirical 

^surface 

Sase^^"^ 

Junction depth (pm) 

L*f , "VL 
min med 

60 

1.0 

SRV (cm/sec) 

10^ 

Anti-reflection layer 

absent 


•'•Value actually~used in calculations. 
Table 4. Parameters of surface diffused layer. 


Parameter 


Cell Mumber 



D-1 

D-2 

D-3 

D-4 

(pm) 

1.0 

0.72 

0.57 

0.2 

C^(cm“^) 

4.xl0^° 

4x10^° 

2x10^° 

2x10^° 

Cg (cm ^) 

8x10^® 

8x10^® 

2x10^^ 

2.xlQ^° 

Cpm) 

0.4 

0.20 

0.07 

0.0 


Cell number 


D-2 

D-3 

D-4 

0.1 n*cm 

10 f2* cm 

16 n»cm 

10.2 

10.5 

6.5 

NO 

NO 

YES 

empirical 

empirical 

90 

2x10 erfc 

L . 
min 

105 

■5:(L • +D ) 
2 min med 

230 

• +1, 

2 min mec 
160(460*) 

0.72 

0.57 

0.2 

10^ 

2x10*^ 

5x10® 

absent 

absent 

59 5A Ta^Og 







by triangles.^ material except those marked 



Figure 2. Dark I-V characteristics of Cell D-1, 





Figure 3. Spectral response of Cell D-1. 



Figure 6. Spectral response of Cell D-2. 





nh theory and experiment for 

Srv2 f I-'" characteristics of Cell 0^2. 
Curve<? n ^ parameters of Table 1. 

Fv. F include a shunting current and 

Franz-Keldysh effect. ana 
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Figure 8. Dark I-V characteristics of Cell D-3. 



Figure 11 , Dark I-V characteristics of Cell D-4. 





Figure 10. Comparison of theory and experiment for 
the photovoltaic I-V characteristics of Cell D-3. 



Figure 13. Comparison of theory and experiment for 
the photovoltaic I-V characteristics of Cell D-4, 
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Appendix 10.2 a two dime nsional analysis of sheet pxsistajice ajjd 

CONTACT RESISTANCE EFFECTS IN SOUR CELLS^ - 
C. R. Fang and J. R. Hauser 


North Carolina State University 
Raleigh., NC 27650 


ABSTRACT 

Most studies "of contact resistance and sheet 
resistance effects in solar cells have modeled 
these effects m terms of ’a lunroed resistance in 
series with an ideal solar cell' The two-dimen- 
sional nature of the distributed sheet resistance 
from the smface layer, however, makes the accu- 
rate modeling of a solar cell more involved than 
asimple series resistance. If an ideal one-dimen- 
sional solar cell has a dark current which varies 
as e:fp(qV/kT), the terminal current of a cell 
including sheet resistance will have a voltage 
dependence which approaches exp(qV/2kT) at large 
currents due to the sheet resistance ohmic volt- 
age drops. This effect cannot be accurately 
modeled by a lumped resistance. In this work an 
analytical model and detailed numerical calcula- 
tions have been studied for describing the sheet 
resistance and contact resistance effects in 
solar cells . 


INTRODUCTION 

The important effects of contact resistance 
and the sheet resistance of the surface layer have 
been recognized for some time. Series resistance 
effects ^e especially important in determining 
peak efficiency and cwve factor of a cell and are 
of lesser importance in determining short circuit 
current, open circuit voltage and spectral 
response. Analytical studies have been conducted 
by several workers [1,2] developing exoression for 
the series resistance to be used in modeling solar 
cells. ^ As discussed later these annroaches all 
have limited application because the two-dimen- 
sional nature of the sheet resistance cannot be 
accurately modeled by a lumped resistance unless 
one is considering only small changes in voltage 
about some operating point. 

“ igure 1 shows a typical solar cell structure 
with a conventional contact finger geometry. The 
iront surface metal consists of a main bus bar to- 
which metal fingers which collect the current are 
attached. Such a structure has the following 
major sources of resistance: 

1. Bus bar metal resistance. 


’‘This section .'has heen published in 
the Proceedings of the Thirteenth 
Photovoltaic Specialists Conference. 


2. Finger metal resistance, 

3. Front surface contact resistance between 
metal grid and semiconductor, 

4. Sheet resistance of the semiconductor 
layer at the siu'face, 

5. Base layer bulk resistance, 

6. Back surface contact resistance between 
metal and semiconductor. 

In p^ticular solar cells any or all of these could 
provide the major source of series resistance, but 
for solar cells operating at 1 sun, the two major 
sources of resistance are usually contact resis- 
tance and surface layer sheet resistance. This 
assumes that the contact metal grid structure has 
been designed with sufficient thickness of metal 
that voltage ^ops along the contact fingers are 
smal_. This is readily done at 1 sun intensity. 

For multi-sun operation it becomes more and 
more difficult as the intensity is increased to 
ac^eve a good contact grid design. The spacing 
between the grids must be decreased as well as 
^eoreasing the width of the contact stripes. Ohmic 
^ps along the contact fingers tend to become more 

^portant unaer multi-sun conditions. The base 

ayer resistance can also become more important. 
However for BSF (back surface field) ceUs in which 
the base layer is less than a diffusion length, 
conductivity^ modulation effects in the base tend 
to minimize base layer resistivity effects, even 
under multi-sun conditions . 


:he ma^or emphasis of this work has been on 
the contact resistance and sheet resistance as they 
aominated the solar cells studied. The two-dimen- 
sional computer analysis, however, was designed to 
include all of the effects discussed above except 
for any bus bar resistance, i_.e., the bus bar was 
assumed to be at a constant potential. 

SOLAR CELL MEASUREMENTS 

There are three fundamental sets of terminal 
l-J measurements which are influenced in different 
ways ny the series resistance and two-dimensional 
e.xects. In^the first method the terminal I-V 
ch^acteristic.is measured under a steady illumin- 
ation. These measurements may be made for AHC , 
spectrum and ’under 1 sun or 
m_ xi -sun conditions . In any case by varying the 


load resistance the power quadrant of active 
solar cell operation is obtained. The important 
Pieters detemined from this are sLS^ScSt 
f^or voltage, efficiency, fill 

course the roost iropor- 
, solar cell performance. For 

e first order model of an ideal solar cell in 
series With a ceU resistance, this curve is 
described By the equation 

“!)!’^,?sc! ^o.“^ ^ constants for the cell 

and illuminafion conditions. 


for example Equation (2), if one assumes that 
Rg IS a function of current level. A series 
varies with current level has 
indeed been found and reported by several workers 
— However, the almost exact factor of 2 in 
the slope of the two curves leads to the suspicion 
that there is some fundamental physical process 
responsible for this and that the effect is not 
simply a current dependent resistance. This is 
indeed the fact as shown in the next section where 
the two-dimensional nature of the sheet resistance 
is shown to account exactly for this factor of 2 
in -slope. 
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In the second technique the dark forward I-V 
cnaracteristics are measured by applying a ter- 
minal voltage and observing the current. If a 
sol^ cen can be described as an ideal diode in 
senes with some resistance, the current in this 
measurement can be expressed as 


where 


X=I^{expCq(V-IR^ )/kT]-l} , 
^o \ constants. 


( 2 ) 


The third method referred to here as the 
fsc oc a varying illumination level 

and measures open circuit voltage and short 
circuit current. A plot of the corresponding 
points then gives an I-V curve, which according 
to -ae model of an ideal cell in series with a 
senes resistance can be expressed as 




(3) 

independetly by Heeger 

[3J , Wolf Cl] and Queisser [4] . The advantage of 
this technique is that the effects of series resis 
uance_OTe eliminated, according to first order 
tneones. 

Figure 2 shows experimental data measured on 
one particular solar cell by the two technicues 
discussed above. Attention should be directed to 
the regions of the curves for which voltage >0.S 
volts. Below about O.S volts the curves are” 
dominated by depletion region currents and shunt- 
ing resistance currents and the ideal diode theorv 

^sc-'^oc 

? expected (See Ecuation (3)) 

exp(qV/kT) dependence as indicated by the solid 
curve. The standard dark current curve, however, 

IS seen to follow almost exactly an exp(qV/2kT) 
^ependence. This cannot be described by Eauation 
I ; -or any constant value of series resistance. 

The type of behavior seen in Figure 2 has 
^een observed on a wide range of solar cells. The 
:sc",oc. currents follows almost the 

la.^ diode behavior. Deviation from the ideal 
exp(qV/kT) behavior are always observed at low 
currents and are again observed at high currents 
wnere .-jigh injection effects begin to occur. In 
t..e_ region where the standard dark current data 
begins to deviate from the I-V data the slope 
o. .he curve on a semilog grlgh ?inds to be almLt 
exactly pne-half that of the I -V data. This 
type of behavior can only be explained bv the 
model of an ideal diode in series with a” resistor. 


DISTRIBUTED SHEET RESISTANCE 

The sheet resistance of the sin?face laye;? of 
a solar cell is one of the most important, if 'not 
the most important, contribution to the resis- 
tance effect in solar cells. Figure 3 shows a 
cross sectional view of a solar cell illustrating 
the distributed nature of the surface layer sheet 
resistance. Consider the dark I-V characteristics 
of the cell generated by applying a voltage to the 
teimiinals as in Figure 3. The current which 
crosses the p-n junction must flow in a direction 
transverse to the p-n junction to reach parts of 
the junction away from the contact. This causes 
a ^tential drop along the surface layer V(x). 

This in turn causes regions of the junction away 
from the contact be be less heavily forward 
biased. For large voltages and high terminal 
currents,, the current density is crowded into an 
area near the contacts and uniform current density 
no longer occurs. 

The current crowding problem described above 
for sol^ cells is almost identical to that which 
occurs in bipolar transistors due to base current 
flow transvers to the p-n junction. In fact for 
analysis purposes the solar cell looks like a 
transistor with zero cuirent gain. The current 
crowding problem has been extensively studied for 
transistors and the results can be readily applied 
,to solar cells [6] . 


The current crowding problem can be thought 
of as giving rise to some effective width, L _j. in 
Figure 3, over which one can consider uniform 
injection at the junction to occur. The current 
on each side of a finger can then be expressed as 

Ip=JoL„L^^^exp(qV/kT), (4) 

where ^ is the length of the finger and L „ is 
the effective distance over which current flow 
occurs. Of course at low voltages L /2, 
where Sp is the spacing between the lingers. 
Equation (4) neglects the current due to the area 
under the fingers and bus bar. This is usually 
small due to the small metal coverage factor. 

As tFifi current density increases, the effec- ‘ 
tive width decreases . When severe current crowd- 
ing occurs, it has been shown that C6] 


^eff 


2kTLp 



(5) 
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Khere pg„ is the sheet resistance of the surface 
layer. This shons that the effective width 
decreases ^versely with total current. Combining 
Equations (4) and CS) and solving for the total 
terminal current gives 


I =/aj 


exp(qV/2kT), 


where 




- 2kT h** 

~ q p 


ST 


( 6 ) 


(7) 


with h equal to the perimeter of the contact grid 
pattern. Since current crowding can occur around 
the entire metal surface, it is clear that h 
should be the total perimeter of the contact grid 
which includes both the contact fingers and the 
bus bar. 

Equation (6) has the correct voltage depen- 
dence to describe the high current experimental 
results as shown in Figure 2. The difference 
between the 1 -V data and the conventional dark 
Client data is seen in this model to be due to 
the sheet resistance, but the effect cannot be 
simply expressed as a constant resistance in 
series with an ideal solar cell. The current I 
has a simple physical interpretation in terms or 
the measured 1-V characteristics. If one solves 
the ideal current equation 


I = AJ exp(qV/kT), 


( 8 ) 


and Equation (6) simultaneously, one finds that 
Ij is the intersection current of the two curves. 
This provides a simple means of determining sheet 
resistance from experimental I-V data. Two curves 
are drawn as shown in Figure 4. One curve fitting 
‘sc"'^OC ® second curve with 1/2 the 

slope of fhe I -V data to describe the conven- 
tional dark I-v data. The intersection of the two 
curves then determines and from this the sheet 


resistance is calciaated as 
_ 2kT 

q 


ST 


I,A> 


(9) 

where h is again the total perimeter of the metal 
area Md A is the total solar cell area. This 
technique has been used to estimate the sheet 
resistance of a number of solar cells with very 
pod results. The values obtained in this manner 
have pen used in a detailed two-dimensional 
numeriop calculation with good results as dis- 
cussed in detail in the next section. 

A large metal sheet resistance can lead to 
a second type of current crowding along the length" 
of the metal fingers instead of perpendicular to 
“t e metaX fingei’s as occurs whan t;fie saiijXconduclro 3 ? 
surface layer is the limiting factor. This could 
occiir for an improperly designed contact structure 
or perpps in cells operating under multi-sun 
cpdpions. This distributed cutoff effect along • 
the ringers will also lead to a dark current which 
varies as exp(qV/2kT), The presence of this 
pfect can be predicted quite readily by calculat- 
ing the expected voltage drop along the contact 
fingers at any given current level. If the 
expected voltage drop assuming a uniform current 


collection along the fingers exceeds kT/q volts 

** distributed resistice 
effect due to finger resistance. 

In ppnciple, a solar cell could have distri 
iuted resistance effects due to both the semiS^ 

i V expected to have a slope of 1/4 

This would lead to a 
ppent depenaence of exp(qV/4kT). This tvne of 

_ After accounting for the distributed sheet ' 
experimental data of Figure 
IfJ- ! ^ ‘^^^^'^ation at large currents from the 
predicted exp(qV/2kT) voltage dependence/ tSs 
pcess voltage is identified as R--l in 4 

For an the cens studied in thi^wLk! 

rat^5T be a“u- 

wt^Ww^ ^ constant resistance in series 

wi^ the teimiinal current. There are several 
potentip. sources for this resistance such as bus 

resistance, and fj^nt 

pd back surface contact resistance. For the 

r studied ,in this work, the metal-semi- 
Mndptor contact resistance has been identified 

DetSLr^^ s^'^rce of this resistance. 

Ssed InTlaf"®''" resistance are dis- 

cv^sea in a later section. 

An anpytical model describing both the low 

h.T'" “51. =»r.„t ai,tra„t.d rSiSlnce ' 

JS ^ oeU cm b. constructed from 

the two limiting cases of Equations (6) and (8) 
Solution of the equation v y 

^ ^^T = ^o^T®=‘P'^'l^''“\'T)/kT], (10) 

c!l^rent°at ™ approximation to the 

tSr?Iofifa r f®' ^"'^Tuded in Equa- 

resistance R' to 

?Jon-distrihuted resistance s§ch as 
contacu resistance. 

should be pointed out that although a 
s^le analppal model can be used to describe 
the tempp dark I-V characteristic, no such sim- 
Lted^-L ^^^^bablG for describing the distri- 

Snf Thl illuminated 

cells. The model of a diode in parallel with a 

employed for the complete 
naii resistance effects are domi- 

n^t. The only known technique for accurate calcu- 
lations undeJ. these conditions is a two-dimJnsSnal 
numerical technique. 'dimensional 

COMPUTER CALCULATIONS 

A ^«siled two-dimensional computer model has 
L""” distributed- resistance 

e— ects _n solar cells. The model which has been 
ueec consists of a f.o-dimen.sional array oTidlll 
e dimensional solar cells interconnected by 
series resistance. An array of NX and NY points 
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between two grid fingers are used to locate an ele- 
ment of the two-dimensional array. The number of 
grid points in the X-direction is NX and the corre- 
sponding number' in the Y-direction is NY. Figure 5 
shows the equivalent circuit of the incremental two- 
diirenslonaX modcX wh^re the current source end 
diode is that of a one-dimensional model including 
both optically generated current and forward bias 
current. For the present work this current was 
calculated fj^om the tabulated current vs. voltage 
, values obtained from a one-dimensional computer 
analysis. Rg_ is the sheet resistance of the semi- 
conductor surface layer and R is the contact resis- 
tance between the metal-semiconductor interface. 

The collecting metal is also assumed to have a 

resistance Rjj which is included between arrav 
elements located on the edges of the grid fingers . 


The distributive resistance elements can be 


a. Sheet resistance: 

^H-l'^tr'^ST AY* 

(10) 

b. Sheet resistance: 

^1-NX“\+NX"^STAX’ 

(11) 

c. Contact resistance: 

= Ri = V^c* 

(12) 

d. Metal resistance: 


(13) 


Liic apdcings oetween gria points 

in the X- and Y-direction respectively, R-„ is the 
sheet resistMce in fJ/square, R is the contact 
resistivity in fl*cm'^, and R is the metal resistiv- 
ity in cm, Aj, is the incremental area of the 
metal contact and T is the thickness of the metal 
contact grid. 


Any bulk resistance R_, in Figure 5 which 
arises from the bulk resistivity of the base layer 
need not be included here, sinde it is already in- 
cluded in calculating the one-dimensional I-V char- 
acteristics. Also the distributive base resistance 
Rj 2 Figure 5 can be usually neglected, since 
most solar cells are covered with an ohmic contact 
over the entire back surface . 

The bus bar is assumed to be at a constant 
potential which equals the terminal solar cell vol- 
tage. The voltage of all other grid points can 
then be calculated from the simultaneous solution 
of the node voltage equations at each grid point. 
Details of the solution method are discussed else- 
where C73. The input data to the two-dimensional 
computer program consists of the one-dimensional 
1-V data plus structural data on the contact fin- 
ger arrangement. An iterative solution technique 
based on a modified Newton-Raphson technique has 
been used to solve the nonlinear coupled node equa- 
tions. Calc^ations have typically Been made until 
the voltage is accurate at each array point to less 
than 10~*t volts. 

The two-dimensional program provides calcula- 
ted values of voltage and current density at each 
node point between the grid fingers as well as pro- 
vi^ng terminal I-V calculations including the dis- 
tributed resistance effects of the solar cell. 


COMPARISON OF THEORY AND EXPERIMENT 

The results of the two-dimensional calculations 
have been compared with experimental data for a num- 
ber of different solar cells. Results are reported 
here for three specific solar cells with widely 
varying sheet resistance values. 

The first requirement for modeling the distri- 
buted resistance effects is to have an accurate 
model or representation for the intrinsic properties 
of a one-dimensioncil solar cell neglecting these 
effects. In this work this was obtained from a 
computer modeling of the one-dimensional solar cell. 
Details of this work are reported elsewhere [0], 

For the pilose of this work it can be stated that 
the one-dimensional calculations were found to agree 
very well with the experimental I -V data such 
as shown in Figure 2. Thus for t^e purposes of this 
present work, the ^g-“V data can be taken as the 
one-dimensipnal solAr cell 1— V characteristics used 
in the two-dimensional computer calculations. 

Figure 6 shows measured I-V data for a good 
solar cell in which there is little difference in 
the measured “V data and the conventional dark 
current I-V data. This cell had a junction depth 
of 0.7 pm and a base layer of 0.1 £5* cm resistivity. 
The sheet resistance and contact resistance values 
which give the best agreement between theory and 
experiment are 21 fl/a and 0.12 H, respectively. 

The calculated points in Figure 6 show very good 
agreement between theory and experiment. 

Figure 7 shows the corresponding agreement 
between theory and experiment for the same cell 
under AMO, 1 sun illumination. The agreement is 
again very good. Because of the small resistance 
values, good agreement can also be obtained for 
this cell using a lumped series resistance model. 

Data for a solar cell with a larger value of sheet 
resistance has been previously shown in Figure 2. 

This' particular cell had a junction depth, of about 
0,3 ym, resulting in a sheet resistance of about 
940. n/t 3 . When this was used along with a contact 
resistance of 0.2 n good agreement between theory 
and experiment was obtained from the two-dimensicnal 
calculations. Figure 0 shows the good agreement 
Between theory and experiment in this case. 

Finally data for a solar cell with a very 
large sheet resistance and contact resistance are 
shown in Figure 9. This data can only be explained 
by severe current crowding. This cell had a junc- 
tion depth of 0.07 ym and the value of sheet resis- 
tance which best describes the experimental data is 
8200 O/b. The sheet resistance is somewhat larger 
than would be expected from a junction depth of 
0.07 ym. However, with such a shallow junction, 
the junction depletion region may occupy a large 
part of the layer accounting for the large sheet 
resistance. The contact resistance at the front 
and/or back surface was also found to be large at 
a value of about 2.84 fj. Some of this large value 
may arise from the back surface contact to the 
10 n*cm base layer of this particular cell. 

Although not shown, even with such large 
values of sheet resistance and contact resistance. 
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the agreement between theory and experiment for the 
illuminated characteristics are very good. For 
such large resistances it is not possible to obtain 
good agreement ‘between theory and experiment with- 
out the use of the detailed two-dimensional compu- 
ter analysis of the distributed resistance effects. 

SUMMARY AND CONCLUSIONS 

In this work the two-dimensional nature of the 
sheet resistance and contact resistance effects in 
solar cells has been discussed. One of the major 
results of this work is the demonstration* both 
■Dheoretically and experimentally, that the distri- 
buted nattire of the semiconductor sheet resistance 
causes the terminal dark I-V characteristics to 
exhibit an exp(qV/2kT) type dependence even when 
the one-dimensional characteristics of the cell 
exhibit an exp(qV/kT) type voltage dependence. The 
analytical model developed also provides an easy 
method for estimating the sheet resistance of a 
solar cell from the terminal I-V data. 

When the distributed nature of the sheet resis- 
tance is important in the terminal I-V character- 
istics it is not possible to accurately model a 
solar cell by an ideal diode and current source in 
series with a fixed resistance. A two-dimensional 
computer analysis using an array of grid points can 
be used to accurately model the solar cell. The 
results of such an analysis are in agreement with 
the . analytical model developed here for the dark 
I-V characteristics of a cell. The computer analy- 
sis has been found to be in very good agreement 
with experimental results of sheet resistance and 
contact resistance effects for a wide variety of 
solar cells. By comparing theory and experimental 
results it has been possible to determine semicon- 
ductor sheet resistance and metal-semiconductor 
contact resistance for a number of silicon solar 
cells. 
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Figure 1, Solar* cell with conventional contact 
finger arrangement. 



Figure 2. Terminal I-V characteristics of cell 
D-1. 
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Figure 3. Distributed nature of the surface layer 
sheet resistance, 
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Figure 4. Experimental method for determining 
sheet resistance. / 
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Figure 6. Terminal 1-V characteristics of cell D-2. 



Figure 8.^ Comparison of the photovoltaic charac- 
teristics of cell D-1. Various curves are: (a) 

one— dimensional with total area^ (b) one— dimen- 
sional with active area; (c) two-dimensional 
with no series resistance and; (d) two-dimen- 
sional with series resistance. 



Figure 5. Incremental model of two-dimensional 
solar cell. 



Figure 7. Comparison of the photovoltaic I-V char- 
act^isties of cell D— 2. One-dimensional calcu- 
lauions- are witE total area Cdashed curve). Two- 
diisensitonal calculations include series resis- 
tance . 




Figure S. Terminal I-V characteristics of cell D-3 
With very large sheet resistance. 
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Appendix 10.3 

. Impurity Gradients and High Efficiency Solar Cells* 

C. R. Fang and J. R. Hauser 
North Carolina State University 
Raleigh, NC 27607 

ABSTRACT 

One potential means of improving the efficiency of solar cells 
especially after space irradiation is to incorporate built-in fields 
into the device through the use of impurity doping gradients. 

Previously published papers have indicated an improved minority 
carrier collection efficiency and improved efficiency when doping 
gradients are present. In this work a detailed numerical calculation 
of solar cell performance has been used to study various types of 
doping gradients. In general the predicted improvements in performance 
have been less than previously reported due to various device effects 
such as high injection and the dependence of .lifetime on doping density. 


'This section has been accepted for publication in Solid-State Electronics. 



I. INTRODUCTION 
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There have been a number of theoretical investigations into the 
efficiency of drift field solar cells in recent years [1-4.1. In these 
earlier investigations the emphasis has h?en on the calculation of 
collection efficiency as a result of including drift fields, with some- 
what different approximations in lifetime, mobility and other important 
parameters. Wolf [1] was the first to demonstrate that a surface drift 
field helps the short wavelength response and to show that a base field 
can considerably reduce the effects of radiation damage and increase 
the useful life of solar cells. In 1967 Bullis and Runyan [3] found 
that there exists an optimum field width which is approximately twice 
the diffusion length or about 25 y whichever is shorter in n'^-p cells 
for maximum collection efficiency. Van Overstraeten, et al. [4] have 
shown that the advantage of the drift field is mainly determined by the 
layer close to the collection junction. To obtain higher collection 
efficiency, the magnitude of the drift field should be high, at its optimum 
width and the doping density should be as low as possible. 

In considering impurity doping gradients, it is useful to distinguish 
two somewhat different types of doping gradients. First, if the doping change 
IS very rapid, space charge regions form. On the other hand if the change 
in doping density is very gradual, a state of quasi-neutrality exists at 
every point and a built-in electric field is produced throughout the bulk 
of a solar cell. A gradient sufficient to give rise to a space charge 
region will be referred as a high-low junction while a gradual 

change in doping will be referr-^Vto as a drift field region. 
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Despite the theoretical predictions of the drift field enhancing 
the collection efficiency, the advantages of drift fields are still 
unconfirmed experimentally, especially with regard to total conversion 
efficiency [5]. The purpose of this paper is to explain the physical 
mechanism governing the operation of drift field celJsand the discrepancies 
of those earlier results. The emphasis is on maximising the total conver- 
sion efficiency as a result of including the base field. In this paper 
. the major mechanisms which tend to limit the conversion efficiency due to 
the incorporation of drift field in Si solar cells are discussed. 

The use of a high-low junction in solar cells is a more recent 
innovation than the drift field concept. However, the advantages of such 
abrupt doping gradients have been experimentally demonstrated and incor- 
porated into the design of high efficiency solar cells. Basically the 
high-low junction near the back surface prevents minority carriers from 
reaching the back ohmic contact and thus acts as a minority carrier 
reflecting boundary [6]. 

II. FIRST ORDER THEORY 

The structure, of the basic solar cell to be considered here is 
shown in Figure 1. The n^^ surface layer is a thin (0.1-0.5 ,m typically) 
heavily doped layer which may have a doping gradient. The p-type base 
layer is shown with a high-low junction near the back ohmic contact and a 
wide, drift field region. Because of the major role of the base layer 
properties in determining solar cell performance, the discussion here 
concentrates on the base layer. 

Within the drift field region of the base, excess minority carrier 
electrons (n) are controlled by the basic drift-diffusion and continuity 
equations which when combined into a single equation for steady state gives 
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1 T 

~ ]n + G = 0, 
n 


(J ) 


where G Is the optical pair generation rate. For a general doping profile 
Equation (1) must be solved by numerical techniques. Also the dependence 
of mobility and lifetime on doping density complloates the solution. 

The special case of a constant electric field (exponential doping 
profile) with constant and has been frequently used for drift cells 
and provides certain insight into drift cell operation. If the simpler 
equation for these conditions 

2 - _ _ 

n d n , _ dn n 

^dx^ 

is considered, an exact solution can be obtained of 



n = exp(x/L^) + exp(-x/L") + t^G, 

(3) 

where 

1/L“ = /(1/L^)^ + (qE/2kT)^ (qE/2kT). 

(4) 


The basic theory behind the drift field solar cell is that the electric 
field increases the collection depth for optically generated minority 
carriers within the base layer since l'*'>L 

n ' 

The value of a constant drift field can be written as 


E 


kT 1 _ kT 

q dx ~ qW 



(b) 


where W is the width of the drift field region and N^, are the doping 
densities at the ends of the drift field region. The ratio N^/N^ is 
limited in value to around 10^ in practical cases and this sets an upper 
limit to the field which can exist. From Equation (5) it is seen that the 
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largest field occurs when W is as small as possible. However, if W becomes 
less than L , Equation (4) begins to lose meaning since the field doesn't 
exist over the assumed dimensions. Thus by setting W larger than or equal to 
and combining Equations (4) and (5) a limit on can be established as 


i^lTjlnCN^/N^). 


( 6 ) 


5 + 

With a ratio of 10 for N /N this gives L <3.54 L . Thus an enhancement in 

^ X n 

collection depth by a factor of about 3 for the conditions stated above 
is the most that can be realized due to doping density limitations. 

Previous studies of drift field cells have concentrated on the 
improved collection efficiency . Equally important in determining solar 
cell efficiency is dark forward current and open circuit voltage which, 
according to first order device models, are related as 


^dark = ^o ®^P(qV/kT) 


(7) 


oc 


^Jln(J /J ), 
a sc o’ 


( 8 ) 


where J is short circuit current density and J is a function of device 

sc O 

geometry, and material parameters. For a drift field cell with constant 
mobility 

2 


° L /L«n(N„/N, ) "l 


n. 

1 

N. 


( 9 ) 


n 


2 ' !• 


For a large open circuit voltage, should be minimized. A large drift 
field requires that be large and be small, while Equation (9) shows 
that this is far from the condition for minimum dark current where should 
be as large as possible. Thus large drift fields and minimum dark current 
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are incompatible objectives. For many potential drift field designs, 
it appears that increases in dark current almost completely eliminate 
any increased collection efficiency due to the drift field. 

The previous considerations have been for cells with a base layer 


thickness large compared with a diffusion length. These considerations 
must be modified considerably if the cell thickness is comparable or 
less than the dlffusl,|n length. For such thin cells, the diffusion 


length is not the major' factor determining collection efficiency and a 

back surface high-low junctjjon has been found to be effective in increasing 
the efficiency of such cells. 1 ) 


The purpose of a high-low junction such as shown in 


r,. 

Figure i is to 


prevent minority carriers created either by light or forward injection 
from recombining at the back surface. The high-low junction interface 
can be characterized as a low surface recombination boundary to the base 
layer. The high-low junction is most effective for solar cells in which 
the base width is less than a diffusion length. As discussed elsewhere 
[6-8], the high-low junction can be very effective in reducing back surface 
losses, giving a collection depth approximately equal to the total cell 
thickness. An ideal high-low junction minimizes dark current and maximizes 
short circuit current and open-circuit voltage simultaneously. 


III. NUMERICAL CALCULATIONS 

In order to investigate in detail the ideas discussed above, a series 
of calculations have been made of solar cell performance with various types 
of drift fields and high-low junctions. The analysis consists of a detailed 
numerical solution of the semiconductor device equations. Details of the 


p(x) DRIFT FIELD REGION 
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modeling and solution techniques have been discussed elsewhere [9], 

The solutions include all types of nonlinear device effects such as 
high level infection as well as doping and field dependences of device 
parameters . 

In any comparison of solar cell performance with and without drift 
fields, the type of cell and the parameters which are held constant in 
the analysis are very important. In this work the approach has been 
as follows. Since a uniformly doped- base layer solar cell with a back 
surface high-low junction has so far demonstrated the highest experimental 
, the approach has been to take this cell as a reference cell 
and to see if other types of doping gradients improve the efficiency. 

The analysis has been performed for a cell with the device parameters 
listed in Table I. The density of 2x10 /cm^ is chosen as a constant 
surface density, since this is approximately the highest doping density 
obtainable after considering heavy doping effect which doesn't produce 
any retrograde surface field [9,10], 

SiO is used as the antireflecticn layer in its optimum thickness of 
800 A where it allows 45,4 mA/cm current density to be available for 
collection. Despite the short wavelength absorption in SiO, it makes 
little difference in the collection efficiency of the base layer since 
most of the short wavelength photons' are absorbed in the surface layers. 

To avoid problems in interpreting the results, heavy doping effects were 
not considered throughout this work. However, such effects as bandgap 
reduction would not significantly change the results because of the low doping 
density used for the surface layer in the calculations. The irradiance condition 
used is AMO and total optical reflection at the back ohmic contact was assumed. 


TABLE I. DEVICE RTPUCTURE PARAMETERS 
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Overall cell thickness 

150 ym 

n*" surface thickness 

0.2 ym 

p base thickness 

144.8 ym 

thickness i; 

5 ym 

n+ surface doping ;l 

2 X 10 /cm 

p base doping 

Optimized 

p+ doping 

10^^/cm® 

Lifetime model 

lies Cio] 

Surface recombination velocity 

10^ cm/sec 

Antireflection layer 

SiO, 800 A 

Irradiance 

AMO 


TABLE II. RESULTS OF SOLAR CELL CALCULATIONS FOR DIFFERENT 
DIFFUSION LENGTHS AND DIFFERENT DRIFT FIELDS 




N^(cm 

N2(cm“3) 

J 

sc 

(mA/cm2) 

V 

oc 

(Volts ) 

n 

(%) 

CF 

J 

o 

(mA/cm^) 

2 

(A/cm ) 

0 

Max 

5. 6x10^6 

5.6x.1016 

42.40 

0.690 

18.19 

0.841 

4.44x10-13 

2.87x10-10 

10 

Max' 

3.24x10.17 

9.9x1014 

42.79 

0.693 

17.79 

0.812 

1.45x10-12 

3.73x10-1° 

20 

Max 

6 . 49x1017 

6. 04x10^2 

42.82 

0.698 

17.60 

0.798 

6.24x10-11 

1.11x10-3 

0 

Med 

9x1016 

9x1013 

35.82 

0.623 

13.70 

0.830 

7.43x10-13 

2.96x10-3 

10 

Med 

3 . 6x1017 

l.lxlolS 

40.83 

0.561 

13 .71 

0.810 

7.19x10-11 

1.92x10-3 

20 

Med 

7.2x1017 

6.7x1012 

42.09 

0.524 

11.79 

0.724 

4.29x10-3 

7.7x10-8 

0 

Min 

5x1017 

SxlO^"^ 

22.63 • 

0.613 

8.39 

0.818 

4.69x10-12 

2.77x10-3 

10 

Min 

2. 9x10 13 

8.8x1013 

27.-13 

0.501 

7.91 

0.787 

4.12x10-1° 

3.39x10-7 

20 

Min 

5.8x1018 

5. 4x10 13 

31.29 

0.385 

6.35 

0.713 

7.53x10-3 

1.83x10-3 
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Another very important factor in any solar cell study is the assumed 
relationship between diffusion length, or lifetime, and doping density. 
Three different relationships as shown in Figure 2 have been used in this 
work. These curves are based upon experimental data by’ lies, et al . [11] 
with the top and bottom curves representing limits near the top and bottom 
of measured diffusion length data, while the middle curve represents values 
near the center of the measured range. These curves are convenient for the 
150 pm cell thickness since the three curves represent, at light doping 
densities, cases where the diffusion length is much larger, approximately 
equal and much less than the total cell thickness. The behavior of the 

cell and the importance of drift fields depends on which of these three 
cases is considered. 

For a given cell thickness and given diffusion length curve, there 
is an optimum uniform base layer doping density which results in maximum 
efficiency, and this condition was first determined. Table II (E^= 0 case) 
lists the optimum uniform doping density for nmximum efficiency and the 
corresponding device parameters for the three diffusion length cases. Also 
shown in Table II are perfornance parameters for cells with drift fields 

(Ej) of 10 and 20 V/cm. For each different case the average doping densitj 
remains constant as the field increases. 

and Nj are the doping density near the high-low junction and p-n 
junction respectively. It is clear that the collection efficiency as 
measured by I^^ indeed increases after building in the drift field with 
the larger Increase for the lower base layer diffusion length. For the 
mnxmum ease increases from 22.6 mA to 31.3 mA 
increases. 


as the drift field 



Figure 2. Model of the diffusion length as 
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Figure 3 shows the minority carrier current density as a function of 
position in the base region for the calculations of Table II. It can be 
seen in the table that the increase in collection efficiency is very 
limited for the maximum diffusion length used here which is already 
larger than the base width. Figure 3 shows that the collection distance 
can be effectively expanded all across the base layer with the appropriate 
drift field for the case of the medium diffusion length. The curves for 
the minimum diffusion length indicates a much smaller increase in collection 
depth. Since collection occurs only over part of the base, this suggests 
that an inproved cell can be obtained by grading only over part of the base 
layer near the p-n junction. This has been verified and is discussed later. 

With the inclusion of the base field, the doping density near the p-n 
junction (N^ of Table 3) is inevitably reduced assuming that the average 
doping density remains constant. This reduction in doping density has 
profound effects on the cell operations. In analyzing the calculated results 
the dark current was approximated by an equation of the form 

J = Jq exp(qV/kT) + exp(qV/2kT) . ( 20 ) 

The and values which give the best fit to the calculated dark current 

curves are listed in Table II. The forward injection current density as 

evidenced by the term increases as a result of lowering the injecting 

barrier, and this typically results in a reduced open circuit voltage as also 

seen in Table II. The depletion layer expands which leads to a larger space 

charge recombination current density as evidenced by the increased J term. 

R 

Third, and perhaps most importantly is the fact that high injection may 

easily occur at a voltage well below the operating voltage of the maximum 
power point. 



20 


Figure 3 


LENGTH ( uM ) 

Minority carrier denaitier. in the base region f-'or 
different built-in drift ti eld and ba«e dirfur.Lon 


rnodely 
lc:ngth . 
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In Table III the transition voltage V* between the low and high 
injection regions is calculated from the equation 

V = 2 ^Jln(N_/n.). 

n q 2 1 

The table also shows the transition voltage Vg where the depletion region 
current equals the ideal injection current density. It's clear that the 
cell is operating at a voltage well above the high injection limit for 
the case of 20 V/cm for both the maximum and medium diffusion length cases. 

The high injection effect is the major reason for the reduction of the curve 
factor CF in Table II and the overall drop in conversion efficiency. 

Although there is a small reduction in CF due to the increased depletion 
region recombination current density, the major reduction is due to the 
deterioration of the open circuit voltage as a result of a much higher 
forward injection current density. The increased dark current which is 
predicted by the approximate model of Equation (9) due to a decreased 
is definitely verified by the computer calculations . 

From the results of Table II several conclusions can be drawn with 
regard to solar cells with the same average base doping levels. First* a 
drift field enhances the collection efficiency with the largest changes 
occurring when the diffusion length is much less than the base layer thickness. 
Second, almost all of the increased collection efficiency is offset by 
increases in dark current and a reduced open circuit voltage. Third, for 
large drift fields, operation near the p-n junction tends to be in the 
high injection region which further tends to reduce efficiency due to a 
reduced curve factor. 
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Table IV shows calculations-) for cells with the same electric field 
but with different average base layer doping densities. In comparing 
this data with that of Table III, it is seen that the highest efficiency- 
results from a compromise between retaining a large open circuit voltage 
and achieving an enhanced collection efficiency due to the drift field. 

The calculations have shown that for highest efficiency the doping density 
near the p-n junction should be kept near the value which gives highest 
efficiency for a uniformly doped base layer. 

A series of calculations have been made to determine the optimum drift 
field conditions for highest effieicny for the three different lifetime 

cases. The general impurity profile considered is shown in Figure 4. The 

+ 19 3 

n- surface layer is taken as uniformly doped at 2x10 /cm . The back surface’ 

p"^ layer is also uniformly doped at 10^^/cm^. The base region was assumed 

to have a constantly doped region near the back p"*" region of varying doping 

density N^. The graded doping region was varied in width and type of doping 

profile to study the effects on overall efficiency. The types of doping 

profiles studied were 1) constant, 2) exponential and 3) erfc. The erfc 

profile gives the largest electric field at the p-n junction while the 

exponential profile gives a uniform electric field. 
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TABLE III. COMPARISON OF TRANSITION VOLTAGES V* AND V* WITH THE 

OPEN CIRCUIT VOLTAGE V AND THE MAXIMUM POWER VOLTAGE V 

OG m 


(V/cm) 


(volts ) 

V* 

H 

(volts) 

(volts ) 

V 

oc 

~'Cvolts) 

0 

Max 

0.334 

0.780 

0.600 

0.690 

10 

Max 

0.286 

0.572 

0.597 

0.693 

20 

Max 

0.148 

0.309 

0.594 

0.698 

0 

Med 

0.428 

0.805 

0,450 

0.623 

10 

Med 

0.288 

0.577 

0.444 

0.561 

20 

Med 

0.149 

0.315 

0.413 

0.524 

0 

Min 

0.448 

0.893 

0.309 

0.613 

10 

Min 

0.346 

0.685 

0.427 

0.501 

20 

Min 

0.165 

0.422 

0^301 

0.385 


Max L 


'D 


Eb=3 

Bb =3 


Min L 


D 


Eb=10 

Eb=10 

E.,=10 


TABLE IV. CALCULATED PARAMETERS FOR CELLS WITH 
DIFFERENT AVERAGE DOPING DENSITIES 


(cm"^) 


3.14x10 


3 . 02x10 


17 

17 


1.37x10 


20 


6.87x10 


18 


1.37x10 


21 


N. 


(cm 


1.15x10 

5.6x10^^ 


15 


5x10 


17 


2.5x10 

5x10^® 


16 


sc 


(mA/cm ) 


42.63 

42.28 


23.28 

26.61 

18.53 


V 

oc 

( volts ) 


0.689 

0.694 


0.617 

0.530 

0.674 


CF 


0.839 

0.844 


0.818 

0.790 

0.824 


(%) 


18.21 

18.29 


8.69 

8.23 

7.60 
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Table V shows calculated results for cells with the maximum values 
for different doping profiles. For these large values of the diffusion 
length is larger than the cell thickness and therefore the drift field 
width was taken as equal to the cell thickness. The field is seen to have 
little effect on cell efficiency giving an increase from 18.19% for a BSF 
cell to only 18.29%. Larger fields than those shown in the table were 

found to give lower efficiencies because of reductions in open circuit 
voltage. 

Calculated results using the medium diffusion length case are shown in 
Figure 5. In this series of calculations the doping density was varied 
exponentially between 9xl0^®/cm^ at the p-n junction to 10^®/cm^ at the 
edge of the drift region (W). The width of the drift region was varied 
from 20 ym to 100 ym with the results shown in Figure 5. The peak 
efficiency is seen to occur when the drift field exists over about 40 ym. 

This peak value of 14.64% is slightly better than the 13.70% obtained without 
the drift field. 

Similar calculations to those shown in Figure 5 have been made with the 
graded region doping varying between 9xl0^®/cm^ and 10^^/cm^. The results 
as a function of width of graded region are similar to those of Figure 5. 

The major differences are a peak efficiency of 14.70% for a graded region 
width of about 80 ym. 





TABLE V. COMPARISON OF CELLS WITH DIFFERENT 
doping PROFILES (W=144.8 )jm, MAX L^') 



(cm ) 

(mA/cm^) 

(volts) 

5.6x10^® 

5.6x10^® 

42.40 

0.690 

3.02x10^'^ 

5.6x10^® 

42.28 

0.694 

1.15x10^^ 

5.6x10^® 

42.20 

0.695 

10^7 

5.6x10^^ 

42.73 

0.692 

CO 

o 

5.6x10^® 

41.72 

0.699 


TABLE VI. MAXIMUM CALCULATED EFFICIENCIES 


E 

(V/cm) 


(mA/cm ■) 


V 

oc 

Cvol-f-R 


Calculations for the medium diffusion length case using an erfc 
doping profile between the same doping limits as used for the exponential 
doping have also been made. The results are very similar to those shown 
in Figure 5. For an erfc doping profile between 9xl0^®/cm^ and 10^®/ cm^ 
the peak efficiency was calculated as 14.62% which is slightly less than 
the 14.64% value obtained with the exponential doping. For the erfc 


doping profile between 9xlO^Vcm^ and lO^^/cm® the peak efficiency was 
found to be 14.76% at a graded region width of about 100 pm. This efficiency 
IS slightly larger than the corresponding value for the exponential profile. 
Calculations made with the minimum diffusion length curve of Figure 2 


have given results similar to those discussed above. The largest calculated 
efficiencies were found to occur when the doping density was graded between 
5xl0^'^/cm‘' at the p-n junction to 10^^/cm® at some distance W from the 
junction. Optimum values of W were found to be about 10 pm for the 
exponential grading and about 5 pm for the erfc grading with peak calculated 
®^^^^^®ncies of 9.50% and 9.45% respectively. 

In all of the calculations, the maximum efficiency has been observed 
to occur when the width of the drift region is approximately twice the 
diffusion length. This appears to be the best compromise between enhanced ' 
collection efficiency and increased dark current. 


The Hiaxinmm efficiencies in various drift field type cells calculated 
ill the present study are shown in Table VI. Shown for comparison purposes 
(in parenthesis) are the maximum efficiencies for cells with uniformly 
doped base layers. These values may not represent absolute maximum values 
since optimization studies were not done on all possible parameters. However 
they should be close to the efficiency enhancement to be expected in drift 
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field cells. The largest improvements are observed in cells with low 
diffusion lengths which is to be expected. However the improvement 
which can be achieved is fairly limited. This is consistent with the 
first order models of Section II which show that the collection depth 
can only be improved by about a factor of three. The major limitation, 
however, is the increased dark current and lowered open circuit voltage 
which tends to accompany any attempt to build in a large drift field 
through the use of doping gradients. 

IV. RADIATION RESISTANCE 

One of the major reasons for considering cells with built-in fields 
is the potential for improved radiation resistance. This is consistent 
with the results of the previous section where it was shown that the drift 
field had the largest Improvement in cells with low diffusion lengths. 

The degradation in lifetime with radiation fluence (p is normally 
modeled by an equation of the form 



o 


where is the initial lifetime and K is a damage coefficient, which may 
change with doping density and radiation dose. In the present calculations 
for electron irradiation, K has been taken as a constant and of value 
3.2x10 ^ cm^/sec. 

To Investigate the expected degradation in efficiency, calculations 

have been made on three different cells with the basic device parameters 

of Table I. In the first cell the base region was uniformly doped at 
16 3 

5.6x10 /cm . In the second cell an exponentially graded region from 
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5.6x10 /cm to 10 /cm was taken over a distance of 54 jam, while in the 
third cell the graded region width was taken over a distance of 18 ym. 

The calculated changes in maximum efficiency with electron fluence are 
shown in Figure 6. It is seen that the drift fields do improve the 
efficiency for large radiation doses. The graded width of 54 ym and 
18 ym.were selected as approximately twice the diffusion length after 
total doses of 10^^ and 10^®/cm^ respectively. The cell with a 54 ym 

graded region is seen to be best at '10 /cm and the’ 18 ym cell is best 

16 2 

at 10 • /cm as expected from the selection of graded region widths . 

Other potential devices have also been studied and in all cases an 
enhanced radiation resistance could only be achieved at the expense of 
the initial cell efficiency. 

V. CONCLUSIONS AND SUMMARY 

In this work a detailed numerical computer analysis program has been 
used to study the efficiency of silicon solar cells with various internal 
drift fields achieved by use of doping gradients. In agreement with 
earlier works, it has been found that a drift field can be used to 
significantly enhance the short circuit current of solar cells with short 
diffusion lengths.' However, these improvements are to a large extent offset 
in terms of peak efficiency by increases in dark current and by reductions in 
curve factor when a drift field is present. 

For a given ratio of doping density across the drift field region, 
there is an optimum width for the drift field region. The optimum width 
was found to be on the order of twice the diffusion length when the 
diffusion length is less than the cell thickness. When the diffusion 
length is larger than the cell thickness, very little improvement was 
found in efficiency due to the drift field. This conclusion holds only 


EFFICIENCY (%) SHORT CIRCUIT CURRENT 

( mA /Cm2 ) 





for cells with a high-low junction at the back surface and larger 
improvements are observed for n^-p cells with a back surface ohmic 
contact. 

The calculations indicate that drift field cells can have somewhat 
higher efficiencies after electron irradiation than similar cells with- 
out the drift field. However, this is obtained only at the expense of 
lower initial efficiencies. 
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Appendix 10.4 Two-Dimensional Model of a Solar Cell 
10.4,1 Introducric^ 

A solar cell is in general at least a two-dimensional device. One 
dimensional is parallel to the p-n junction where the light-generated cur- 
rent flowsjwhile the other dimension is perpendicular to the p-n junction. 
Although the optical current may be uniformly generated over the junction 
area, the surface current density is non-uniformly distributed over the 
active area. The sheet resistance and contact resistance which are tra- 
versed by the surface current density are thus functions of the relative . 
locatxon of a given area to the finger contact. All these considerations 
add to the complexity of a solar cell and require the use of a two-dimen- 

i’ 

sional model for an accurate analysis of the terminal properties of a 
solar cell. 

The equations of the two-dimensional carrier flow and its associated 
photovoltaic potential have been previously developed in analytical form 
[1-3], Since the equations are non-linear functions of the current density 
and s,eries resistance, they can only be solved in closed form under very 
stringent assumptions such as low light levels or specimens with special 

• 7 

contact shapes. Hence the use of closed form equations is very limited. 

In this chapter a general two-dimensional model for solar cell analysis 
is proposed and developed. A comparison of the calculated results to 
experimental data is also presented. 

10.4.2 Distributed Resistance and Current Density Model of a Solar Cell 

The dxstributed resistance and current density model considered here 
IS actually a two-dimensional array of ideal one -dimensional solar cells 
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interconnected by series resistance. Shown in Figure lO.lis an array of 
NX and NY points between two grid fingers which are used to locate an 
element of the two-dimensional array. The number of grid points in the 
X-direction is NX and the corresponding number in the Y-direction is 
NY. Figure 10.2 shows the equivalent circuit of the two-dimensional model. 
Figure 10. 3’ shows the distributive resistance and current model at each grid 
point where the current source is that of an ideal one- dimensional model 

including both optically generated current and forward bias current. 

1 / 

for the present work, this current is calculated from the tabulated current 
vs. voltage values obtained from the one-dimensionaJ computer analysis. R is 

O i 

the sheet resistance on the . surface and i.s the contact .resistance between 
the metal-semiconductor interface. The collecting metal is also assumed 
to have a finite resistance which is included between array elements 
located on the edges of the grid fingers . 

The distributive resistance elements can be calculated from the 
following equations . " 


a. 

Sheet resistance - 

Vi ° 

^+1 " ^ST AY’ 

f (lai) 

/ 

(1Q2) 

b. 

Sheet resistance 


^+NX "" ^ST AX’ 

c. 

Contact resistance 

R’ = 
c 


(10.3) 

d. 

Metal resistance 

M M 

T 

Ac ’ 

(10.4) 


where AX and AY are the spacings between grid points in the X- and Y- 

direction respectively, R^^ is the sheet resistance in Q./Q, R^ is the 

.... 2 

contact resistivity in S^*cm , and R^^ is the metal resistivity in J^*cm. 
is the incremental area of the metal contact and T is the thickness of 


the metal contact grid. 
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Figure 10.1. Array of Grid Points used for Two-Dimensional Solar 
Cell Calculation. 
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Any bulk resistance in Figure 10,2 which ai''ises from the bulk 
resistivity of the base layer need not be included hei'e, since it is 
already included in calculating the ona~dimensional I-V charaatoristics . 
Also the distributive base resistance in Figure 10.2 can be usually 
neglected, since most solar calls are covered with an ohmic contact 
over the entire back surface . 

The bus bar is assumed to be at a constant potential which equals 

the terminal solar cell voltage. The voltage of all other grid points 

then can be calculated from the simultaneous solution of the node 

voltage equations at each grid point. 

The node voltage equation can be written fox' each point of the array 

in the active area (Equation 0.0.5)), for points under the grid contact 

(Equation CL0.S))and for the points on the grid pattern (Equation ( 10.7) ) 

F[V(N)>V(N) [GCN~.l.)+G(Ntl)+G(N+NX)+G(N-NX)>V(N-l)*G(k-l) 

-V ( N+1) • G(Ntl) - V(N-NX )‘ G(M-NX ) - V( N+NX )C (N+NX>-Ij^C V (N) >0 (10.5) 

FCV(M)]=V(M)CG(M*-l)+G(M-NX)/2tG(MtNX)/2+GC(M)]-V(M-l).G(M-l) 

-V(M-NX)*G(M-NX)/2-V(MtNX)‘G(M+RX)/2-GC(M)*V(p)-I^jCV(M)3=0, (10.6) 

FCV(P)D=V(F'!)*|l3C(M)tGF-V(P‘-l)+GF •V(P+l)]-^V(P~l) •GF-V(Ptl) • 

GF~GG(M).V(M)“0. (10.7) 

The task now becomes a problem of solving a system of (NX+2)(NY+1) simul- 
taneous equations in the same number of unknotnis. 
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10.4.3 Numerical Algorithm 

The computer algorithm which has been used to solve this system of 
simultaneous equations is an iterative solution technique based on tile 
modified Newton-Raphson method. A brief discussion of this technique 
follows . 

The Newton-Raphson method can be derived from a Taylor series 
expansion. For a single function f(x) = 0, the algorithm used is 
= x^-f(x^)/f ’ (x^) where is the approximate value at the 

(i+l)th iteration. For two coupled equations with two unknowns f^(x,y)=0 
and f 2 (x,y)= 0 , the algorithm can be written as 


^i+1 


1 

X. - -r- 

1 J 




9 y 




^1+1 = - j 


where J is the Jacobian 


J = 


3 X 


3x 




, f^Cx^.y^) 


3 f (x. ,y . 3f (x. ,y, ) 
1 I’-'i 1 


3x 


9y 


Sf 2 (x.,y.) 3fj(x.,y.) 


9 X 


9y 


For a solution to occur , J must not be zero. 


( 10 . 8 ) 


(10.9) 


( 10 . 10 ) 
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The requirement of a non-zero Jacobian is difficult to check before 
running the program and the computation of the matrix operations is time 
consuming, hence a simpler modifier Newton-Raphson method is generally 
used. This consists of applying the single-variable Newton-Raphson method 
n times, once for each variable in a system of n simultaneous equations. 
Each time we do this, we assume that the other variables are kept constant. 

Considef^ as an example two equations with two unknowns such as 


fj^(x,y) = 0, 


( 10 . 1 * 1 ) 


Taking x^ and y^ as the initial guesses new values are obtained as 


( 10 . 12 ) 


0 


X, = x - 




1 


(10.13) 


9x 





(10.14) 


9y 

The algorithm is then repeated until the desired degree of accuracy is 
achieved. An important question is which variables should be used to calculate 
the next approximate solution and in what order. 

One simple example given below will illustrate this point. 


f^(x,y) = 0.2X t 0.1 y-0.01, (10.15) 

f^(x,y) = O.IX + 0.2Y-0.01. 


(10.16) 
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When we use fj^(x,y) to calculate x and f^Cxjy) for y, convergence to an 
-4 . 

error of 10 is achieved in 14 iterations. While the choice of f2(x,y) 
to calculate x and f^(x,y) to calculate y gives a fast divergence. In 
general, it can be easily shown that the function with the steeper slope 
at the solution point with respect to variable x should be chosen to 
calculate the next approximate x, and similarly for y. 

The question of convergence for the modified Newton-Raphson method 
is a touchy one, since one cannot always guarantee a solution. For n 
simultaneous equations with n unknowns, there are n! ways of picking the 
variables and order of execution and sometimes, only one of these choices 
will converge [4]. 

Sometimes the modified Newton-Raphson method does not converge but 
instead oscillates back and forth around the solution. This raises the 
question of when to stop the iterations. In this work the iteration has 
typically been continued until the maximum changes of the variable is 
belov7 some selected small value, but what this may mean with respect to 
the answer is another question. In fact, it is possible that the differ- 
ence between two successive calculations may be very small even though 
the values are nowhere near the right answer in the case of very slow 
convergence. To overcome this difficulty, the so-called under-and-over 
relaxation method has been used as a weighting parameter in the variable 
correction equation. 

In this work, the unknown variables are always chosen from the equa- 
tion which shows the steepest slope with respect to that particular 
variable at the solution point. The order of evaluations are 
arranged to assure fast convergence with an appropriate relaxation weight. 
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Suppose V(N)^ is the j'th iteration value of V(N). The next 
correoted value by using the one-variable Newton-Raphson method is then 


V(N) = V(K). - 


FCVCN)j] 

F^[V(N) J’ 
D 


where F'CV(N)^] is the partial derivative with respect to V(N) . . 
explicit form of F'[V(N) .] is ^ 


(10.17) 


The 


F’[V(N)^]=G(N+L)+G(N-1)+G(N-NX)+G(N+NX) - (10.18) 

for the grid points on the -active area. Ij^CV(N)] is the current density 
of the dark current density superimposed on the optically-generated 
current density. The current derivative can be accurately calculated, if 
the injected dark current density is assumed to be an exponential function 
of the potential at each particular grid point. 

A flow chart of the two-dimensional analysis program is shown in 
FigurelO.4. The input data consists of the one-dimensional I-V data plus 
structural data on the contact finger arrangement. The complete two- 
dimensional I-V characteristic of the solar cell is then calculated at 
specified terminal voltage points using the modified Newton Raphson method. 
Calculations have typically been made until the voltage is accurate at each 
array point to loss than lO-** volts . -With the voltage known at each array 
point the total solar cell current can then be evaluated ]§ summing the 
current contributions from each node in the array. 



Figure 10.4. Flow Chart of Two-Dimensional Analysis Program 
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10.4,4 Discussions 

Several questions may arise concerning the accuracy of this model . 

The first question concerns” the number of grid points used in the cal- 
culation. Theoretically we can use as many points as we want to improve 
the precision in calculating. In this work, an array of 20x20 grid points 
has generally been used unless mentioned otherwise. This results in errors 
of only a few percent even in very severe conditions such as high sheet 
and contact resistance. 

The second question concerns the validity of the one dimensional I-V 

characteristxcs used in the two-dimensional analysis. Since the typical 
built-in field in the diffused surface is in the range of 10®-10^ volts/ 

cm,v7hich is several orders of magnitude larger than the transverse field, 

the development of a transverse field in the two-dimensional analysis has 

a negligible influence upon the minority current density obtained from 

the one- dimensional model. 

The third question is the assumption that the optically generated 
carriers in the base region will only be collected at the junction under 
the illuminated area. It is possible for the generated carrier in the 
base region to diffuse to the junction under the grid pattern, especially 
for the case of a thin collecting metal grid and a long diffusion length ■ 
of the base minority carrier. For "good" cells with low sheet and contact 
resistance , this effect is found to be relatively unimportant, since the 
total collected optical current is the same in both cases. 

The fourth assumption concerns modeling of contact ’resistance 
through a distributed surface contact resistor. In good solar cells, the 
surface contact resistance is usually very small because of the high 
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surface doping density. The base contact resistance is limited by the 
substrate doping density for non BSF cells and may not be negligible. 

The base contact resistance can however be transformed into an effective 
surface contact resistance. The calculated photovoltaic potential is 
the potential difference across the p-n junction irrespective of the 

combination of contact resistance from the surface or base. The relative 
voltage drop across the surface or base coiitact resistance can be deter- 
mined experimentally. 

i 

10.4.5 Conclusions 

(A) A general two-dimensional program has been developed which can simu- 
late a practical solar cell with any arbitrary grid pattern and series 
resistance. 

(B) A general two-dimensional program is a good tool for the optimum 
design of grid patterns and the prediction of the non-linear series 
resistance effects at the maximum power output. 
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Appendix 10,5 Fill Factor and Diode Factor of Solar Cells. 


10.5,1 Abstract 

The interrelation between fill factor, diode factoid series resistance, 
saturation current and space charge current densities for solar ■ 
cells has not beenderlniteljr identified before. In the literature, 
a simplified dark I-V characteristic is usually assumed as follows: 



AkT 


where and A are taken as two independent parameters. Howeyer, this 
assumption cannot be physically justified, and the extensive adoption 
of this equation may sometimes lead to erroneous conclusions. In this 
paper, a bettef physical model of the cell current-voltage characteristics 
is used and the fill factor dependence on the diode factor, series 
resistance, saturation current and space charge current densities can 
therefore be more accurately predicted. 


10.5.2 Introduction 

In recent years, considerable effort has been made to raise the 
conversion efficiency of photovoltaic devices. The efficiency of a 
solar cell can be conveniently represented by 


n 


V I CFF-V *I 

M M ^ '^OC SC 


IN 


■IN 


Cl) 


where the fill factor CFF is a measure of the sharpness of the photo- 
voltaic I-V characteristics of a solar cell. It is clear that high cell 
efficiency can be achieved by raising the values of CFF, V and I 

oc sc 
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However, it is also clearly recognized that increasing one component 
of the above parameters will adversely affect other components. The 
interelationships between fill factor and other cell character- 

istics have not always been definitely identified. 

In previous analyses, there are some discrepancies about the 
dependence of CFF on the above cell characteristics. Lindmayer II] 
suggests that the reduced fill factor of a practical solar cell is due 

( ' - V I 

to space charge recombination Current. However the effects of aeries 
resistance and saturation current density are not considered in this 
analysis. Pulfrey [2], on the other hand, shows that the fill factor 
of a solar cell is principally determined at a given series resistance 
by the saturation dark current density instead of the diode A factor. 

Ihe fill factor Is also found to. ln.provo with an increasing value of 
diode A factor for a -constant series resistance. 

On the contrary. Hovel {3J and Green f4J point out that, at a given 
series resistance, the fill factor is reduced with a higher value of 
diode A factor. However, it is very doubtful that a complete independence 
between the saturation current density or open circuit voltage and the 
diode A factor exists in a practical solar cell as was arbitrarily 
assumed in the previous analyses [1-4J . 

This paper will calculate the fill factor of a solar cell in a 
more general analysis without the constraints and assumptions of the 
prevxous references [1-4J . The dependence of fill factor on the cell 
characterxstics, namely the saturation current density, the space charge 
current density, the diode A factor, and the series resistance is pre- 
sented and examined. 
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10.5.3 Calculation of Fill Factor and Diode Factor 

It has long been recognized that silicon p~n junction behavior 
does not follow ideal diffusion theory. In addition to injection 
into the neutral region of a p-n junction, electrons and holes can 
recombine via the localized recombination centers in the space charge 
region without surmounting the potential barrier. In the quantitative 
treatment of reference [5] , which assumes a Shockley-^ead recombination 

ir 

center, it is shown that 




where the first term represents normal injection; and the second term 

is the space charge current,. The series resistance R can be represented 

o 

by the following equation in the first order model; 

R ^ /CB 

S 12-AA 3-T.WF-AA AR (AR-AA) (3) 

where the meanings of the above parameters are explained in Table 1. 


It is also shovm that the n value of Equation (2) varies between 1 and 
2, where the value of n is 2 if the recombination peaks in the vicinity 
of the center of the space charge region [6]. 

If the superposition of the light and dark currents of an 
illuminated solar cell is assumed to be valid, then the photovoltaic 
I-V characteristic becomes 

^LT “ ^DK “ ^SC ’ 
or 

^LT “ Io2-(^PtqCV-Rg-Ij,^)/2kTj-^l)-I^^,C5) 



213 


Where „-2 has been aas„ed i„ the space charge region Irent. The 

Short circuit current and the open circuit voltage can be inter-, 
related by setting I^^=o in Equation ( 5 ), 

^SC ^01*^®^P^qVQ^/kTj-l) + i^^.(expXqV^^/2kxj-l) . 

The .axi^u. peer can be calculated by setting d(I^^,v)/dV to cero add 

calculating the corresponding and I„. The following equation is 
found: 




q(Vjj-Rg.lj^j,)/2kI. The fill factor CFP can be calculated as 


( 7 ) 


CFF = 


V I 
MM 

V I 
OC SC 


(. 8 ) 


instead of using Equation ( 2 ) to describe the dark l^V characteristics 

01 a solar cell, the following simple equation is frequently used in the 
solar cell literature [1-4] 

^DK “ ), 

Where A is a dimensionless number which is usually found equal to or 

greater than 1 in a practical silicon solar cell. Also A is not a 

constant value, but changes with the terminal voltage or current level 

of the device. In addition to the uncertain value of A, has no 

physical meaning ercept as a parameter fitting to Equation C9) . This 

value Of 1^ does not correspond to the actual reverse saturation current 

density, and it is also different from the prediction of simple diffusion 
theory. 
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From Equation (.2) and C9)i the empirical A factor can he calculated 


at th.e maximum power point as: 


kT 

q 


■ ’ 

logj AfT 


CIQ) 


where and are two voltages in the vibinlty 

10.5.4 Relationships Between Fill Factor, Saturation Current Density, 
Series Resistance and A Factor 

Hnmerical calculations of Equations Cl) to (10) have been performed 

for cells of Tables 1 and 2. The results are shown as Figures 1 to 3. 

Figure 1-a shows that the fill factor decreases with the series 

resistance and the space charge current density. The diode A factor 

can be calculated from Equation (10) for a specific value of the series 

resistance and the space charge current density. The interrelation 

between fill factor and the corresponding diode A factor is thereupon shown 

as Figure 1-b. This figure clearly shows that the fill factor is a 

monotonically decreasing function of the diode A factor for a constant 

series resistance. Moreover, the dashed line shows that the fill factor 

will decrease with the series resistance at a given value of the space 

charge current density. At a zero series resistance, the fill factor is 

-3-4 

about 0.80 or 0.75 at a current ratio 

factor is about 1.08 or 1.58 In this case. Howeyer, with the series 

resistance of model c in Table 2, the fill factor has a value of 

0.77 or 0.72 and the corresponding A factor is 1.2 or 1.92 at a 
i -3 -4 

current ratio respectively. 



Table 1 


Area 

Active Area 
Sheet Resistance 

No. of Fingers 

Width of Fingers 

Space of Fingers 

Length of Fingers 

Thickness of Fingers 

Length of Bus 

Width of Bus 

Base Bulk Resistance 

Surface Contact Resistance 

Base Contact Resistance 

Saturation Current Density 

Space Charge Current Density 

Short Circuit Current 


.Baseline Celig 
AR 


(cm^) 
AA (cm^) 

^ST 

NF 


WF 

SF 

LF 

T 

LB 


(cm) 

(cm) 

(cm) 

(cm) 

(cm) 


_ ion* cm 

in* cm 

0 * pT 

2 

2 

— — — — 

2 

1.734 

■1.734' . 

1.734 

120 

120 

120 

10 

10 

10 

0.015 

0.015 

0.015 

0.00019 

0.00019 

0.00019 

0.9375 

0.9375 

0.9375 

0.0635 

0.0635 

0.0635 
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WB 

(cm) 

0.0625 

0.0625 

RB 

(n* cra^) 

0.25 

0.025 

‘‘cs 

(n* cm^) 

lO""^ 

lO"'^ 

\b 

(n* cm^) 

Varied 

Varied 

^01 

(A/ cm^) 

6. 94xl0~^^ 

6 . 83x10 

^02 

(A/ cm^) 

Varied 

Varied 

'sc 

(mA/ cm^) 

43.0 

40.0 


-12 


2 

0.0625 
0.0025 
10""^ 
Varied 
1 • 76x10 
Varied 
35.0 


-13 


Table 2. 


Models of the Spt-tqo d 

■factor of a solar ^ calculation of Fill 


Models 

a 

b 

c 

d 

e 

f 


Rqb cm^) 


0 

10~^ 

0.5 

1,5 

3.0 

5.0 


10 cm 

Rs(fl) 

0 

0.32 

0.57 

1.07 

1.82 

2.82 


1 fl* cm 

Rg(«) 

0 

0.21 

0.46 

0.96 

1.71 

2.71 


0.1 fi • cm 

%(^) 

0 

0.20 

0.45 

0.95 

1.70 

2.70 


FILL FACTOR 



LO 6,^ 101/102) 

Figure 1-a. Fill factor dependence on the space charge current density and series resistance a to f 
for 10 n-cm cell. 
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FILL FACTOR 


0.9 


DIODE FACTOR 

Figure 1-b. Fill factor dependence on the diode A factor and series resistance a to f for 10 fl-cm 
cell. logj(jQ^/lQ 2 )=- (-• ), -4(™). 


FILL FACTOR 



8T5 


FILL FACTOR 




FILL FACTOI 



Figure 3-a. 


current density and series resistance 

“ to X tor u.l s2— cm cell. 




FILL FACTOR 



222 




It is known that a conventional n"^p silicon cell has a fill factor 
of about 0.72. On the other hand, the theoretically predicted fill 
factor is about 0.81. This discrepancy has been identified above, as a 
result of a high space charge current density and a high series 
resistance in a practical silicon solar cell. Similar results are shown 
in Figures2-a, 2— b, 3— a and 3—b for a 1 fl*cra or 0.1 0* cm solar cell. From 
the comparison o( parts a and b of Figures 1, 2 and 3, it is clearly seen 
that the fill factor is higher for a low resistivity solar cell at a 
specific series resistance and space charge current density. Therefore 
a low resistivity cell has the advantage of a higher fill factor than a 
high resistivity cell, 

10.5.5 Discussions and Conclusions 

^kis section discusses the discrepancies in previous theoretical 
analyses [1-4] . As has been shown, the previous analyses of fill factor 
began with the simple Equation (9) . In this representation, the values 
of the two empirical parameters of A and 1^ need to be varied along with 
the voltage and current levels in order to fit the dark I-V characteris- 
tics of a practical solar cell. Therefore, it is doubtful that one can 
arbitrarily assume independent values of A and I^. However, 

it is still a good representation if it is used to calculate the cell 
characteristics at or around only pne particular point such as the 
maximum power density. Moreover the values of A, 1^ and Rg are closely 
interrelated and they can not he treated as independent parameters. 
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In this paper, Equation (2) instead of Equation (9) has been used 
in the analysis. There is a definite physical meaning for each 
parameteir of Equation (2). In our calculation of the fill factor, the 
space charge current density i^ allowed to vary due to the physical 
condition in the space charge region. The diode A factor is calculated 
at the maximum power point; therefore, the fill factor and diode factor 
can be correlated without the limitation of the assumptions of the 
previous analyses . 

It can therefore be concluded that the fill factor of a solar cell 
is a monotonically decreasing function of the series resistance, the 
space charge current density and the diode A factor. In order to 
obtain a high value of fill factor, it is important to design a solar 
cell with a low series resistance and space charge current density 

A low resistivity cell is also found to give a higher fill factor than 
a high resistivity solar cell. 
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